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Introduction 



Superconductivity 

The mysterious effect of superconductivity (zero electric resistance) in a metal at low 
temperatures was discovered by H. Kamerlingh Onnes in 1911 Q]. For a long time, the 
microscopic explanation of this phenomenon was lacking, while a number of successful 
phenomenological approaches (the London theory 0, the Ginzburg-Landau theory 
jH]) were employed to study various properties of superconductors. 

During this time, the effect of superfluidity (flowing without dissipation) in the 
liquid helium *He was experimentally discovered by P. L. Kapitza ^ and theoretically 
explained by L. D. Landau 0- The basic mechanism behind the superfluidity is the 
Bose-Einstein condensation of helium atoms. Although at first sight, Hquid helium 
does not have much common with solid superconductor, actually the phenomena of 
superfluidity and superconductivity are very close to each other. While the superflu- 
idity is the flowing of helium atoms without dissipation, the superconductivity is the 
flowing of electrons without dissipation. Thus superconductivity can be imagined as 
the superfluidity of electronic Hquid. At the same time, there is a crucial difference 
between the two phenomena. While a macroscopic number of bosonic heHum atoms 
undergo the Bose-Einstein condensation and occupy the lowest energy level, the elec- 
trons, which are fermionic particles, cannot do this as they are not allowed to share 
the same state. 

A key step in unravelling the mystery of superconductivity was made by 
L.N. Cooper in 1956 jQ, who showed that if electrons attract then the electronic 
system becomes unstable against forming the so-called Cooper pairs (the attraction 
that can overcome the Coulomb repulsion can be produced due to interaction of 
electrons with lattice vibrations). The newly formed pairs of electrons are bosonic 
objects, hence they can undergo the Bose-Einstein condensation and form a state 
necessary for the superfluidity! This idea became the keystone of the microscopic 
theory of superconductivity which was finally built in the late fifties by J. Bardeen, 
L.N. Cooper, and J.R. Schrieffer (BCS) and N.N. Bogolyubov 0. Shortly there- 
after, L. P. Gor'kov formulated the BCS theory in the language of the Green functions 

m 

Later, the possibility of unconventional superconductivity, which is characterized 
by nonzero spin or orbital momentum of the Cooper pairs, was discussed. The ordi- 
nary case corresponds to Cooper pairing between electrons with opposite spins, the 
singlet pairing. The other, unconventional case, corresponds to pairing between elec- 
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trons with the same spin projection, the triplet pairing. Another nontrivial possibiUty 
is a nonzero orbital momentum. The superconducting state is characterized by the 
order parameter A, which describes the "strength" of the superconductivity and de- 
termines the energy gap in the single-particle density of states. In the ordinary case, 
the energy gap A is the same for all directions of electron's momentum. However, 
if the orbital momentum of the Cooper pairs is nonzero then the superconductivity 
is anisotropic, which means that A depends on the direction of motion inside the 
superconductor (an example of anisotropic superconductivity is established in the 
high-temperature superconductors jlOl Illj V 

Extensive investigations demonstrated that zero electric resistance is only one of 
the many manifestations of superconductivity. Below we mention two of them which 
are particularly relevant for the subject of the thesis and appear in its title. 

When a superconductor contacts a normal metal, a number of phenomena known 
as the proximity effect takes place. The two materials influence each other on a 
spatial scale of the order of the coherence length in the vicinity of the interface. In 
particular, the superconducting correlations between quasiparticles are induced into 
the normal metal, because the Cooper pairs penetrating into the normal metal have a 
finite lifetime there. Until they decay into two independent electrons, they preserve 
the superconducting properties. Alternatively, the proximity effect can be viewed as 
resulting from the fundamental process known as the Andreev reflection . Imagine 
a low-energy electron impinging from the normal metal onto the interface with the 
superconductor. A single electron can penetrate the superconductor only if its energy 
is larger than the superconducting energy gap A, while below it only Cooper pairs can 
exist. Thus the low-energy electron cannot penetrate the superconductor. If there is 
no potential barrier at the interface, the electron cannot be reflected back either, due 
to the momentum conservation law. The only way out of this contradiction is the 
Andreev reflection: a Cooper pair goes into the superconductor, while a hole ("the 
antiparticle" for the electron) goes back into the normal metal in order to conserve 
charge. 

Deep understanding of the nature of the superconducting state lead also to the 
discovery of one of the most spectacular superconducting phenomena, the Joseph- 
son effect 53]. This occurs when two superconductors are connected via a weak 
link, which can be either a nonsuperconducting material (insulator as in Ref. ^3] or 
metal) or a geometrical constriction (the variety of Josephson junctions is reviewed in 
Ref. |1.4,). The superconducting condensate in each of the two weakly coupled super- 
conductors is described by its wave function and the corresponding phase. Josephson 
found that if the phases are different and the difference is (/?, then the supercurrent 
(nondissipative current) I = Ic sin tp arises across the junction in the absence of volt- 
age. The quantity Ic is called the critical current. 



Motivation 

From the most fundamental aspects of superconductivity let us now switch to the 
particular issues that became the subject of the thesis. 
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Although the investigation of the proximity effect in SN systems^ was started 
about forty years ago JllElEI; the technology allowing to produce and measure 
experimental samples of mesoscopic dimensions was achieved relatively recently. In 
particular, it became possible to study SN structures consisting of thin layers (having 
thickness smaller than the coherence length). Such structures behave as a single 
superconductor with nontrivial properties. The most basic of them has already been 
studied mostly for the case of ideally transparent interface. At the same time, the 
experimental progress requires the corresponding advances in theory, especially taking 
into account arbitrary interface transparency. This crucial parameter determines the 
strength of the proximity effect and at the same time is not directly measurable. From 
practical point of view, the SN proximity structures can be used as superconductors 
with relatively easily adjustable parameters, in particular, the energy gap and the 
critical temperature. The parameters of the proximity structures can be tuned, e.g., 
by varying the thicknesses of the layers. This method has already found its application 
in superconducting transition edge bolometers and photon detectors for astrophysics 
(see, e.g., Refs. 

The physics of SF systems is even richer. In contrast to the SN case, the su- 
perconducting order parameter does not simply decay into the nonsuperconducting 
metal but also oscillates. This behavior is due to the exchange field in the ferromag- 
net that acts as a potential of different signs for two electrons in a Cooper pair and 
leads to a finite momentum of the pair (similarly to the Larkin-Ovchinnikov-Fulde- 
Ferrell state in bulk materials [SD!.). This oscillations reveal itself in nonmonotonic 
dependence of the critical temperature Tc of SF systems as a function of the F layers 
thickness, both in the cases of SF multilayers [HIES] and bilayers j22l EHl El] . At the 
same time, in most of the papers investigating this effect, the methods to calculate 
Tc were approximate. An exact method was proposed in Ref. (22] for the limiting 
case of perfect interfaces and large exchange energies. An exact method to calculate 
Tc at arbitrary parameters of the system was lacking. The need for such a method 
was also motivated by the experiment 12^] that did not correspond to the previously 
considered approximations and Hmiting cases. 

Another interesting effect in SF systems takes place if the magnetization of the 
ferromagnet is inhomogeneous. Then the triplet superconducting component can 
arise in the system ^^Ej- Physically, this effect is similar to generating the triplet 
correlations in magnetic superconductors jJHI (with Cooper interaction only in the 
singlet channel). Recently, it was demonstrated that the triplet component also 
arises in the case of several homogeneous but noncollinearly oriented ferromagnets 
|29j . However, the conditions at which the superconductivity is not destroyed in this 
system were not found. The simplest system of the above type is an FSF trilayer. The 
answer to the question about the conditions for the superconductivity to exist can be 
obtained when studying the critical temperature Tc of the system. At the same time, 
a method to calculate Tc in a situation when the triplet component is generated, was 
lacking. A possible practical application of FSF structures is a spin valve |^^()[ \'M\ . a 
system that switches between superconducting and nonsuperconducting states when 
the relative orientation of the magnetizations is varied. Although the superconductive 



^The notations: S — s-wave superconductor, D — d-wave superconductor, N — normal metal, F 
— ferromagnetic metal, I — insulator, c — constriction. 
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spin valve is not yet experimentally realized, the work in this direction has already 
started jS2l- 

The Josephson effect in the systems containing ferromagnets (e.g., structures of 
the SFIFS, SIFIS or SFcFS type) also has a number of pecuHarities. Among them are 
the transition from the ordinary (0-state) to the so-called 7r-state (in other words, the 
inversion of the critical current sign or the additional tt phase shift in the Josephson 
relation — this effect was theoretically predicted in Refs. I33t l34 j and experimentally 
observed in Ref. [SlI) when the two ferromagnets are ahgned in parallel, the enhance- 
ment of the critical current by the exchange field in the case of antiparallel orientation 
[36 , and nonsinusoidal current-phase relation. Although such effects have been stud- 
ied before, it was often done in the simplest models and at the simplest assumptions 
about system parameters. To achieve better understanding of this phenomena, one 
should study them at various conditions and determine the physical mechanisms be- 
hind the effects. The interest to SFS junctions with nontrivial current-phase relation 
is in particular due to their possible employment for engineering logic circuits of novel 
types, both classical and quantum bits (see, e.g., Refs. \Sll l3Bj). 

Another interesting type of nonuniform superconducting systems is a junction 
between superconductors of nontrivial symmetry. The superconductors with the d- 
wave symmetry of the order parameter are widely discussed because this symmetry 
is realized in the high-temperature superconductors |1()[ 111) . A possibility to imple- 
ment a so-called qubit (quantum bit) based on d-wave junctions was proposed in 
Refs. I^iiilji^l. Quantum bit is, simply speaking, a quantum mechanical system 
with two states (which can be imagined as spin 1/2). While a classical bit can be 
either in one state or in the other, a qubit can also be in a superposition of the two 
states. If a quantum computer is built of such qubits, it would have the advantage 
of natural computational parallelism that can enormously speed up certain types of 
computational tasks. A possibility to implement a qubit based on DID junction stems 
from the fact that the energy of such a junction as a function of the phase difference 
can have a double-well form with two minima. This degeneracy of the ground state 
arises due to the nontrivial symmetry of the superconductors. Due to tunneling be- 
tween the wells, the ground state splits, and the two resulting levels effectively form 
the quantum-mechanical two-state system. At the same time, the gapless nature of 
the d-wave superconductors leads to appearance of low-energy quasiparticles which 
can destroy the quantum coherence of the qubit and hence hamper its successful func- 
tioning. Calculation of the corresponding decoherence time is necessary for estimating 
the efficiency of the proposed qubits. 



Outline of the thesis 

In Chapter 1, the theory of superconductivity in thin SN sandwiches (bilayers) in 
the diffusive limit is developed, with particular emphasis on the case of very thin 
superconductive layers, ds ^ dw. The proximity effect in the system is governed 
by the interlayer interface resistance (per channel) pB- The case of relatively low 
resistance (which can still have large absolute values) can be completely studied 
analytically. The theory describing the bilayer in this limit is of BCS type but with 
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the minigap (in the single-particle density of states) <C A substituting the order 
parameter A in the standard BCS relations; the original relations are thus severely 
violated. In the opposite limit of an opaque interface, the behavior of the system 
is in many respects close to the BCS predictions. Over the entire range of Pb, the 
properties of the bilaycr arc found numerically. Finally, it is shown that the results 
obtained for the bilayer also apply to more comphcated structures such as SNS and 
NSN trilayers, SNINS and NSISN systems, and SN superlattices. 

In Chapter 2, we propose two exact methods to calculate the critical temperature 
Tc of dirty SF bilayers at arbitrary parameters of the system. The methods are apphed 
to study the nonmonotonic behavior of the critical temperature versus thickness of 
the F layer. Comparing our results with experimental data, we find good agreement. 
Then we study the critical temperature of FSF trilayers, where the triplet supercon- 
ducting component is generated at noncollinear magnetizations of the F layers. We 
reduce the problem to the form that ahows us to employ the exact numerical methods 
developed earher, and calculate Tc as a function of the trilayer parameters, in par- 
ticular, mutual orientation of magnetizations. Analytically, we consider interesting 
limiting cases. Our results determine conditions which are necessary for existence of 
the odd (in energy) triplet superconductivity in SF multilayers. 

In Chapter 3, the quantitative theory of the Josephson effect in SFS junctions 
(with one or several F layers) is presented in the dirty limit. Fully self-consistent nu- 
merical procedure is employed to solve the Usadel equations at arbitrary values of the 
F-layers thicknesses, exchange energies, and interface parameters. In SFIFS junction, 
at antiparallel ferromagnets' magnetizations the effect of the critical current enhance- 
ment by the exchange field is observed, while in the case of parallel magnetizations 
the junction exhibits the transition to the 7r-state. In the limit of thin F layers, we 
study these peculiarities of the critical current analytically and explain them qualita- 
tively; the scenario of the O-tt transition in our case differs from those studied before. 
The effect of switching between and n states by changing the mutual orientation of 
the F layers is demonstrated. Also, various types of the current-phase relation I{lp) 
in SFcFS point contacts and planar double-barrier SIFTS junctions are studied in the 
limit of thin ferromagnetic interlayers. The physical mechanisms leading to highly 
nontrivial I{ip) dependence are identified by studying the spectral supercurrent den- 
sity. In particular, these mechanisms are responsible for the O-tt transition in SFS 
Josephson junctions. 

In Chapter 4, we study the Josephson junction between two d-wave superconduc- 
tors, which is discussed as an implementation of a qubit. We propose an approach 
that ahows to calculate the decoherence time due to an intrinsic dissipative pro- 
cess: quantum tunneling between the two minima of the double- well potential excites 
nodal quasiparticles which lead to incoherent damping of quantum oscillations. In 
DID junctions of the mirror type, the contribution to the dissipation from the nodal 
quasiparticles is superohmic and becomes small at smah tunnel splitting of the energy 
level in the double-weh potential. For available experimental data, we estimate the 
quality factor. 
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Green functions in the theory of superconductivity 

The quantum field theory methods turned out to be very powerful also in the theory 
of solid state systems ^^I- The theory of superconductivity was formulated in this 
language, the language of the Green functions, by L.P. Gor'kov 0. The Gor'kov 
equations were derived directly from the BCS Hamiltonian. They can be conveniently 
written in the matrix form: 



xG(i;,ri,r2) = i5(ri-r2), (1) 



where the retarded matrix Green function G contains the standard Green function 
G and the anomalous Green functions F and F that describe the superconducting 
correlations: 

-Ef is the Fermi energy, V{v-y) is the impurity potential. 

In real systems, solving the Gor'kov equations can be a formidable task. At the 
same time, the information contained in those is often redundant. In particular, the 
Green functions entering the Gor'kov equation describe rapid oscillations on the scale 
of the Fermi wave-length. Of course, there are situations when these oscillations 
play crucial role and must be carefully described. At the same time, if characteristic 
scales in the problem are large compared to the interatomic distance, then the rapid 
oscillations are averaged and only the slow part of the Green functions determines 
the physical properties. For those cases, the Gor'kov equations can be simplified from 
the very beginning. The resulting formalism is known as the quasiclassical method 
jini^SBll, and below we briefly outline the main steps leading to it. 

The Green functions entering the Gor'kov equations depend on two coordinates 
and energy. We can do the Wigner transform over the coordinates, which means 
that we introduce the center of mass coordinate r = (ri + r2)/2 and the relative 
coordinate p — ri — r2, and make the Fourier transform over p resulting in the 
relative momentum p. The above transformations are exact, they are made simply 
for convenience. The next step is the quasiclassical approximation itself, the idea 
of which is that the dependence of the Green function on the relative coordinate 
is fast, while the dependence on the center-of-mass coordinate r is slow (the small 
parameter of this approximation is Tc/^-f, with being the critical temperature). 
The quasiclassical approximation implies that we average the Green functions over 
the rapid oscillations described by the relative momentum p: 

g{E,r,n)^^ f d^GiE,r,p). (3) 



The integration over p, the absolute value of the relative momentum, is standardly 
rewritten as integration over ^ = p^/2m — Ep. The unit vector n = p/p points in the 
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direction of the relative momentum. Finally, we should average the Green function 
over impurities. 

The equation for the resulting quasiclassical function was derived by G. Eilen- 
berger (32] (see also Ref. (331); it reads 



or 



-i{E + iO) a3 + A + ^,g 
Zt 



0' ^-[a* o')' 



where the square brackets denote the commutator and the angular brackets denote 
averaging over directions of n (angular averaging), vp is the Fermi velocity, txa is 
the third Pauli matrix, and t is the time of the mean free path. The Green function 
entering this equation obeys the normalization condition = 1. 

The Eilenberger equation can be further simplified in the diffusive (dirty) limit. 
The physical reason is that due to frequent scattering on impurities, the Green func- 
tion becomes isotropic. The equation for the isotropic function 

giE,r) ^ {g{E,r,n)) (5) 

was derived by K. D. Usadel |l5] : 



D— (- — - 1 + 



liE + iO) as - A, g =0, g' = 1. (6) 



All forms of the equations discussed above must be complemented with the corre- 
sponding self-consistency equation for the order parameter A. For brevity, we write 
down the self-consistency condition only for the Usadel equation: 



" dEts^nhf^] f{E,r), (7) 



2T 



where / is the anomalous component of the matrix Green function, T is the tempera- 
ture, A is the coupling constant (effective parameter of electron-electron interaction), 
and the integration is cut off at the Debye energy ojd- In particular, this equation 
yields the zero-temperature gap of the bulk superconductor A = 2a;D exp(— 1/A). 

The quasiclassical approximation breaks down near interfaces between different 
materials, when the interfaces are sharp on the scale of the Fermi wave-length (atomic 
scale). In order to use the quasiclassical equations for describing nonuniform systems, 
one should complement them with the effective boundary conditions. The boundary 
conditions for the Eilenberger equation were derived by A.V. Zaitsev and for 
the Usadel equation — by M.Yu. Kupriyanov and V. F. Lukichev [32|. The latter 
boundary conditions have the form 

. dgi . dgr 1 , , , , 

a. = a. g.— = ^[.9^.9.], (8) 

where the indices / and r refer to the left- and right-hand side of the interface, re- 
spectively, a is the conductivity of metals, Rb is the total resistance of the interface, 
and A is its area. The x axis is normal to the interface. 
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Chapter 1 

Superconductivity in thin SN bilayers 



1.1 Introduction 

It is well known that the majority of metallic superconductors is well described by 
the classical BCS theory of superconductivity Jj. One of the main qualitative fea- 
tures of the BCS theory is a simple relation between the superconductive transition 
temperature Tc and the low-temperature value of the energy gap for s-wave super- 
conductors: A(0) = 1.76 Tc. Experimentally, violations of this simple relation are 
considered as a sign of some unusual pairing symmetry or even of a non-BCS pair- 
ing mechanism. Recently, an evident example of such a violation of the BCS theory 
predictions was found in experiments by Kasumov et al. |2] , who studied the current- 
voltage characteristics of a carbon nanotube contact between two metallic bilayers 
(sandwiches) made of ordinary metals, tantalum and gold. The observed value of 
the low-temperature Josephson critical current is 40 times larger than the maximum 
expected (Ambegaokar-Baratoff) value /c = 7rA(0)/2ei?tubo, where the energy 
gap of the bilayer A(0) is estimated from its transition temperature. The source of 
such discrepancy is not clear at present. The most recent experiments jl| demon- 
strate the existence of intrinsic superconductivity in carbon nanotubes. However, the 
discrepancy could also be due to unusual superconductive properties of the bilayers. 

Although the proximity effect in SN bilayers is rather well studied, mostly the 
limit of the perfect interface (which is called "the Cooper limit" in the thin bilayer case 
jHllEI) or the opposite Hmit of opaque interface has been explored. At the same time, 
as we demonstrate below, between the two limiting cases the characteristics of the 
system and the relations between them, depending on the interface resistance, behave 
quite nontrivially. The aim of the present chapter is to study the superconducting 
properties of a thin SN bilayer depending on the interface resistance. 

An essential feature of the experiment [5] was that the superconductive layer in 
the bilayer was very thin: ds/dN = 5 nm/100 nm = 1/20. In the present chapter, we 
investigate such a bilayer both analytically and numerically, calculating quantities 
characterizing the superconductivity in this proximity system: the order parameter 
A, the density of the superconducting electrons n, the critical temperature Tc, the 
(mini)gap Eg in the single-particle density of states (DoS), the critical magnetic field 
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He parallel to the bilayer and the upper critical field Hc2 perpendicular to the bilayer. 
In our calculations, the parameter controUing the strength of the proximity effect is 
the (dimensionless) resistance of the SN interface per channel ps, which is related to 
the total interface resistance i?B by the Sharvin formula 

- (1.) 

where Rq = 7r/e^ is the quantum resistance, and A^ch = -4/(Af/2)^, with Xp being 
the Fermi wave-length, is the number of channels in the interface of area A. We 
choose Xp referring to the S layer. 

The values of the interface resistance can be divided into three ranges: (a) at large 
resistance, many characteristics of the superconductor (A, n, He, Hc2) are almost 
unaffected by the presence of the normal layer — this is the BCS limit (however, 
we note that Eg does not coincide with the order parameter, and even vanishes 
as pb increases); (b) at low resistance, the theory describing the bilayer is of BCS 
type but with the order parameter A substituted by the minigap Eg (for instance. 
Eg ~ 1.76 Tc, whereas Eg <C A); the original BCS relations are thus severely violated; 
(c) at intermediate resistance, the behavior of the system interpolates between the 
above two regimes. 

The violation of the BCS relations, corresponding to the parameters of the ex- 
periment P], appears to be insufficient to explain the observed value of the critical 
current. Probably, the intrinsic superconductivity of the carbon nanotubes 4^ plays 
an essential role. Therefore, the experiment [2] is mainly a motivation for the theo- 
retical study described below. The practical outcome is an estimate for the interface 
resistance, obtained from the experimental values of the critical temperature and the 
parallel critical magnetic field (see Sec. 11.81 belowV 



1.2 Method 

1.2.1 Usadel equation 

Equilibrium properties of dirty systems are described [J] by the quasiclassical retarded 
Green function g(r, E), which is a 2 x 2 matrix in the Nambu space satisfying the nor- 
malization condition = 1. The retarded Green function obeys the Usadel equation 

m 

D^ig^g) + i^Eas - A, g] = 0. (1.2) 

Here the square brackets denote the commutator, D — vl/3 is the diffusion constant 
with V and I being the Fermi velocity and the elastic mean free path, E is the energy, 
whereas as (the Pauli matrix) and A(r) are given by 

The order parameter A(r) must be determined self-consistently from the equation 
A(r) = -y dE tanh (^^^ /(r, E), (1.4) 
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where /, the anomalous Green function, is the upper off-diagonal element of the 
g matrix, A is the effective constant of electron-electron interaction in the S layer 
(whereas we assume A = and hence A = in the N layer), and integration is cut 
off at the Debye energy ujjj of the S material. 

Equation (11.211 should be supplemented with the appropriate boundary conditions 
at an interface, which read [2| 

m {gNngi) = (^r (grV„.gr) = [gu gA , (1-5) 

where the subscripts I and r designate the left and right electrode, respectively; a is 
the conductivity of a metal in the normal state, and gs = Gb/A (with Gb = ^/Rb) 
is the conductance of the interface per unit area when both left and right electrodes 
are in the normal state. V„ denotes the projection of the gradient upon the unit 
vector n normal to the interface. 

We use the system of units in which the Planck constant and the speed of light 
equal unity: h — c — 1. 



1.2.2 Angular parameterization of the Green function 

The normalization condition allows the angular parameterization of the retarded 
Green function: 

. _ / 008 6* e''^ sin 6I\ 

^ ~ l^e~''^sin6' -cose ) ' ^ 

where 6 — 9{r, E) is a complex angle which characterizes the pairing, and — v?(r) 
is the real superconducting phase. The off-diagonal elements of the matrix g describe 
(3 the superconductive correlations, vanishing in the bulk of a normal metal (0 ^ Q). 
The Usadel equation takes the form 



— ^ iy^)^ cos( 



sin6l+ |A|cos6i = 0, (1.7) 

V (sin^ V(/7) = 0. (1.8) 
The corresponding boundary conditions are 

cri'^ndi = gs [cos{ipr — (pi)cos6i sin 0^ — sin6'i cos^r] , (1-9) 

CrVn^r — 5B [cOS 0; SmOr — COs{(pr ~ '■Pi) Sa\9i COS 0^] , (1-10) 

(Ti sin^ 6i Vn'ySi — (Tr sin^ Or Vn^^r — 3b sin{ipr — '■Pi) sn\6i sin6'r. (l-H) 
The self-consistency equation for the order parameter A(r) takes the form 



lAI-A 



j dE tanh^^^ Im [sing] . (1.12) 



The above equations are written in the absence of an external magnetic field. To 
take account of the magnetic field, it is sufficient to substitute the superconducting 
phase gradient in the Usadel equations H1.7|l - lll.8|l by its gauge invariant form 2mv = 
V(/5 -I- 2eA, where A is the vector potential and v denotes the supercurrent velocity. 



22 



Chapter 1 



Physical properties of the system can be expressed in terms of the pairing angle 
6'(r, E). The single-particle density of states v{y, E) and the density of the supercon- 
ducting electrons n(r) are given by 

1/ = 1^0 Re [cos 61], (1.13) 
2m(j 



n = 



dEt&nh(^^^lm[sin'0], (1.14) 

where m and e are the electron's mass and the absolute value of its charge, and 
lyQ = m^v/TT^ is the normal-metal density of states at the Fermi level. The total 
number of single-particle states in a metal is the same in the superconducting and 
normal states, which is expressed by the constraint 

/•OO 

/ dE [i^{r, E) -i^o]=0. (1.15) 
Jo 

Below, it will be sufficient to consider only positive energies, E > 0. 
1.2.3 Simple example: the BCS case 

The simplest illustration for the above technique is the BCS case, when the order 
parameter A(r) = Abcs is spatially constant. Its phase can be set equal to zero, 
ip = 0. Then the Usadel equations (ll.7|l - lll.8|l are trivially solved, and we can write 
the answer in terms of the sine and the cosine of the pairing angle: 



sinOBcsiE) = ^^M^=, (1.16) 

V^' - ^BCS 

cosOBcsiE) ^ ^=£==. (1.17) 



2 

BCS 



An infinitesimal term iO should be added to the energy E to take the retarded nature 
of the Green function g into account, which yields 

Im [sin^ OBcsiE)] = ^Abcs S{E - Abcs)- (1-18) 



The usual BCS relations are straightforwardly obtained from Eqs. (I1.12|l - lll.l4ll 
(for simpHcity, we consider the case of zero temperature): 




at < Abcs 
at E> Abcs 



(1.19) 

(1.20) 
(1.21) 



The critical temperature must be determined from Eq. Ijl.l2ll with vanishing A{Tc); 
the result is 

Abcs (0) = 1.76 Tf^^. (1-22) 
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"T" 




Figure 1.1: SN bilayer. The N and S layers occupy the regions — djv < 2: < and 
< X < ds, respectively. 



1.3 Usadel equations for a thin bilayer 

Let us consider a SN bilayer consisting of a normal metal (— djv < a; < 0) in contact 
(at a; = 0) with a superconductor (0 < x < ds) — see Fig. II. II We choose the phase f 
of the order parameter A equal to zero in the S layer. At the same time, we suppose 
that electron-electron interaction is absent in the normal layer: A = 0, hence A = 0, 
although the superconductive correlations {9 ^ 0) exist in the N layer due to the 
proximity effect. The Usadel equations (ll.7|l - lll.8|l take the form 



Dm 92 



^ ^+tEsmeN = 0, (1.23) 



2 

's 



2 fe2 



i£;sin6ls + Acos6ls = 0, (1.24) 



where 9n and 9$ denote the pairing angle 9 ai x < Q and a; > 0, respectively. The 
boundary conditions H1.9|l - lll.ll|l reduce to 

fjjv-^ = 0-5-^ = gssm(6'5 - t'w)- (1-25) 



We assume that the layers are thin, hence the order parameter A can be regarded 
as constant in the superconductive layer. The standard condition that the layers 
are thin is obtained if we compare their thicknesses to the coherence lengths: <C 
Dm / ^BCS and ds <C y/Ds/ Abcs- However, in the case of low interface resistance, 
the condition can become stronger. The thin-layers assumption will be discussed in 
more details in Sec. 11.81 
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Equations HI .2311 - 111 .2411 can be integrated once, yielding 
Dn fdON 



^ y , iEcos0N = bN, (1.26) 
D / do \ 

~T\~di) -^^^'^^^3 + AsinOs ^bs. (1.27) 

The functions biq^E) and bs{E) are determined from the boundary condition 86 /dx ~ 
at the nontransparent outer surfaces of the bilayer, which give 

bN{E) = —iEcos9Ni—dN,E), 

bsiE) = -iEcosesids,E) + Asmesids,E). (1.28) 

Let us denote On{E) = 9N{—dN, E), 9s{E) ~ 9s{ds, E). Due to small thickness of the 
layers, the functions 9n{x,E) and 6s{x,E) are nearly spatially constant. However, 
in order to determine them, we should take account of their weak spatial dependence 
and make use of the boundary conditions at the SN interface. Substituting 

9n{x,E) ^ 9n{E)+56n{x.E), 

9s{x,E)^9s{E)+d9s{x,E) (1.29) 

intoEqs. Ijl.2(i|l - lll.27|l and linearizing them with respect to \S9n(x. E)\. \69s{x, E)\ <^ 
1, we fmd the solution. Finally, the boundary conditions at the SN interface lead to 

-iTNEsm9N{E) = iTsEsm9s{E) + TsAcoa6s{E) = sm[9s{E) - 9n{E)] , (1.30) 

where we have denoted 

ia^dN 2asds 

TN = 7^ , Ts = — . (1.31) 

Physically, rjv and ts are the escape times from the respective layers (see Ap- 
pendix EJ • Using the definition of the interface resistance per channel lll.l|l , we 
can represent these quantities as 

VAfrfjV ds , 

TN = 2-K — ^Pb, ts = 27r — pB, (1-32) 

with vjv and vg being the Fermi velocities in the N and S layers. The ratio tjv/ts = 
VNdN /^sds, which is independent of the interface properties, can also be interpreted 
as the ratio of the total densities of states (per energy interval) in the two layers: 

Tn AdNl'ON ,^ „„x 

— = —T-, . (1.33) 

Ts AdsvQS 

Having solved the boundary conditions H1.3flll . we can determine all equilibrium 
properties of the system (because knowledge of 6*5, 6* at implies knowledge of the re- 
tarded Green function q). 
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A useful representation of the boundary conditions l|1.3fl|l is obtained as follows. 
Excluding 9n{E) from Eq. (|1.3nil . we arrive at a single equation for the function 
9s{E), which can be written, in terms of Z = e-xjp{i9s), as a polynomial equation 



iCeZ'' + C^Z" + iCiZ" + Cj,Z^ + iC^Z"- + C^Z + iCo = 



(1.34) 



with real coefficients 



C2 

Co 



-tnE 



TS_ 

TN 



2 r 



1 - (tnE) 
-tnE 



H 

TN 



A' 

'^E 



TS 
TN 



C.=2-2 



1 - [tnEY 



TS 
TN 



A 
) 

^ E 



-tnE 



TS 
TN 



1 - [tnEY 



TS 
TN 



-tnE 



TS 
TN 



A' 
E 



A 
) 

' E 



A' 

'-E 



(1.35) 



During further analysis, the choice between the boundary conditions in the 
forms lll.3fl|l and H1.34II will be a matter of convenience. 



1.3.1 Numerical results 

The solution of Eq. ifOll can be found numerically. To this end, we solve the system 
of two nonhnear equations for the functions Re Z{E) and Im Z{E), using the modified 
Newton method with normalization. 

The solution depends on the bilayer's parameters: the thicknesses of the layers, 
characteristics of materials constituting the bilayer, and the quality of the SN inter- 
face. This dependence enters Eqs. Ill.34|l . H1.35|l via tn and rg. For numerical calcula- 
tions, we assume the characteristics of the bilayer to be the same as in the experiment 
by Kasumov et al. The superconductive layer is made of tantalum, = 5 nm, 
and the normal layer is made of gold, (In — 100 nm. Approximate experimental 
values of the conductivities are jTH] as — 0.01 ^fl^^ cm^^ and crjv = 1 /ii^^^ cm^^. In 
order to calculate the Fermi characteristics of tantalum and gold, we use the values of 
the Fermi energy Ep{Ta.) = 11 eV, Ep{Au) — 5.5 eV, and the free-electrons model. ^ 

^Certainly, the free-electron model does not describe the details of the electronic structure of 
tantalum. However, this simplified model is used in the very derivation of the Usadel equation. 
Also, we do not expect drastic dependence of our results on the Fermi characteristics. To check this, 
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Figure 1.2: Dependence of the order parameter in the S layer A, of the supercon- 
ducting electrons' density in the S layer ns at zero temperature, and of the bilayer's 
critical temperature Tc on the interface resistance per channel ps- All the quantities 
are normalized by the corresponding BCS values. The discrepancy between the curves 
impHes a violation of the BCS relations between A, ns, and Tc. The choice of the 
bilayer's parameters correspond to the experiment so that tsAbcs = 0.016 ps, 
tnAbcs = 0.23 Pb- 
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Figure 1.3: Zoomed part of Fig. 11.21 In the shown range of relatively small resistance 
Pb, the BCS relations between A, ns, and Tc are severely violated. The upper and 
lower graphs differ only in the scaling of the ordinate axis (normal and logarithmic, 
respectively). 
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Figure 1.4: Temperature dependence of A and ns at pB = 90 and 110. The tem- 
perature is normalized by the critical value Tc, which depends on ps; A and ns are 
normalized by their zero-temperature values. The choice of the bilayer's parameters 
implies the relations ts^bcs — 0-016 ps, tnAbcs — 0.23 ps- For comparison, the 
same dependence is also plotted for the BCS case. 
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As a result, we obtain the relations tsAbcs — 0.016pB, tnAbcs — 0.23 pB (so 
that tn/ts ~ 14) , after which the solution of Eq. Ill.34|l depends only on the interface 
resistance pB- Having found the function Z{E) [which is equivalent to finding 0s{E)], 
we start from the case of zero temperature, T = 0, and study the dependence of the 
order parameter A and of the superconducting electrons' density in the S layer ns 
[Eqs. H1.12II . H1.14II ] on pB- The results are plotted in Fig. 11.21 where we also show the 
ps-dependence of the critical temperature Tc, determined by the formula of Ref. 

The suppression of A, ns, and Tc, in comparison to their BCS values in the S 
layer, is a natural consequence of proximity to the normal metal. At the same time, 
there is a possibility of BCS-like behavior, which implies the BCS relations between 
the suppressed quantities and the coincidence of the three curves plotted in Fig. 11.21 
However, the curves split, and the difference between them is largest for relatively 
small values of pB- Figure [T7^ presents the range 80 < < 150 on a larger scale. 

Figure 11.41 shows the temperature dependence of the order parameter A and of 
the superconducting electrons' density in the S layer ng. Although the smaller pB 
the further it is from the BCS limit (corresponding to pB — > oo), we observe that at 
Pb = 90 the curves are closer to the BCS behavior than at pB — 110. An explanation 
of this feature is given in Sec. 11.41 

The DoS in the thin SN bilayer was studied by McMillan Hl^ and Golubov [T3j. 
They demonstrated that (in the case of thin S layer) there is a minigap Eg that is much 
smaller than Abcs (in the case of a bulk superconductor and the perfect interface, 
the minigap was found in Ref. ^H). The gap in the DoS is a property of the bilayer 
as a whole and does not depend on the coordinate, while the energy dependence of the 
DoS is different in the S and N layers. The results of the self-consistent calculation 
of the minigap versus the interface resistance are shown in Fig. 11.51 

The minigap is nonmonotonic, and this fact is already contained in the results of 
McMillan ^iSj, who obtained 

Eg = -^^A at « 1, (1.36) 

tn + ts ts + Tn 

Eg = ^ at » 1. (1.37) 

Tn Ts + Tn 



Formula lll.3(i|l does not determine the dependence of Eg on the interface resistance 
self-consistently; however, if we know that A is suppressed as the interface resistance 
is lowered, then we can conclude that the gap decreases with decreasing pB in the 
corresponding range. At the same time, formula 111.3711 explicitly contains (according 
to the definition of rjv) the inverse ps-dependence. Therefore, Eg reaches a maximum 
at an intermediate pB corresponding to tstnA/{ts + tn) ^ 1, which yields pB ^ 140. 
This estimate agrees with the numerical results (see Fig. ll.5|l . 



we have reproduced all the calculations with slightly different (within 10-15%) values of the Fermi 
energies, and found that the changes in the results amount to, roughly speaking, a rescaling of the 
interface resistance pb within 10%. 
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Figure 1.5: Minigap Eg in the single-particle density of states versus pB- The minigap 
is normalized by the BCS gap value Abcs- Eg is a nonmonotonic function oi ps, 
reaching its maximum aX ps — 160. The inset shows Eg{pB) on a wider (logarithmic) 
scale over pB- The choice of the bilayer's parameters impHes the relations ts^bcs = 
0.016 pi3, tatAbcs = 0.23 

1.4 Anderson limit 

In the Hmit of relatively low interface resistance (following McMillan ^2], we call it 
the Anderson Hmit), the theory describing the bilayer can be developed analytically. 
The condition defining this limit is A, tjv A <C 1. 

Previously, Tc was calculated in this Hmit ^SElEl, and the relation 111.3611 
between the minigap and the order parameter in the S layer A was found, as well 
as the expression for the DoS that has the standard BCS form with the gap Eg 
was obtained [T^. At the same time, the relation between Eg (or A) and the order 
parameter of the isolated S layer /S.bcs was not found. Below we find the minigap 
and other superconductive characteristics of the bilayer as functions of the interface 
resistance. 

First of aH, we need to determine 9{E) [or Z{E)\ solving Eq. H1.30|l [or Eq. H1.34|i ] 
over the entire range of energies E. 

In the region E > A, the solution of Eq. H 1.3411 can be written as Z = 1 + SZ, 
with \SZ\ <^ 1. Keeping terms up to the second order in SZ, we obtain 



5Z^^ 



A{1 -iTNE) 



(1.38) 




This result is valid for arbitrary values of ps. 
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At E < A, the same calculation as for the minigap leads to the result 

iE 

sin 6>g = sin gjv = , (1.39) 

- El 

with the minigap Eg given by Eq. H1.36|l . [To avoid confusion, we note that under 
the less strict limitations for rgA, t^A. used in Eq. Ill.3fi|l . the BCS-Hke result 111.3911 
is valid only up to energies of the order of Eg .] 

Now Im [sin 515] is readily calculated, and in the case of zero temperature, A can 
be found from the self-consistency equation Ijl.l2|l . then the relation ljl.36|l yields 
a formula for Eg. In the limit of the perfect interface (the Cooper limit), which is 
determined by the condition TsTjyuj d / [ts + tat ) <c 1 , we reproduce the classical result 
of Cooper and de Gennes (HE]: 



Eg{pB ^ 0) = 2cji3cxp ( - — ) , (1.40) 



1 

(A) 

with the effective pairing constant 



(A>-^^A. (1.41) 

If the interface is imperfect, tstjmujd/{ts + Tjy) >• 1, then the condition of the An- 
derson limit nevertheless determines a wide range of the interface resistances, where 
we obtain 



Eg 

'^BCS 



^ ^ A -\tn Its 



We emphasize that the Anderson limit does not reduce to the Cooper limit with 
small corrections. On the contrary, due to the relation A <ti ujd, the Cooper limit's 
condition is not satisfied over the most part of the Anderson limit's validity range; 
therefore, the minigap Eg and the quantities calculated below differ drastically from 
the Cooper limit expressions. 

Now we proceed to calculate the density of the superconducting electrons in the 
S layer ng. The main contribution to the integral Ijl.l4ll that determines ng, comes 
from the vicinity of Eg . Inserting the infinitesimal imaginary part of the energy into 
Eq. HI. 3911 . we obtain 

Im [sin^ Os] = ^Eg 5{E - Eg), (1.43) 

which immediately yields the density of the superconducting electrons ng at zero 
temperature: 



ns _ Eg 



Abcs' 



(1.44) 



where Eg is given by Eq. H1.42|l . 

The critical temperature Tc of the bilayer in the Anderson limit [at (rg -|- 
tn)/tstn ^ Tc] was found by Khusainov ^J. Employing our result (11.4211 for 
the minigap, we obtain the following relation: 
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Now we can discuss the general structure of the theory describing the bilayer in 
the Anderson Umit. In the hmit pB ^ 0, our results for the pairing angle 9 (which 
is constant over the entire bilayer, 9 = 9s = 9n) yield expressions which can be 
obtained from the BCS ones II1.16|I - (I1.18|I if we substitute the BCS order parameter 
^BCS by the bilayer's minigap Eg. At pB > 0, corrections to this simple result 
are small while the Anderson limit's conditions are satisfied. Therefore, we obtain a 
BCS-type theory with Eg substituting Abcs in all formulas. 

The results of this section immediately explain the numerical results in the limit of 
relatively small pB, shown in Fig. llJ^I As we have found, the Anderson limit implies 
the following relations between the quantities under discussion: 



Eg = 1.76 T„ 
mas 



(1.46) 
(1.47) 



which substitute Eqs. Ill.22|l and (jl.21|l . For the Ta/Au bilayer to which the numer- 
ical results refer, the Anderson limit is valid at pB < 80 (we see that the values of 
Pb can be large although they are relatively small). Therefore, approaching pB = 80, 
the curves ns/rig'^'^ and Tc/T^'~'^ tend to coincide, while A/ Abcs exceeds them by 
the large factor (r^ + tn)/ts ~ 15. 

The temperature dependence of A and ns, shown in Fig. 11.41 is quite different at 
Pb — 90 and pB ~ 110; at pB = 90, the curves are much closer to the BCS behavior. 
This is also explained by approaching the Anderson limit, where the curves coincide 
with the BCS ones. 

According to McMillan's results QJ], the DoS in the S and N layers coincide in 
the Anderson Hmit (this also follows from Eq. Ill.39|l l. having a BCS-Hke square-root 
singularity at E = Eg. 

Over the whole range of the Anderson limit (including the region tsT]\jlu]j/ (ts + 
Tjv) ~ 1, where the crossover to the Cooper limit takes place), we obtain the result 



E, 



^BCS 



Abcs L , f tstnujd\' 
2lod V 



\Ts +Tn J 



tn/ts 



(1.48) 



which reproduces Eq. Ill.42|l at tstmijJd / {ts + tn) ^ 1 and the Cooper - de Gennes 
result H1.4n|l &t TsTMijJo / [ts+tm) <C 1. In the crossover region T5r7vti'i3/(T5 -I- rjv) ~ 1 
this formula should be considered as an interpolation, because the Debye energy lod 
only determines the order of magnitude for the cutoff energy in the BCS theory. 

A similar calculation for other superconductive characteristics shows that the re- 
lations lll.44|l - lll.47ll between them and the minigap lll.48|l are valid over the whole 
Anderson limit. 



1.5 Parallel critical field 



We proceed to calculate the critical magnetic field He directed along the plane of the 
bilayer. As it was mentioned in Sec. 11.2.21 in the presence of an external magnetic 
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field, the superconducting phase gradient in the Usadel equations H1.7|l - lll.8|l must be 
substituted by its gauge invariant form, which can be expressed via the supercurrent 
velocity v. The spatial distribution of v in the bilayer can be found as follows. 

Let us direct the y axis along the magnetic field H. The supercurrents j — —env 
are directed along the bilayer and perpendicularly to H, i.e., j = {0,0, j{x)) and 
V = (0, 0, v{x)). Near He, the magnetic field inside the bilayer is uniform, so the vector 
potential can be chosen as A = (0, 0, —xH). The supercurrent velocity distribution is 
determined by the equation V x v = eH/m. Another essential point is the continuity 
of V at the SN interface, which follows from the continuity of the superconducting 
phase ip [see the boundary condition The result is 

v{x) = Wo X, (1-49) 

m 

where wo is the supercurrent velocity at the interface, which must be determined from 
the condition that the total charge transfer across the bilayer's cross-section is zero: 

ds 

j{x)dx = 0, (1.50) 



leading to 

2m ) nsds + nj^dN 

The density of the superconducting electrons is constant in each layer (ns and um)- 
Near He, the superconducting correlations are small, 16*1 ^ 1, and the Usadel 
equation (ll.7|l for the paring angle 0(x,E) can be linearized: 



2 dx^ 

Ds 



{iE-2m^DNV^)0N = 0, (1.52) 



-^^ + {tE-2m^Dsv^)9s + \A\ = 0. (1.53) 

At the same time, the second Usadel equation H1.8|l is trivial: its l.h.s. is proportional 
to 

V (sin^ v) = sin 29Vev + sin^ 6 Vv, (1.54) 

where both terms vanish due to the fact that V0 is directed along the x axis whereas 
V is parallel to the z axis. 

The pairing angle 6 is almost spatially constant in each layer; this allows us to 
average each of Eqs. Ijl.52ll - ljl.53|l over the thickness of the corresponding layer, 
obtaining 



d0N 


2dN 


dx 


x=0 Dn 


de^ 


_ 2ds_ 


dx 


x=0 Ds 



where 



iEN-iE)9N, (1.55) 
[{iE-Es)9s + \A\], (1.56) 

En = 2m^DN{^Hx))^, Es = 2m^Ds (v^x)) ^ (1.57) 
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are i/-dependent energies. Using Eqs. Ill.49|l . H1.51|l . we express them via and the 
densities of the superconducting electrons: 



^ 6 



{nsds + UNdN) 



(1.58) 



and En is obtained by the interchange of all the S and N indices. 

Substituting H1.55ll - H1.56|l into the boundary conditions H1.25|l (which should be 
Hnearized), we find 

Zsj^ 

TsEs + tnEm + tstnEsEn - tstnE'^ - iE [ts + tn + tstn {Es + En)] ' 

Os = {1 + tnEn 'tTNE)eN. (1.59) 



The order parameter A cancels out from the self-consistency equation Ill.l2|l . How- 
ever, the resulting equation alone does not suffice for determining Hc{T) because it 
contains Eg and En, which are functions oi un/us- Therefore, to obtain a closed 
system, we must consider the self-consistency equation together with the equation 
determining the ratio UN/ns; the latter equation is obtained from Eq. II1.14|I . The 
resulting system of two nonlinear equations for the quantities He and un /ns is 

riN _ (^N /o°° dE tanh(^) Img^ 

ns ~ crs dE tanh{^)lni0y ^' ^ 

with 9n and 9s given by Eqs. Ijl.59ll . The first equation of the system, Eq. Ijl.(i()|l . 
can be written via the digamma functions, thus taking exactly the same form as Eq. 
Ijl.85ll below (which determines the perpendicular upper critical field) if we denote 
£s^Es + l/rs, £n = En + I/tn- 

In the limit pB oo, Eqs. H1.6flll - H1.61|l lead to the BCS result. In this case, the 
layers uncouple, the density of the superconducting electrons in the N layer vanishes, 
nN /ns 0, and Eq. Ill.(i()|l finally yields 

which determines the parallel critical field Hf'~^^{T) of a thin superconducting film. 



The system of equations lll.6fl|l - Hl.fil|l can be solved numerically at arbitrary val- 
ues of the temperature T and the interface resistance ps; the results for He are 
presented in Figs. 11.61 11.71 

A remarkable feature of the function Hc{pb) at zero temperature (Fig. I1.6|l is 
the steep behavior of He at pb — 120-123. This feature is due to rearrangement of 
the supercurrents inside the bilayer, which occurs in the following way. The super- 
current velocity changes across the thickness of the bilayer according to the simple 
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Figure 1.6: Parallel critical field He, normalized by the BCS value, versus pB at zero 
temperature. The upper and lower graphs differ only in the scaling of the ordinate 
axis (normal and logarithmic, respectively). The nature of the steep behavior of He at 
Pb = 120-123, which is best seen from the lower graph, is explained in the text. The 
inset shows Hc{pb) on a wider scale over pB- The choice of the bilayer's parameters 
implies the relations tsAbcs — 0.016 ps, t^Abcs — 0.23 pB- 
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Figure 1.7: Temperature dependence of the parallel critical field He at pB — 121.2. 
The experimental value of i?c(0), analyzed with the use of the results shown in 
Fig. 11.61 suggests that this value of ps corresponds to the experiment by Kasumov et 
al. j2l The critical field is normaHzed by its zero-temperature value, and the temper- 
ature is normalized by the corresponding Tc. The choice of the bilayer's parameters 
implies the relations ts^bcs = 0.016 ps, tisi/S.bcs — 0.23 ps- For comparison, the 
same dependence is plotted for the BCS case. 
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Figure 1.8: Position of the stationary point xq of the supercurrent distribution versus 
Pb at zero temperature. The coordinate xq is normalized by the S layer thickness ds- 
The fast shift in xq from the center of the S layer at large pb to (nearly) the center 
of the N layer at small pb corresponds to the steep drop in He, shown in Fig. 11.61 
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linear law lll.49|l . This supercurrent distribution may be characterized by the position 
of the stationary point a;o, where the supercurrent velocity is zero: v{xo) — 0, hence 
Xq = mvo/eH. At large values of the interface resistance pB, the density of the super- 
conducting electrons in the N layer is very small, un/us <C 1, and the supercurrents 
circulate only in the S part of the system; this case corresponds to 

xo = y. (1.63) 

Then, while decreasing pB, a shift in xq occurs. Now the supercurrents in the S 
layer are not compensated (in the sense of the charge transfer); therefore, they must 
be compensated by the supercurrents in the N layer, which are enhanced due to 
significant increase in un- This situation corresponds to the beginning of the drop in 
He- The ratio of the superconducting electrons' densities grows rapidly, approaching 
the Anderson limit value Hn/us — (Jn/c^s (see Sec. 11.5.11 belowl: simultaneously, xq 
tends to 

-0 = (1.64) 

and the steep drop in He finishes. For the bilayer to which the numerical results 
refer, ds ^ d^ and as ^ ctn ^ hence Eq. Ill.(i4|l yields xq ~ —d^/'^- 

This scenario is illustrated by Fig. 11.81 which has been obtained numerically. 

The analytical solution of Eqs. Ijl.6()ll - ljl.(il|l at zero temperature in the Anderson 
limit is presented below. 



1.5.1 He at zero temperature in the Anderson limit 

In the zero-temperature Anderson Hmit (defined by the conditions tsEs, t^En <C 1), 
the ratio of the superconducting electrons' densities (jl.(il|l becomes independent of 
the magnetic field, un/hs = (Jn/o's, and the self-consistency equation Ijl.fiOII yields 



TsEs + TnEn _ Abcs 

TS + TN 2 




tjv/ts 

(1.65) 



which determines He (at tstnujd/{ts + tn) ~ 1 this formula should be considered 
as an interpolation — see the discussion in the end of Sec. I1.4|l . This result can be 
compared to the BCS case, which corresponds to the Hmit pB oo. In this case, the 
density of the superconducting electrons in the N layer vanishes, nj^ jns — > 0, and 
the self-consistency equation yields 

Er' - (1.66) 

where Eg'-^^ is given by Eq. Ill.58|l with riM = 0. Finally, we reproduce the result of 
Maki [B] : 

tBCS 



TT iscs ds ' Y ^^^^ ' ( ) 
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where $o = Tr/e is the flux quantum, and S,bcs is the correlation length in the dirty 
limit. 

Remarking that the r.h.s. of Eq. Ill.65|l is identical to Eg/ 2 with the minigap 
Eg given by Eq. (I1.48|l . we see that equation lll.65|l . determining the parallel critical 
field of the bilayer in the Anderson limit, is obtained from the BCS equation lll.66|l if 
we substitute the order parameter Abcs by the minigap Eg (in accordance with the 
results of Sec. ll.4|l and the i/-dependent energy Eg'~'^ by the corresponding averaged 
quantity {tsEs + tnEn)/{ts + tn). 

The explicit result for the parallel critical field of the bilayer, obtained from Eq. 
Ijl.65ll . can be cast into a BCS-like form: 

H, = ^. ,1.08) 

The bilayer's correlation length ^ is the characteristic space scale on which the order 
parameter (or the pairing angle 9, or the Green function) varies in the absence of the 
magnetic field. In the Anderson limit (under discussion), the explicit formula for ^ 
is a natural generalization of the BCS expression [see Eq. Ijl.67|l ] which implies that 
Ds must be substituted by the averaged diffusion constant {D) and Abcs must be 
substituted (in accordance with the results of Sec. 11.411 by the bilayer's characteristic 
energy scale, the minigap Eg [Eq. (|1.48|l ]: 



\ Eg Ts + Tn 



The effective thickness of the bilayer in Eq. (Il.fi8|l is 



\2 



{asds + ONdN) {<ysd% + UNd]^) + iasONdsdN {ds + dw) 
^— ^ . (1.70) 

In the case of equal conductivities, as — <Jn, the effective thickness is simply the 
geometrical one: dos = ds + d^- This case corresponds to a uniform density of 
the superconducting electrons, ns = nj^f, which implies a continuous distribution 
of the supercurrents, centered at the middle of the bilayer [this can also be seen 
from Eq. Ill.fi4|l which yields xq = {ds — dN)/2 in the case as = CAf]. However, 
in a more subtle situation when the conductivities are different, the density of the 
supercurrent experiences a jump at the SN interface; this nontrivial supercurrent 
distribution results in the nonequivalence of dcff to the geometrical thickness of the 
bilayer. 



1.6 Perpendicular upper critical field 

Now we turn to calculating the upper critical field Hc2 perpendicular to the plane of 
the bilayer. 

As in the case of the parallel critical field, we start with discussing the supercur- 
rent distribution, which is now a function of the sample boundaries in the yz plane. 
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perpendicular to the magnetic field H (the magnetic field is directed along the x 
axis) . The infinite bilayer under consideration can be thought of as a disk of a large 
radius; let us assume y = 0, z = at the axis of the disk. Then the supercurrent 
distribution is axially symmetric, and, with the gauge chosen as A = [Hr] /2, the 
superconducting phase must be constant, ip = 0, which yields a simple result for the 
supercurrent velocity: v = eA/m. 

Near Hc2, the superconducting correlations are small, \6\ <^ 1, and the Usadel 
equations can be linearized: 

+ 2eAfe + iE9 + A = 0, (1.71) 
Ave = 0. (1.72) 

The second of these equations is trivially satisfied because 0{r) is axially symmetric. 

Thus, the Usadel equations reduce to the single Eq. Ijl.71ll for the pairing angle 
6{r, E). Introducing the cylindrical coordinates r ^ (x, p) and denoting P = — iVp + 
2eA(p), we rewrite this equation as 

DmB^O^ ^^p20^+^i^0^ = O, (1.73) 



2 2 



Ds d^Os Ds 
2 dx^ 2 



-p'^Os+iEOs + A^O. (1.74) 



We cannot solve these equations straightforwardly because near the upper critical 
field, the order parameter A(p) is a nontrivial unknown function of the in-plane 
coordinate p (while the x-dependence is absent due to the small thickness of the 
bilayer). In this situation, we employ the following approach. 

Averaging each of Eqs. Ijl .7811 - 11 1.7411 over the thickness of the corresponding layer, 
we obtain 





2dN 


dx 




dO^ 


_ 2ds 


dx 


x=0 



P^9n -iEON] , (1.75) 



^P^Os + iEOs + Ay (1.76) 

The averaged pairing angles entering the r.h.s. of Eqs. Ill.75|l - (ll.76ll are 

1 

9Nip,E)^— dx9N{x,p,E), (1.77) 

esip,E)^^ [ ' dx9six,p,E). (1.78) 
d-s Jq 

Substituting Eqs. Ijl .7511 - 111.7611 into the boundary conditions lll.25|l (which should be 
linearized), we obtain a system of two differential equations for the function 6{p, E): 

TN (^P^On - iE9n] = Ts (~^P^9s + iE9s + a\^9s- 9^. (1.79) 
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From the vicinity of the superconductive transition it follows that the pairing 
angle 9 depends on the order parameter A linearly: 



Os{p,E) 



A(p) ^ 

aN (E) ' 
A(p) 

as (E) ' 



where the functions aN{E) and as{E) are spatially independent. Then Eqs. 
can be rewritten as 



^P^A(P) 



Ds^2 



P^A(p) 



iE 



iE+ — 

TN 

1 

TS 



aN [E) 
as (E) 
as{E) 
aN (E) 



A(p), 



as{E) 



A(p). 



(1.80) 
(1.81) 

(1.82) 
(1.83) 



We see that the order parameter must be an eigenfunction of the differential 
operator P^. Moreover, in order to obtain the largest value of Hc2, we should choose 
the eigenfunction corresponding to the lowest eigenvalue (in complete analogy with 
Refs. ^1E|)- The solution of the emerging eigenvalue problem is readily found 
thanks to its formal equivalence to the problem of determining the Landau levels of 
a two-dimensional particle with the "mass" 1/D and the charge — 2e in the uniform 
magnetic field H directed along the third dimension. The lowest Landau level is 
DeH; the function as{E) is straightforwardly determined, 



as{E) ^ DseH - iE 



tn{Dn eH -iE) 
TS [1 + TN {Dn eH - iE)] ' 



(1.84) 



and we substitute 9s{p,E) into the self-consistency equation Ill.l2|l . The order pa- 
rameter A(p) cancels out, and the resulting equation, which determines Hc2{T), can 
be cast into the form 



In- 



T, 



BCS 



T 

£s—£n 



where 



y^{£s-£N)''+i/TSTN 

£s—£n 

V(£s-£jv)H4/rsr„ 




i£s - £n)" + 



1 

4irT 



N 



{£s - £. 



N) 



£s = Ds eH,2 + — 



1 

TS ' 



£. 



N 



Dn eHc2 + — 

TN 



(1.85) 



(1.86) 



are i?-dependent energies. The logarithmic term in the r.h.s. of Eq. Ijl.85ll takes 
account of the fmiteness of the Debye energy wd; it becomes important only in the 
limit of the perfect interface (the Cooper limit), i.e., when tst^wd / {ts + tn) ^ 1. At 
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Figure 1.9: Perpendicular upper critical field Hc2, normalized by its BCS value, versus 
Pb sX zero temperature. The upper and lower graphs differ only in the scaling of the 
ordinate axis (normal and logarithmic, respectively). The inset shows Hc2{pb) on 
a wider scale over ps- The choice of the bilayer's parameters implies the relations 
TsAscs = 0.016 Ps, tatAbcs 0.23 ps. 



tstn^d/{ts +tm) ^ 1 the logarithmic term should be considered as an interpolation 
— see the discussion in the end of Sec. 11.41 

In the limit pB ^ co, Eq. (jl.85|l reproduces the classical result of Maki and 
de Gennes ^HI for the BCS case (see also the book |20jl: 



T " V2 27rr J " V2, 

which is valid for bulk superconductors and superconductive layers of arbitrary thick- 
ness (when the magnetic field is directed perpendicularly to them). 

Equation (jl.85ll can be solved numerically at arbitrary values of the temperature 
T and the interface resistance ps; the results for Hc2 are presented in Figs. ll.9lll.10l 

The analytical solution of Eq. (I1.85|l at zero temperature in the Anderson limit is 
presented below. 



42 



Chapter 1 



1.0 - 












0.8 - 


■ \\ 






\\ 






■ 






\\ 




0.6 - 


■ \\ 






■ \\ 
X\ 




0.4 - 




x\ 




Pg = 80; BCS 


\\ 




■ N> 


0.2 - 


Pb = 110 


■ ■ X 




■ ■ . X 




Pb = 125 


0.0 - 







0.0 0.2 0.4 0.6 0.8 1.0 



Figure 1.10: Temperature dependence of the perpendicular upper critical field Hc2 at 
Pb = 80, 110, 125, and in the BCS case. In each case, the critical field is normalized by 
its zero-temperature value, and the temperature is normalized by the corresponding 
Tc. According to the results of Sec. ll.4L the curves in the BCS and Anderson (ps = 80) 
limits coincide. At intermediate values of pB, the curves can lie both above (ps = 110) 
and below {pb — 125) the BCS curve. The choice of the bilayer's parameters implies 
the relations tsAbcs = 0.016 ps, tnAbcs = 0.23 ps. 
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1.6.1 Hc2 at zero temperature in the Anderson limit 

In the zero-temperature Anderson limit (defined by the conditions DseHc2 ^ l/^'s, 
DNeHc2 < I/tat), Eq. ltn?K|l yields 



Hc2 (ts + Tn) Ds 



2lOd V \TS +TN J 



(1- 



(at tstnujd/{ts + Tn) ~ 1 this formula should be considered as an interpolation — 
see the discussion in the end of Sec. I1.4II . Here the zero-temperature BCS value of 
the upper critical field, as follows from Eq. H1.87|l . is 



„BCS _ ^BCS _ '^0 an\ 
^c2 - r,^n - 77-72 ■ (1-^9) 



It is instructive to rewrite the perpendicular upper critical field of the bilayer (jl.88ll 
in the standard BCS-like form 

H.2 = (1.90) 

where ^ is the bilayer's correlation length given by Eq. Ijl.69ll [the physical interpre- 
tation of this result for ^ precedes Eq. H1.69II ]. 

In the Cooper limit {i.e., at tstni^d / {ts + tn) ^ 1) Eqs. H1.88II . (I1.9fl|l reproduce 
the result of Refs. [HIISI- 



1.7 SNS, NSN, SNINS, NSISN, and superlattices 

Our results for A, ns, Tc, Eg, and Hc2 {i.e., all the results except He) can be directly 
appHed to more complicated structures such as SNS and NSN trilayers, SNINS and 
NSISN systems, and SN superlattices. 

Let us consider, for example, a symmetric SNS trilayer consisting of two identical 
S layers of thickness ds separated by a N layer of thickness 2dN- The SN interfaces 
can have arbitrary (but equal) resistances. As before, the x axis is perpendicular to 
the plane of the structure. This trilayer can be imagined as composed of two identical 
bilayers perfectly joined together along the N sides. Indeed, the pairing angle 9 has 
zero x-derivative on the outer surfaces of the bilayers, thus producing the correct 
(symmetric in the x-direction) solution for 9 in the resulting trilayer. Consequently, 
the symmetric SNS trilayer has exactly the same physical properties [A, ns, Tc, Eg, 
Hc2] as the SN bilayer considered in the present chapter. The only point where the 
above reasoning fails is the calculation of the parallel critical field He- In this case, 
the combination of the supercurrent distributions in the two bilayers does not yield 
the correct distribution in the resulting SNS trilayer, which implies that the Usadel 
equations for the two systems are different. 

Evidently, the above reasoning, based on the formal equivalence of the outer- 
surface boundary condition for the bilayer to the symmetry-caused condition in the 
middle of the SNS trilayer, also holds for symmetric NSN trilayers (N layers of thick- 
ness dN, S layer of thickness 2ds, identical SN interfaces) and SN superlattices (N 
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layers of thickness 2dN, S layers of thickness 2^5, identical SN interfaces). Moreover, 
the same applies to systems composed of two bilayers in nonideal contact with each 
other: SNINS and NSISN (where I stands for an arbitrary potential barrier), because 
the presence of a potential barrier does not violate the applicability of the symmetry 
argument. Thus, all the results obtained for the bilayer (except He) are also valid for 
these structures. 



1.8 Discussion 

Now we turn to a possible experimental application of our results. Our results provide 
a method for determining ps, a very important parameter of the bilayer which is not 
directly measurable. By analyzing the experimental ^ values Tc — 0.4 K and 
He = 0.1 T, we get pb ~ HI and ps ~ 121, respectively. Within the experimental 
accuracy of the bilayer's parameters, the two estimates for pB should be considered 
close. Interestingly, the value pB ~ 121 extracted from the measured value of He 
corresponds to the extremely narrow region of the steep drop in Hc{pb) (see Fig. ll.6|l . 

An essential property of the bilayer used throughout the chapter is its small thick- 
ness. Now we shall argue that the bilayer employed in the experiment by Kasumov 
et al. [2 (and to which our numerical results refer) can be considered thin. The Us- 
adel equations lll.23|l - ljl.24ll imply that the characteristic spatial scale of the bilayer's 
properties variation is \/Dn\s/Eo for the N and S layers, respectively. Here Eo is 
the characteristic energy scale for the self-consistency equation. As the interface re- 
sistance decreases from infinity to zero, first increases from Abcs to (r^^ +'''N^)y 
and then — to ujd (the crossovers between these regimes occur at such pb that the 
corresponding expressions are of the same order) . The above estimates show that the 
experiment corresponds to the case Eq ~ Abcs- Therefore, the thicknesses of the 
layers (c?Ar = 100 nm and ds = 5 nm) must be small compared to D^^s / ^bcs, 
which equals 194 nm and 16 nm for the N and S layer, respectively. We thus conclude 
that the condition of thin bilayer is approximately satisfied. 

Finally, we wish to remark on a pecuHarity of real systems which can be relevant 
when one compares our findings with an experiment. The point is that during the fab- 
rication of a bilayer, the interface between S and N materials cannot be made ideally 
uniform. In other words, the local interface resistance possesses spatial fluctuations. 
At the same time, as we have shown, the bilayer's properties are highly sensitive to 
the interface quality, which could lead to complicated behavior not reducing to the 
simple averaging of the interface resistance embodied in p^- One possibility could be 
a percolation-Hke proximity effect. We leave the study of inhomogeneity effects for 
further investigation. 

1.9 Conclusions 

In this chapter, we have studied, both analytically and numerically, the proximity 
effect in a thin SN bilayer in the dirty limit. The strength of the proximity effect is 
governed by ps, the resistance of the SN interface per channel. 
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The quantities calculated were A, the order parameter; ns, the density of the 
superconducting electrons in the S layer; Tc, the critical temperature; Eg, the minigap 
in the density of states; He and Hc2, the critical magnetic field parallel to the bilayer 
and the upper critical field perpendicular to the bilayer. 

These quantities were calculated numerically over the entire range oi ps- For 
this purpose, the characteristics of the bilayer were assumed to be the same as in 
the experiment by Kasumov et al. (Ta/Au bilayer, ds/dN = 1/20). In the Hmit 
of an opaque interface, A, ns, Tc, He, and Hc2 approach their BCS values. At the 
same time, Eg does not coincide with the order parameter A, and Eg ^ when 
Pb ^ oo, although in general, the energy dependence of the DoS in the S and N 
layers approaches the BCS and normal-metal results, respectively. 

The minigap Eg demonstrates nonmonotonic behavior as a function oi pB- Analyt- 
ical results for the two Hmiting cases of small and large pB show that in the Anderson 
limit. Eg increases with increasing pB, whereas in the limit of an opaque interface. 
Eg tends to zero. Thus, Eg reaches its maximum in the region of intermediate pB- 

Also in the region of moderate resistances, a jump of the parallel critical field (due 
to a redistribution of supercurrents in the bilayer) is discovered. 

The most interesting case of relatively low interface resistance (the Anderson 
limit) has been considered analytically. The simple BCS relations between A, ns, 
Tc, He, Hc2 are substituted by similar ones with Eg standing instead of A. The 
relation between the minigap Eg and the order parameter A in this limit is expressed 
by Eq. implying that in the case where rg < tn, the BCS relations are 

strongly violated (by more than the order of magnitude for the above-mentioned 
Ta/Au bilayer). The DoS in the S and N layers coincide, showing BCS-like behavior 
with the standard peculiarity at E = Eg. It should be emphasized that absolute 
values oi pb corresponding to the Anderson limit can be large; for the experiment [2] 
this Hmit is already valid at pB < 80. 

All the results (except He) obtained for the bilayer also apply to more compli- 
cated structures such as SNS and NSN trilayers, SNINS and NSISN systems, and SN 
superlattices. 
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Proximity effect in SF systems 



2.1 Nonmonotonic critical temperature in SF bilay- 
ers 

2.1.1 Introduction 

Superconductivity and ferromagnetism are two competing orders: while the former 
"prefers" an antiparallel spin orientation of electrons in Cooper pairs, the latter forces 
the spins to align in parallel. Therefore, their coexistence in one and the same material 
is possible only in a narrow interval of parameters; hence the interplay between 
superconductivity and ferromagnetism is most conveniently studied when the two 
interactions are spatially separated. In this case the coexistence of the two orders is 
due to the proximity effect. Recently, much attention has been paid to properties of 
hybrid proximity systems containing superconductors (S) and ferromagnets (F); new 
physical phenomena were observed and predicted in these systems PQ Ql IS ^ ID lEl ■ 
One of the most striking effects in SF layered structures is highly nonmonotonic 
dependence of their critical temperature Tc on the thickness dp oi the ferromagnetic 
layers. Experiments exploring this nonmonotonic behavior were performed previously 
on SF multilayers such as Nb/Gd 0, Nb/Fe jEj, V/V-Fe [HI, and Pb/Fe but 
the results (and, in particular, the comparison between the experiments and theories) 
were not conclusive. 

To perform reliable experimental measurements of Tc{dp), it is essential to have 
dp large compared to the interatomic distance; this situation can be achieved only 
in the Hmit of weak ferromagnets. Active experimental investigations of SF bilayers 
and multilayers based on Cu-Ni dilute ferromagnetic alloys are carried out by several 
groups [111 112) . In SF bilayers, they observed nonmonotonic dependence Tc{dp). 
While the reason for this effect in multilayers can be the O-tt transition P], in a 
bilayer system with a single superconductor this mechanism is irrelevant, and the 
cause of the effect is interference of quasiparticle, specific to SF structures. 

In the present chapter, motivated by the experiments of Refs. ^iJEl we theoreti- 
cally study the critical temperature of SF bilayers. Previous theoretical investigations 
of Tc in SF structures were concentrated on systems with thin or thick layers (com- 
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pared to the corresponding coherence lengths); with SF boundaries having very low or 
very high transparencies; the exchange energy was often assumed to be much larger 
than the critical temperature; in addition, the methods for solving the problem were 
usually approximate OElEIinS, 14, 13^10]. The parameters of the experiments 
of Refs. |ll[ll2j do not correspond to any of the above limiting cases. In the present 
chapter we develop two approaches giving the opportunity to investigate not only 
the Hmiting cases of parameters but also the intermediate region. Using our methods, 
we also confirm different types of nonmonotonic Tcidp) behavior, found previously 
OEIEI: the minimum of Tc and the reentrant superconductivity. Comparison of 
our theoretical predictions with the experimental data shows good agreement. 

A number of methods can be used for calculating Tc- When the critical temper- 
ature of the structure is close to the critical temperature Tcs of the superconductor 
without the ferromagnetic layer, the Ginzburg-Landau (GL) theory applies. How- 
ever, Tc of SF bilayers may significantly deviate from Tcs, therefore we choose a 
more general theory valid at arbitrary temperature — the quasiclassical approach 
[TtI [Tsl [Tfl] . Near Tc the quasiclassical equations become Hnear. In the literature 
the emerging problem is often treated with the help of the so-called "single-mode" 
approximation |14l 1151 0] , which is argued to be qualitatively reasonable in a wide 
region of parameters. However, this method is justified only in a specific region of pa- 
rameters which we find below. Moreover, below we show examples when this method 
fails even qualitatively. Thus there is need for an exact solution of the linearized 
quasiclassical equations. The limiting case of perfect boundaries and large exchange 
energy was treated by Radovic et al. [3;- 

Based on the progress achieved for calculation of Tc in SN systems (where N de- 
notes a nonmagnetic normal material) , we develop a generaHzation of the single- 
mode approximation — the multimode method. Although this method seems to be 
exact, it is subtle to justify it rigorously. Therefore we develop yet another approach 
(this time mathematically rigorous), which we call "the method of fundamental so- 
lution". The models considered previously El ^1 ^1 ^1 \M correspond to 
Hmiting cases of our theory. 



2.1.2 Model 

We assume that the dirty-limit conditions are fulfilled, and calculate the critical 
temperature of the bilayer within the framework of the linearized Usadel equations for 
the S and F layers (the domain < x < ds is occupied by the S metal, —dp < x < 
— by the F metal, see Fig. I2.1|l . Near Tc the normal Green function is G = sgncjn, 
and the Usadel equations for the anomalous function F take the form 

fsT^Tcs^^ - \uJn\Fs + A = 0, 0<x<ds, (2.1) 
d^F 

ilnTcs-r^ - i\^n\ + ihsgnuJn)FF = 0, -dp < x < 0, (2.2) 



dx'^ 



Aln^-T5:(^-F,) (2.3) 
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"T" 




Figure 2.1: SF bilayer. The F and S layers occupy the regions 
< X < ds, respectively. 



'dp < X < and 



(the order parameter A is nonzero only in the S part). Here = \/ Ds/'^t^Tcs, 
^p = Dp /2-kTcs are the coherence lengths, while the diffusion constants can be 
expressed via the Fermi velocity and the mean free path: D = uZ/3; w„ = 7rT(2n + 1) 
with n = 0, ±1, ±2, . . . are the Matsubara frequencies; h is the exchange energy; and 
Tcs is the critical temperature of the S material. Fg^^p^ denotes the function F in the 
S(F) region. We use the system of units in which Planck's and Boltzmann's constants 
equal unity, h = = 1- 

Equations H2.1ll - H2.3|l must be supplemented with the boundary conditions at the 
outer surfaces of the bilayer: 

dFsjds) ^ dFpj^dp) ^ 
dx dx ' 

as well as at the SF boundary pT]: 

dFsjO) dFp jO) ps^s 

is—; =7?F — -j , 7 = — ^, (2-5) 

dx dx Pf£,f 

iFlB^^^FsiO)~FpiO), 1B = ^. (2.6) 

dx Pfc,f 

Here ps, pp are the normal-state resistivities of the S and F metals, Rb is the re- 
sistance of the SF boundary, and A is its area. The above boundary conditions 
were derived for SN interfaces (N is a normal metal); their use in the SF case 
is justified by the small parameter h/E-p <^ 1 {Ep is the Fermi energy). Indeed, 
the interface has atomic (by the order of magnitude) thickness. While the exchange 
energy is small compared to the Fermi energy, the characteristic length for magnetic 
properties is much larger than the atomic scale. Therefore the boundary conditions 
are determined by the properties of the interface itself but not by the properties of 
the contacting metals. In the limit of strong ferromagnets, this condition fails and the 
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Kupriyanov-Lukichev boundary conditions lose its validity. However, in the thesis 
we consider only weak ferromagnets. 

The Usadel equation in the F layer is readily solved: 



Fp = C(aj„) cosh (kh[x + dF]], (2.7) 




and the boundary condition at a; = can be written in closed form with respect to 

- , PsiOl (2.8) 



khip t&ii\i{khdF] 



This boundary condition is complex. In order to rewrite it in a real form, we do 
the usual trick and go over to the functions 

= F{un) ± F{-LO„). (2.9) 

According to the Usadel equations Ij2.1ll - lj2.3|l . there is the symmetry F(— u;„) = 
F*{LUn) which impHes that F^ is real while F~ is a purely imaginary function. 

The symmetric properties of F+ and F~ with respect to a;„ are trivial, so we 
shall treat only positive w„. The self-consistency equation is expressed only via the 
symmetric function F^: 

Ai„^.,ri:(H^-f^), (2^,0) 

and the problem of determining Tc can be formulated in a closed form with respect 
to Fg as follows. The Usadel equation for the antisymmetric function Fg does 
not contain A, hence it can be solved analytically. After that we exclude Fg from 
boundary condition H2.8|l and arrive at the effective boundary conditions for Fg-. 



where 



ur( \ A5(7B+Reg^)+7 
= fcs^stanh(fc5d5), ks^-^y ^" 



is V ""^cS 

The self-consistency equation Il2.1fl|l and boundary conditions (I2.11|I - II2.12|I . together 



with the Usadel equation for Fg- 



ilnT,s-^ - iOnF+ + 2A = (2.13) 
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will be used below for finding the critical temperature of the bilayer. 

The problem can be solved analytically only in limiting cases (see Appendix El- 
In the general case, one should use a numerical method, and below we propose two 
methods for solving the problem exactly. 



2.1.3 Multimode method 

Starting point: the single-mode approximation and its applicability 

In the single-mode approximation (SM A) one seeks the solution of the problem Il2.1fl|l - 
H2.13|l in the form 

F+{x,ujn) = f(c^„)cos (^^^) ' (2-14) 

A(x) = S cos . (2.15) 

This anzatz automatically satisfies boundary condition II2.11|I at x — ds. 
The Usadel equation H2.13|l yields 

fK) = Jlr ' (2.16) 

then the self-consistency Eq. (12.1011 takes the form {S and fl do not depend on ujn) 

■»t-(^?lr)-G)- 

where is the digamma function. 

Boundary condition II2.11|I at a; = yields 

ntan(^n^^ ^W{uJn)- (2.18) 

The critical temperature Tc is determined by Eqs. Ij2.17ll . (I2.18|l . 

Although this method is popular, it is often used without pointing out the limits of 
its appHcability. The explicit formulation of the corresponding condition is as follows 
[13] : the single- mode method is correct only if the parameters of the system are such 
that W can be considered w„-independent [because the left-hand side of Eq. (12.1811 
must be cJn-independent]. 

Appendix [0 demonstrates examples of the SMA validity and corresponding ana- 
lytical results. 

In one of experimentally relevant cases, h/nTcs > 1, d_F ~ Cf, the SMA is appli- 
cable if \/h/TTTcs ^ l/7s (see Appendix [0 for details). 



Inclusion of other modes 

The single-mode approximation impHes that one takes the (only) real root SI of Eq. 
Ij2.17|l . An exact (multimode) method for solving problem (I2.1()|I - (I2.13II is obtained 
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Figure 2.2: Plot of the ■0(1/2 + 2:) — V'lV^) function. At any positive y, the equation 
y = ^^{1/2 + z) — ■0(1/2) has one positive and infinite number of negative solutions. 
The positive solution is employed in the single- mode approximation [see Eq. II2.17II ]. 
while taking into account the negative solutions [which implies imaginary fl in Eq. 
Ij2.17ll ] we obtain the exact multimode method. 



if we also take imaginary roots into account — there is infinite number of these |2fH . 
see Fig. lO 

Thus we seek the solution in the form 

f x-ds\ ^ cosh(n„^^ 
Fg{x,ujn) = fo(w„)cos ^lo^^ + > ^ (2.19) 

V J cosh(l7„,|fj 

/ x-d<,\ A COshf^nM^) 

A(x)=<5ocosk!o^^ ^7 (2.20) 

V J cosh(ri„|fj 

(The normalizing denominators in the cosh-terms have been introduced in order to in- 
crease accuracy of numerical calculations.) This anzatz automatically satisfies bound- 
ary condition H2.11|l at a; = ^5. 

Substituting the anzatz [Eqs. Ij2.19ll - I|2.2flll ] into the Usadel equation Ij2.13ll . we 
obtain 



25o 



foK) = , . , oL^ - ' (2.21) 



fm(Wn) = rj. , m=l,2,. 



then the parameters fl are determined by the self-consistency equation Ij2.1()ll {S and 
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fl do not depend on w„): 



From Eqs. II2.22|I and properties of the digamma function ^] it follows that the 
parameters fl belong to the following intervals: 

J-cS J^cS 

where C w 0.577 is Euler's constant. 

Boundary condition II2.11|I at a; = yields the following equation for the ampli- 
tudes 6: 

W{un) COS {^ods/£.s) - ^0 sin {^ods/ts) 

^0 . ■ .02. 



m—1 

The critical temperature Tc is determined by Eqs. H2.22II and the condition that Eq. 
H2.24II has a nontrivial (w„-independent) solution with respect to S. 

Numerically, we take a finite number of modes: m = 0, 1, . . . , M. To take account 
of Wn-independence of the solution, we write down Eq. Ij2.24ll at the Matsubara 
frequencies up to the TVth frequency: n = 0, 1, . . . , TV. Thus we arrive at the matrix 
equation KnmSm = with the following matrix K: 

^ _ Wjujn) cos {ilpds/^s) - ^0 sin {flods/£.s) 

"° LUn/TrTcS + 

^ _ + tanh (Slr,idslis) 

n = 0, 1, . . . ,7V, 771= 1,2,..., M. 

We take M = N, then the condition that Eq. Ij2.24ll has a nontrivial solution takes 
the form 

deti^ = 0. (2.26) 

Thus the critical temperature Tc is determined as the largest solution of Eqs. 
12221, (l2:2fijl . 



2.1.4 Method of fundamental solution 

By definition, the fundamental solution Q{x,y;LJn) (which is also called the Green 
function) of problem Il2.11ll - lj2.i;^ll satisfies the same equations, but with the delta- 



56 



Chapter 2 



functional "source" |23] : 

S.s'^Tcs '^ ^^^2^^ - uJnG{x,y) = ~6{x-y), (2.27) 

is'^-WMQiO^y), '^^-0. (2.28) 

The fundamental solution can be expressed via solutions vi, V2 of Eq. II2.27|I with- 
out the delta-function, satisfying the boundary conditions at x = and x = ds, 
respectively: 

y[x,y,u;n) ,^^^ksds) + {W/ksis) cosh (ksds) \v2{x)vi{y), y^x'^ ' 
where 

vi(x) = cosh(fcsa:) + (W/ks^s) sinh(fcsx), (2.30a) 
V2ix) = cosh {ks[x - ds]) ■ (2.30b) 

Having found Q{x,y;u!n), we can write the solution of Eqs. H2.llll - H2.13ll as 

F+{x;uJn)^2 f g{x,y;uJn)A{y)dy. (2.31) 
Substituting this into the self-consistency equation Ij2.1()ll . we obtain 



A(x)ln^ = 2^r, J2 



Tc 



A(x) f^"^ 

g{x,y;uJn)A{y)dy 



(2.32) 



This equation can be expressed in an operator form: Aln(Tcs/Tc) = LA. Then the 
condition that Eq. Ij2.32ll has a nontrivial solution with respect to A is expressed by 
the equation 

det - i In = 0. (2.33) 

The critical temperature Tc is determined as the largest solution of this equation. 

Numerically, we put problem Ij2.32ll . I|2.33ll on a spatial grid, so that the linear 
operator L becomes a finite matrix. 



2.1.5 Numerical results 

In Sees. 12.1.31 12.1.41 we developed two methods for calculating the critical temper- 
ature of a SF bilayer. Specifying parameters of the bilayer we can find the critical 
temperature numerically. It can be checked that the multimode method and the 
method of fundamental solution yield equivalent results. However, at small temper- 
atures Tc ^ Tcs, the calculation time for the multimode method increases. Indeed, 
the size of the matrix K [Eq. Il2.25|l ] is determined by the number N of the maximum 
Matsubara frequency ujn, which must be much larger than the characteristic energy 
nTcs', hence N ^ Tcs/Tc- Therefore, at low temperatures we use the method of 
fundamental solution. 
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Experiment 
Theory 



3^ • , • 1 ' 1 • h 

5 10 15 20 
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Figure 2.3: Theoretical fit to the experimental data of Ref. pi]. In the experiment, 
Nb was the superconductor (with ds = 11 nm, Tcs = 7K) and Cuo.43Nio.57 was the 
weak ferromagnet. From our fit we estimate h « 130 K and ~ 0.3. 

Comparison with experiment 

Using our methods we fit the experimental data of Ref. ^I]; the result is presented 
in Fig. 12.31 Estimating the parameters ds = 11 nm, Tcs = 7K, ps = 7.5/^f2cm, 
£,s = 8.9 nm, pp = GOpflcm, = 7.6 nm, 7 = 0.15 from the experiment |2|j, and 
fitting only h and 75, we find good agreement between our theoretical predictions 
and the experimental data. 

The fitting procedure was the following: first, we determine h « 130 K from the 
position of the minimum of Tc{dp); second, we find jb ~ 0.3 from fitting the vertical 
position of the curve. 

The deviation of our curve from the experimental points is small; it is most pro- 
nounced in the region of small dp corresponding to the initial decrease of Tc- This is 
not unexpected because, when dp is of the order of a few nanometers, the thickness 
of the F film may vary significantly along the film (which is not taken into account in 
our theory) , and the thinnest films can even be formed by an array of islands rather 
than by continuous material. At the same time, we emphasize that the minimum 
of Tc takes place at dp « 5nm, when with good accuracy the F layer has uniform 
thickness. 

Various types of Tc{dp) behavior 

The experimental results discussed above represent only one possible type of Tc{dp) 
behavior. Now we address the general case; we obtain different kinds ofTc{dp) curves 
depending on parameters of the bilayer. 

To illustrate, in Fig. 12. 41 we plot several curves for various values of 75 [we recall 
that 7b oc Rb, where Rb is the resistance of the SF interface in the normal state — 
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Figure 2.4: Characteristic types of Tc{dp) behavior. The thickness of the F layer is 
measured in units of the wavelength defined in Eq. II2.41|I . The curves correspond 
to different values of 7b- The exchange energy is /i = 150K; the other parameters are 
the same as in Fig. VI.'M One can distinguish three characteristic types of Tc{dp) be- 
havior OElEl: 1) nonmonotonic decay to a finite Tc with a minimum at particular 
dp (76 = 2; 0.5; 0.1; 0.07), 2) reentrant behavior (73 = 0.05; 0.02), 3) monotonic de- 
cay to Tc = at finite dp (73 = 0). The bold points indicate the choice of parameter 
corresponding to Fig. 12.71 

see Eq. H2.6|l ]. The exchange energy is ft. = 150 K; the other parameters are the same 
as in Fig.lOl 

We confirm three characteristic types of Tc{dp) behavior, found previously in 
limiting cases or by approximate methods in Refs. ^,15, 16j 1) at large enough in- 
terface resistance, Tc decays nonmonotonically to a finite value exhibiting a minimum 
at a particular c?f, 2) at moderate interface resistance, Tc demonstrates the reentrant 
behavior: it vanishes in a certain interval of dp^ and is finite otherwise, 3) at low 
enough interface resistance, Tc decays monotonically vanishing at finite dp. A similar 
succession oiTddp) curves as in Fig. l2.4l ca,n be obtained by tuning other parameters, 
e.g., the exchange energy h or the normal resistances of the layers (the parameter 7). 

A common feature seen from Fig. 12.41 is saturation of Tc at large dp > Xh. This 
fact has a simple physical explanation: the suppression of superconductivity by a 
dirty ferromagnet is only due to the effective F layer with thickness on the order of 
Xh, adjacent to the interface (this is the layer explored and "felt" by quasiparticles 
entering from the S side due to the proximity effect). 

It was shown by Radovic et al. that the order of the phase transition may 
change in short-periodic SF superlattices, becoming the first order. We also observe 
this feature in the curves of types 2) and 3) mentioned above. This phenomenon 
manifests itself as discontinuity of Tc{dp): the critical temperature jumps to zero 
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Figure 2.5: Change of the phase transition's order. This phenomenon manifests itself 
as discontinuity of Tc{dp): the critical temperature jumps to zero abruptly without 
taking intermediate values. Formally, becomes a double-valued function, but the 
smaller solution is physically unstable (dotted curve). For illustration we have chosen 
the curve from Fig. corresponding to 75 = 0.05. 

abruptly without taking intermediate values (see Figs. 12. 41 12. 5|l . Formally, Tc becomes 
a double-valued function, but the smaller solution is physically unstable (dotted curve 
in Fig. ESI- 

An interesting problem is determination of the tricritical point where the order 
of the phase transition changes. The corresponding result for homogeneous bulk 
superconductors with internal exchange field was obtained a long time ago in the 
framework of the Ginzburg-Landau theory However, the generalization to the 
case when the GL theory is not valid has not yet been done. We note that the 
equations used in Refs. OEI were applied beyond their applicability range because 
they are GL results vaHd only when Tc is close to Tcs- 

Comparison between single- and multimode methods 

A popular method widely used in the literature for calculating the critical temperature 
of SF bi- and multi-layers is the single-mode approximation. The condition of its 
vaHdity was formulated in Sec. 12.1.31 However, this approximation is often used for 
arbitrary system's parameters. Using the methods developed in Sees. 12.1. 31 12.1.41 
we can check the actual accuracy of the single-mode approximation. The results are 
presented in Fig. 12.61 

We conclude that although at some parameters the results of the single-mode 
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Figure 2.6: Comparison between single- and multimode methods. The parameters 
are the same as in Fig. 12.41 Generally speaking, the results of the single-mode and 
multimode (exact) methods are quantitatively and even qualitatively different: b), c), 
d), and e). However, sometimes the results are close: a) and f). Thus the single- mode 
approximation can be used for quick estimates, but rehable results should be obtained 
by one of the exact (multimode or fundamental-solution) techniques. 
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and multimode (exact) methods are close (Figs. ITfil a.fV in the general case they 
are quantitatively and even qualitatively different [Figs. 12.61 b,c,d,e — these cases 
correspond to the most nontrivial Tc(dp) behavior]. Thus to obtain reliable results 
one should use one of the exact (multimode or fundamental-solution) techniques. 

Spatial dependence of the order parameter 

According to the general theory of second-order phase transitions, the nonzero order 
parameter characterizes the "ordered" phase, in particular, the superconducting state. 
In the BCS theory, the role of the order parameter is usually played by A. However, 
this choice of the order parameter is not a good one when we consider proximity 
systems, in which a superconductor contacts a normal metal or a ferromagnet. The 
order parameter A, defined as 



is zero in the nonsuperconducting part of the structure, if the pairing constant A is 
zero. Nevertheless, the superconducting correlations are induced into the nonsuper- 
conducting metal, and these correlations are adequately described not by A but by 
the anomalous Green function F . 

Thus it is reasonable to define also a "proximity order parameter": 



where r denotes the imaginary time [in the S metal F{x,t = 0) cx A(a;)]. This 
function is real due to the symmetry relation F(— w„) = F*(w„). The proximity 
effect in the SF bilayer is characterized by the spatial behavior of this proximity 
order parameter. 

We illustrate this dependence in Fig. 12.71 which shows two cases differing by the 
thickness of the F layer dp (and by the corresponding Tc). Although the critical 
temperatures differ by more than the order of magnitude, the normalized proximity 
order parameters are very close to each other, which means that the value of Tc has 
almost no effect on the shape of F{x, t — 0). Details of the calculation are presented 
in Appendix IdI 

Another feature seen from Fig. 12. 71 is that the proximity order parameter in the F 
layer changes its sign when the thickness of the F layer increases (this feature can be 
seen for the dotted curve, although negative values of the order parameter have very 
small amplitudes). We discuss this oscillating behavior in the next section. 

2.1.6 Discussion 

Qualitative explanation of the nonmonotonic Tc{dp) behavior 

A qualitative explanation of the nonmonotonic Tc{dp) behavior in SF bilayers has 
been presented in several papers (see, e.g., Ref. 51]). Below we present another 
interpretation. 




(2.34) 




(2.35) 
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Figure 2.7: Spatial dependence of the proximity order parameter normalized by 
its value at the outer surface of the S layer. Two cases are shown differing by the 
thickness of the F layer dp (and by the corresponding Tc) at jb — 0.05. The other 
parameters are the same as in Fig. 12.41 where the chosen cases are indicated by the 
bold points. Although the critical temperatures differ by more than the order of 
magnitude, the normalized proximity order parameters are very close to each other, 
which means that the value of Tc has almost no effect on the shape of F(a;, r = 0). 
The jump at the SF interface is due to its finite resistance. With an increase of dp 
the proximity order parameter starts to oscillate, changing its sign (this can be seen 
for the dotted curve, although negative values of the proximity order parameter have 
very small amplitudes). 
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Figure 2.8: Four types of trajectories contributing (in the sense of Feynman's path 
integral) to the anomalous wave function of correlated quasiparticles in the ferromag- 
netic region. The solid lines correspond to electrons, the dashed lines — to holes; the 
arrows indicate the direction of the velocity. 



The thickness of the F layer at which the minimum of Tc{dp) occurs, can be 
estimated from qualitative arguments based on the interference of quasiparticles in 
the ferromagnet. Let us consider a point x inside the F layer. According to Feynman's 
interpretation of quantum mechanics ||52], the quasiparticle wave function may be 
represented as a sum of wave amplitudes over all classical trajectories; the wave 
amplitude for a given trajectory is equal to exp(iiS'), where S is the classical action 
along this trajectory. We are interested in an anomalous wave function of correlated 
quasiparticles, which characterizes superconductivity; this function is equivalent to 
the anomalous Green function F{x). To obtain this wave function we must sum 
over trajectories that (i) start and end at the point x, (ii) change the type of the 
quasiparticle (i.e., convert an electron into a hole, or vice versa). There are four 
kinds of trajectories that should be taken into account (see Fig. I2.8|l . Two of them 
(denoted 1 and 2) start in the direction toward the SF interface (as an electron and 
as a hole), experience the Andreev reflection, and return to the point x. The other 
two trajectories (denoted 3 and 4) start in the direction away from the interface, 
experience normal reflection at the outer surface of the F layer, move toward the 
SF interface, experience the Andreev reflection there, and finally return to the point 
X. The main contribution is given by the trajectories normal to the interface. The 
corresponding actions are 

Si = -Qx - a, (2.36) 
S2^Qx- a, (2.37) 
5-3 = ~Q{2dF + x)-a, (2.38) 
S4 ^ Q{2dF + x) - a (2.39) 

(note that x < 0), where Q is the difference between the wave numbers of the 
electron and the hole, and a = arccos(i?/A) is the phase of the Andreev reflec- 
tion. To make our arguments more clear, we assume that the ferromagnet is strong. 
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the SF interface is ideal, and consider the clean limit first: then Q ^ ke — kh = 
^2m(E + /i + E-p) — ^2m{—E — h + Ep) « 2h/v, where E is the quasiparticle en- 
ergy, Ep is the Fermi energy, and v is the Fermi velocity. Thus the anomalous wave 
function of the quasiparticles is 

4 

F{x) (X ^ exp(iS'„) cx cos{QdF) cos (Qldp + x]) . (2.40) 

n=l 

The suppression of Tc by the ferromagnet is determined by the value of the wave 
function at the SF interface: F{0) cx cos^ {Qdp)- The minimum of Tc corresponds to 
the minimum value of F{0) which is achieved at dp = tt/2Q. In the dirty limit the 
above expression for Q is replaced by 

nr 27r 

(here we have defined the wavelength of the oscillations Xh); hence the minimum of 
Tcidp) takes place at 

_,(min) _ 7[ — hh. 

2y h ^ A 

For the bilayer of Ref. Jjj we obtain d^'"' sa 7 nm, whereas the experimental value 
is 5nm (Fig.|2i3|l; thus our qualitative estimate is reasonable. 

The arguments given above seem to yield not only the minimum but rather a 
succession of minima and maxima. However, numerically we obtain either a single 
minimum or a minimum followed by a weak maximum fFig. 12.41 . The reason for this 
is that actually the anomalous wave function not only oscillates in the ferromagnetic 
layer but also decays exponentially, which makes the amplitude of the subsequent 
oscillations almost invisible. 

Finally, we note that our arguments concerning oscillations of F{x) also apply to 
a half-infinite ferromagnet, where we should take into account only the trajectories 1 
and 2 (see Fig. I2.8|l . This yields F{x) cx coa{Qx) (another qualitative explanation of 
this result can be found, for example, in Ref. ^^)- 



Multilayered structures 

The methods developed and the results obtained in this chapter apply directly to 
more complicated symmetric multilayered structures in the 0-state such as SFS and 
FSF trilayers, SFIFS and FSISF systems (I denotes an arbitrary potential barrier), 
and SF superlattices. In such systems an SF bilayer can be considered as a unit cell, 
and joining together the solutions of the Usadel equations in each bilayer we obtain 
the solution for the whole system (for more details see Sec. ll.7|l . 

Our methods can be generalized to take account of possible superconductive 
and/or magnetic 7r-states (when A and/or h may change their signs from layer to 
layer). In this case the system cannot be equivalently separated into a set of bilay- 
ers. Mathematically, this means that the solutions of the Usadel equations lose their 
purely cosine form [see Eqs. (|22||, lITTiji . (I2l5|l . ll2l9jl . ^1^, (I2.30b|l ] acquiring a 
sine part as well. 
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2.1.7 Conclusions 

In the present section, we have developed two methods for calculating the critical 
temperature of a SF bilayer as a function of its parameters (the thicknesses and ma- 
terial parameters of the layers, the quality of the interface). The multimode method is 
a generalization of the corresponding approach developed in Ref. [201 for SN systems. 
However, the rigorous justification of this method is not clear. Therefore, we propose 
yet another approach — the method of fundamental solution, which is mathematically 
rigorous. The results demonstrate that the two methods are equivalent; however, at 
low temperatures (compared to Tcs) the accuracy requirements are stricter for the 
multimode method, and the method of fundamental solution is preferable. Comparing 
our method with experiment we obtain good agreement. 

In the general case, we confirm three characteristic types of the Tc{dF) behavior, 
found previously in limiting cases or by approximate methods: 1) nonmonotonic decay 
of Tc to a finite value exhibiting a minimum at particular dp, 2) reentrant behavior, 
characterized by vanishing of Tc in a certain interval oidp and finite values otherwise, 
3) monotonic decay of Tc and vanishing at finite dp. Qualitatively, the nonmonotonic 
behavior of Tddp) is explained by interference of quasiparticles in the F layer, which 
can be either constructive or destructive depending on the value of dp- 

Using the developed methods we have checked the accuracy of the widely used 
single-mode approximation. We conclude that although at some parameters the re- 
sults of the single-mode and exact methods are close, in the general case they are 
quantitatively and even qualitatively different. Thus, to obtain reliable results one 
should use one of the exact (multimode or fundamental-solution) techniques. 

The spatial dependence of the order parameter (at the transition point) is shown 
to be almost insensitive to the value of Tc. 

The methods developed and the results obtained in this section, apply directly to 
more complicated symmetric multilayered structures in the 0-state such as SFS and 
FSF trilayers, SFIFS and FSISF systems, and SF superlattices. Our methods can 
be generalized to take account of possible superconductive and / or magnetic 7r-states 
(when A and/or h may change their signs from layer to layer). 

In several limiting cases, Tc is considered analytically. 
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2.2 Triplet proximity effect in FSF trilayers 

2.2.1 Introduction 

A striking feature of the proximity effect between singlet superconductors and inho- 
mogeneous ferromagnets is the possibihty of generating the triplet superconducting 
component [27, 28j. Recently, it was shown that the triplet component also arises 
in the case of several homogeneous but differently oriented ferromagnets jSH] • Physi- 
cally, the generating of the triplet component in SF systems |27l 1281 1^ is similar to 
the case of magnetic superconductors I^Oj . 

In Ref. [22j, the Josephson effect was studied having in mind that the supercon- 
ductivity in the system is not destroyed by the ferromagnets. However, this issue 
requires separate study. 

Although the SF proximity effect is rather well studied, the influence of the mutual 
orientation of F layers magnetizations (exchange fields) on Tc of layered SF structures 
has been mostly considered basing on the cases of parallel (P) and antiparallel (AP) 
alignment [H \lMi IH^ \XM \l^4^ \'^^ . At the same time, those are the only cases when the 
triplet component is absent. 

A FSF trilayer with homogeneous but noncollinear magnetizations of the F layers 
is the simplest example of a layered structure in which the triplet component is gen- 
erated. The triplet component (superconducting correlations between quasiparticles 
with parallel spins) arises due to a mechanism that is similar to the one described in 
Ref. j2Hl; with the difference that instead of local magnetic inhomogeneity we deal 
with magnetic inhomogeneity of the structure as a whole. This mechanism can be 
described in terms of the Andreev reflection at FS interfaces. In the case of a single 
FS interface, the Andreev reflection of a spin-polarized electron impinging on the FS 
interface from the F side, generates the singlet superconducting correlations in the 
ferromagnet. At the same time, in the FSF trilayer, we must take into account that 
the Andreev reflection from the FS interface is nonlocal in space: it takes place on a 
scale of the order of the coherence length in the vicinity of the interface. If the S layer 
is thin, this process "touches" the second F layer. If its magnetization is noncoUinear 
with the first one, it acts as a spin splitter, inducing the opposite spin component and 
hence the triplet superconducting correlations between the input and the output of 
the Andreev reflection. 

The critical temperature of the noncollinear FSF system was studied in Ref. 
However, in that work the triplet component was not taken into account. Thus 
calculation of Tc in the noncollinear FSF trilayer is still an open question. 

In this chapter we study the critical temperature of a FSF trilayer at arbitrary 
angle between the in-plane magnetizations (see Fig. I2.9|l . which makes it necessary 
to take the triplet component into account. 

2.2.2 General description 

We consider the dirty limit, which is described by the Usadel equations. Near Tc, 
the Usadel equations are linearized and contain only the anomalous Green function 
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Figure 2.9: FSF trilayer. The system is the same as in Ref. PEl- The thickness of 
the S layer is 2ds, of each F layer — dp- The center of the S layer corresponds to 
X = Q. The thick arrows in the F layers denote the exchange fields h lying in the 
(y, z) plane. The angle between the in-plane exchange fields is 2a. 



F j22|: 

Dd'^F . ,s i 



2 dx^ '"^"'^ ' 2 

_ ( fn /tt 



Lj^\F + Ad3-- sgn w„ ( FH* +HF)^0, (2.43) 
F 



fll /it 

Here D is the diffusion constant {Ds and Dp for the S and F layers), w„ = TTT{2n+l) 
are the Matsubara frequencies, and is the third Pauli matrix. The function F 
is a matrix in the spin space. The /f| and fn components describe the triplet 
superconducting correlations. In the P and AP cases it is sufficient to consider only 
the scalar equation for the singlet component /^j. 

Equation H2.43II is written in the general case when both order parameter and 
exchange field are present. In our system, in the F layers the order parameter is 
absent, A = 0, while 

if = /i ((72 sina + CT3 cosa) (2.44) 

at the exchange field h = ft,(0,sina,cosa). h is the exchange energy, and a describes 
the direction of the in-plane magnetization. 

In the S layer, the exchange energy is zero, while the superconducting order pa- 
rameter obeys the self-consistency equation 

= '^i;(i;^-/u), (2^45) 

where Tcs is the critical temperature of the S material. In the case of a single S layer, 
A can be chosen real. 

The boundary conditions at the outer surfaces of the trilayer are 



(2.46) 
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while at the SF interfaces 



dFs dFp psis .„^ 

^s-j— =iCf-j—, 7 = — ^, (2.47 

dx dx Pfc,f 

± ?F7B^ = i?5 - ^^F, 7B = (2.48) 

dx Pf^f 



Here Cs(f) = \/Ds{F)/^t^TcS and P5(_f) are the coherence lengths and the normal 
state resistivities of the S and F metals, Rb is the total resistance of the SF boundary, 
and A is its area. The ± sign in the l.h.s. of Eq. Ij2.48ll refers to the left and right SF 
interface, respectively. The above boundary conditions were derived for SN interfaces 
|2T] (N is a normal metal); their use in the SF case is justified by the small parameter 
h/Ep <C 1 (Ep is the Fermi energy). 

Our strategy is to reduce the problem to the S layer only, with effective boundary 
conditions. 

We expand the Green function F in the basis of the Pauli matrices txi , (T2 , f?3 , 
and the unit matrix ctq- It can be shown that the solution has the form 

F = /oao + fidi + hds. (2.49) 

The /o component is imaginary, while /i and /a are real. The relations /o(— ^n) = 
-/o(wn), /i(-Wn) = -/i(w„), /3(-Wn) = fsi^^n) make it sufficient to consider only 
positive Matsubara frequencies. 

The /i component describes a special type of triplet condensate |27l I2ij) . odd 
in frequency a;„) = — /i(w„)] and even in the relative momentum of electrons 
in the Cooper pair,^ which is similar to the one proposed by Berezinskii It 
is independence on the momentum direction that allows the triplet condensate to 
survive in the diffusive limit, in contrast to the standard (odd in momentum) case 
|88j . This Berezinskii-type triplet phase is characterized by zero orbital momentum 
of the triplet Cooper pairs. The odd dependence on frequency implies (spontaneous) 
breaking of the time-reversal symmetry. 

Equation (I2.43|l yields three coupled scalar equations (we consider cj„ > 0): 

w„/o - ihh cos a = 0, 



2 dx^ 
Dd^fi 



^n/i+/i/3sina = 0, (2.50) 



Dd^h 



2 dx^ 

- ^nf?, — ihfQ cos a — hfi sin a + A = 0. 



2 da;2 

Analyzing symmetries implied by Eqs. Il2.5fl|l and geometry of the system, we conclude 
that /o(a;) = fo{—x), fi{x) — — x), /3(a;) = /^{—x). Thus we can consider only 
one half of the system, say a; < 0, while the boundary conditions at x = are 

f =0, /i=0, f = 0. (2.51) 



^The relative momentum of electrons in the Cooper pair should not be confused with the mo- 
mentum corresponding to the x coordinate. The latter describes the center of mass of the Cooper 
pair, while the relative momentum is not an argument of the quasiclassical Green functions which 
are integrated over it. 
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Below we shall use the following wave vectors: 

kp = yJlu^nlDp-, kh = \/h/DF, 



kp + 2ikf^, 



ks = y/2LJn/ Ds- 

The solution in the left F layer, satisfying the boundary condition (|2 
form 



Fp = Ci {iaa sin a + ai cosq) cosh [kp (x + ds + dp)] 
+C2 ((To cos a + i(Ti sin a + (T3) cosh kh (x + ds + dp) 
+C3 (ao cos a + iai sin a — (T3) cosh kh (x + ds + dp] 



(2.52) 
has the 



(2.53) 



The matrix boundary condition II2.48|I yields three scalar equations, which allow to 
express the coefficients C'l, C2, C3 in terms of the components /o , /i , /a of the Green 
function on the S side of the FS interface: 



C3 



—ifo sin a + fi cos a 

1 + 7b^f 
/o cos a ~ ifi sin a + fa 

2{l+jBAh) 
/o cos a - i/isina - /a 



(2.54) 



2(1 + 75^;;) 

where we have introduced the following notations: 

Ap = kp^p tanh^kpdp), Ah ~ khS,p tanh^khdp), 
jAp -/Ah 



1 + ■jbAp 

Then the boundary condition (I2.47|l yields 
dfo 



Vh = 



1 + jeAh ■ 



(2.55) 



^S^ = fo{Vp sin^ a + ReVh cos^ a) 



if I {V — Re Vh) sina cos a + i/a Im Vh cos a, 



^ = i/o ( Vf — Re Vh) sin a cos «+ 
ax 



+ /i (l^F cos^ a + Re Vh sin^ a) + /a Im T4 sin a, 



dfi 



£.8-7- = i/o Im cos a - /iImT4sina + /aReV/i. 
ax 



(2.56) 

(2.57) 
(2.58) 



Thus the Green function of the F layer is eliminated, and we obtain equations for the 
S layer only. Moreover, we can proceed further, because in the S layer the unknown 
function A(x) (which must be determined self-consistently) only enters the equation 
for the /a component [the last of Eqs. H2.5flll ]. At the same time, taking boundary 
conditions Ij2.51ll into account, we can write /o = cosh(fc5x), /i = i3i sinh(fc5x). 
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Excluding i?o and Bi from the boundary conditions H2.56ll - H2.58|l . we arrive at the 
effective boundary condition for f^: 

^s^^^h> (2.59) 

where 

>v=RcF.+ ^ (2.60) 

ks^sA{a) + Re Vh 
and the angular dependence is determined by 



ks^s tanh(fc5ds) + Vp [sin^ a + tanh^ (ksds) cos^ a] 
ksis [cos^ a + tanh^ {ksds) sin^ a] + Vp tanh(fc5d5) 



Effectively, we obtain the following problem: 



Aln^ = 27rr ^ (—-/a), (2.62) 

Ds dy-. 



2 dx^ 



LUrJs + A = 0, (2.63) 



^s'^^=y^Mhi^dsl M = (2.64) 

All information about the F layers is contained in a single function W, all informa- 
tion about the misorientation angle — in its part A{a). Knowledge of W is already 
sufficient to draw several general conclusions about the behavior of Tc. First, if the 
S layer is thick, i.e. ds ^ fs, then tanh(fcsc?s) ~ 1 at characteristic frequencies, 
and Tc does not depend on a. Qualitatively, this happens because the effect of mu- 
tual orientation of the F layers is due to "interaction" between the two SF interfaces, 
which is efficient only in the case of thin S layer. Second, Tc does not depend on dp if 
dp ^ ^F- Qualitatively, this is due to the fact that the superconducting correlations 
penetrate from the S to F layer only on the scale ^p, and if dp is even larger, it does 
not affect Tc. 

The triplet component is "nonmonotonic" as a function of a: it vanishes at a = 
and a = 7r/2 (P and AP case, respectively), and arises only between the two boundary 
values. However, the Tc{a) dependence is always monotonic (Tc monotonically grows 
with increasing a). It can be directly proven from the monotonic behavior of A{a), 
and, hence, W. This rigorously derived conclusion disproves the result obtained by 
the approximate single- mode method in Ref ,32|, where it was claimed that Tc in the 
AP configuration can be smaller than in the P case. 



2.2.3 Numerical results in the general case 

At arbitrary parameters of the system, the effective problem (I2.62ll - H2.64|l can be 
solved numerically by the methods developed in Sees. 12.1.31 [2?OI 

Numerical results obtained by the exact methods are shown in Figs. 12.101 12.111 
A question arises: why is there pronounced angular dependence in the case ds > ^s, 
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Figure 2.10: Critical temperature Tc vs. thickness of the F layers dp, which is normal- 
ized on the wavelength of the singlet component oscillations — 2n/kh. Parameters 
ds/^S — 1-2, h/nTcs = 6.8, 7 = 0.15, 7b = 0.02 are close to the experiment {llj . 
The curves are calculated at different angles 2a between the in-plane exchange fields 
in the F layers. A minimum of the Tcidp) dependence and reentrant behavior were 
obtained for coUinear orientations in Refs. [3l [T5l[Tfi| . 

when the S layer is not thin? The answer is that the condition ds ^ = 
Ds/2ttTcs is a sufficient condition of thin S layer, whereas the necessary condition 
is weaker: c?s ^ = \/Ds/2TrTc, since the characteristic energy for a particular 
system is nTc with its own value of Tc. The two conditions become essentially differ- 
ent if Tc is notably suppressed, and in this case Tc can exhibit pronounced angular 
dependence at ds <^ ^, while it is possible to have ds > S,s- 

Experimentally, the conditions for observing the angular dependence of Tc are 
more easily met when Tc is essentially (but not completely) suppressed. Accordingly, 
the effect of a on Tcidp) dependence is most pronounced near the reentrant behavior. 
Experimental detection of the reentrant behavior was reported in Ref . . 



2.2.4 Analytical results for the case of thin S layer 

If ds ^ ^s, then A is constant. The Usadel equation II2.6^^|I can be solved, and the 
equation determining Tc takes the form 



ln^ = 2.r.5: 

Tc ^„ V UJn 



UJ„>0 



■ WnTcsCs/ds 



where W is given by Eq. Il2.fi0|l with simplified function A{a): 

kg^sds + Vp [sin^ a + (ksds)"^ cos^ a] 



A = 



ks^s [cos^ a + (ksds)^ sin^ a] + Vpksds 



(2.65) 



(2.66) 
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Figure 2.11: Tc vs. misorientation angle 2a. The curves correspond to different 
thicknesses of the F layers dp- The parameters are the same as in Fig. l2.1fll 



For the P and AP alignments, under additional assumption of strong ferromag- 
netism {h 3> ttTcs), we obtain: 



1 >ve.ic.\ , A 



where -0 is the digamma function, Vh is determined by Eqs. II2.55|I with kh = {l+i)kfi, 
and in the region of parameters, where Tc ^ [the corresponding conditions can be 
extracted from the results for the critical thickness — see Eqs. Ij2.7()ll . Ij2.71ll below], 
we may write 

W ^ReVh + ^ilmVhf . (2.69) 

Due to symmetry, the result for the P case H2.67|l reproduces that for the SF bilayer 
with S layer of thickness ds [Eq. ljB.2ll ]. In the AP case, if the second terms in the 
r.h.s. of Eq. Ij2.(i9|l can be neglected (e.g., at khdp ^ 1 in the region of parameters 
where Tc ^ 0), then W = ReT4 and we reproduce the result of Ref. PSI. However, 
the second term becomes essential in the Cooper limit, defined by conditions ds ^ 
Ds/^ujd, dp <C vtva\{yj Dp /2ujd^ kf^^), = 0, with ujd the Debye energy of the 
S material. In this case ReV/i — and Eqs. H2.68II . Il2.fi9|l reproduce the result of 
Tagirov g]. 

The critical thickness d^s of the S layer, below which the superconductivity van- 
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ishes, immediately follows from Eqs. Il2.fi7|l . H2.68|l for the P and AP cases:^ 

= 2exp(C) \Vh\ , -f- = 2exp(C) W (2.70) 



is " " is 



at 



dcs 

is 



< 1. (2.71) 



Here C « 0.577 is Euler's constant. Condition H2.71II is necessary for applicability 
of Eqs. H2.7flll . If this condition is not satisfied, then Eqs. only tell us that 

at ds/is <C 1 the superconductivity is certainly absent, i.e., Tc = 0. According to 
the monotonic growth of Tc(a), the function dcsio) decreases monotonically, hence 
dcS > '^fs- 7b = 0, khdp » 1, Eqs. H2.7()ll reproduce the results of Ref. PB] for 
the P and AP cases. 

The Tc{a) dependence can be most easily studied in the Cooper limit. In this case 
a simple analysis can be done already on the level of the Usadel equations, and the 
system is described as a uniform layer with the effective exchange energy^ 

Tp 

hcs — — hcosa, (2.72) 

TS 

where T5(^) = 2ds^p-^RBA/ Ps(f)Ds(f) (similarly to Eqs. The accuracy of 

this result is limited to the first order over h, which becomes insufficient in the vicinity 
ofa = 7r/2. AtQ; = 7r/2, the first-order effect of h vanishes, while a more accurate 
analysis (Ref. and Eqs. H2.68II . (I2.69|l ) reveals the second-order effect of h on Tc- 
Let us now consider the same limit as in Ref. |i36) : 

ds^is, fc/idF>l, /i>7rTcs, 7s = 0, (2.73) 
ikhiF^ < 1. (2.74) 

is 

The condition to have superconductivity at least at some orientations has the form 
^cs" ds is, and in the case under discussion, Eqs. Ij2.7()ll . (12.7111 yield: 

2exp(C)7fc^eF<^«l, (2.75) 
is 

hence condition (I2.74|l becomes redundant. 

Starting from Eqs. Il2.65|l . 112.6011 . Ij2.66|l . we finally obtain the following equation 
for T,: 



Tr V 2 2TrTr / V 2 27rr, / V 2 



^When calculating dcS, we assume that the phase transition is of the second order. However, this 
issue may require a separate study. In principle, under some circumstances, the order of the phase 
transition in SF systems can change from second to first one — see Sec. 12.1.51 

^ Since aj„ was neglected in comparison with h in the Usadel equation, the result of the Cooper 
limit is valid only at rg Tp. 
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Figure 2.12: Critical thickness of the S layer dcs vs. misorientation angle 2a. Dashed 
line is the result of Baladie et al. obtained without account of the triplet com- 
ponent. 



where 

R = l-Q, 



1 ^ sin^ a 

2 2\/ s'm^ a — A cos^ a 
do 



111 2 = -t-ttTcs f 1 + cos^ a ± Vsin-^ a- A cos^ a) , (2.77) 
da = jkh£,F£.s/2. 

In the P and AP cases, where the triplet component is absent, Eqs. H2.76II . 112.7711 
reproduce the results of Refs. jHIlEEl- At the same time, at a noncollinear alignment 
the results are clearly different. 

The critical thickness is found from Eqs. H2.7(ill . II2.77|I : 



dcs{a) . /- / 1 + cos^ a+ \J sin** a - 4 cos^ a \ 2\/ sin* q-4 cos^ a 

— = 4v2exp(6) cos a -j^^^=^^^= 

do \ 1 + cos^ a — v sin^ a — 4 cos^ a / 

(2.78) 

Although the square root in this expression can become imaginary, the whole ex- 
pression remains real (z* is real if \z\ = 1). Figure [2. 121 illustrates the result Ij2.78|l . 



2.2.5 Conditions of existence of the odd triplet superconduc- 
tivity in SF multilayers 

Now we turn to analyze the conditions of appHcability for the results reported in 
Ref. jSni- A noncollinear FSF trilayer is a unit cell of the multilayered structure 
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studied in that work. The main result of Ref. [22l, the Josephson current due to the 
long-range triplet component, requires that the S layer is thin ds ^ £,s, while the F 
layers are thick for the singlet component and moderate for the triplet one: fc^^ <C 
< dp In this case the condition that superconductivity is not completely 

suppressed at least in the vicinity of the AP alignment [Eqs. Ij2.70p . Il2.71|l ] takes the 
form 

4exp(C)7fc.e.7T^f^|2^^ < ^ « 1. (2.79) 
At 7_B = (as it was assumed in Ref. [22|), this yields 

2 exp(C) -fkh^F < ^ < 1, (2.80) 

which is a rather strong condition for 7, since khS,F 3> 1. Finite interface resistance 
relaxes this condition: already at 75 > 1, Eq. (I2.79|l yields 

2exp(C)^ < ^< 1. (2.81) 

IB ?S 

The condition that superconductivity exists at all orientations has the form similar 
to Eq. II2.79|I but with the corresponding expression for d^g instead of d^f in the 
l.h.s. This only leads to a minor difference, since the two critical thicknesses are of 
the same order: d^g = ^/2d^g at 7s — 0, while c?^ = d^g at 73 > 1. 



2.2.6 Conclusions 

In this section, we have studied Tc of a FSF trilayer as a function of its parameters, 
in particular, the angle 2a between magnetizations of the F layers. At noncoUinear 
orientation of magnetizations, we take into account the triplet superconducting com- 
ponent generated in the system due to nonlocality of the Andreev reflection. 

The Tc{a) dependence becomes pronounced when the S layer is thin, and can lead 
to switching between superconducting and non-superconducting states as the angle 
is varied. In the general case, we reduce the problem to the form that allows us to 
employ the exact numerical methods of Sec. 12.11 In the most interesting limiting 
cases, we analytically analyze and the critical thickness of the S layer. 

Our results directly apply to multilayered SF structures, where a FSF trilayer 
is a unit cell. We have formulated the conditions which are necessary for existence 
of recently investigated odd triplet superconductivity in SF multilayers [22l. The 
conditions are severe in the case of transparent interfaces, but become reaHzable at 
moderate interface transparency. 



76 



Chapter 2 



References 

[1] V. V. Ryazanov, V. A. Oboznov, A. Yu. Rusanov, A. V. Veretennikov, A. A. Gol- 
ubov, J. Aarts, Phys. Rev. Lett. 86, 2427 (2001); 

V. V. Ryazanov, V. A. Oboznov, A.V. Veretennikov, A.Yu. Rusanov, Phys. 

Rev. B 65, 020501(R) (2001). 
[2] T. Kontos, M. Aprili, J. Lesueur, X. Grison, Phys. Rev. Lett. 86, 304 (2001); 

T. Kontos, M. Aprih, J. Lesueur, F. Genet, B. Stephanidis, R. Boursier, Phys. 

Rev. Lett. 89, 137007 (2002). 
[3] Z. Radovic, M. Ledvij, Lj. Dobrosavljevic-Grujic, A. L Buzdin, J. R. Clem, Phys. 

Rev. B 44, 759 (1991). 
[4] L.R. Tagirov, Phys. Rev. Lett. 83, 2058 (1999). 
[5] A. Buzdin, Phys. Rev. B 62, 11377 (2000). 

[6] M. Zareyan, W. Belzig, Yu.V. Nazarov, Phys. Rev. Lett. 86, 308 (2001). 
[7] J.S. Jiang, D. Davidovic, D.H. Reich, C.L. Chien, Phys. Rev. Lett. 74, 314 
(1995). 

[8] Th. Miihge, N. N. Garif'yanov, Yu. V. Goryunov, G. G. KhaliuUin, L. R. Tagirov, 
K. Westerholt, LA. Garifullin, H. Zabel, Phys. Rev. Lett. 77, 1857 (1996). 

[9] J. Aarts, J.M.E. Geers, E. Briick, A. A. Golubov, R. Coehoorn, Phys. Rev. B 
56, 2779 (1997). 

[10] L. Lazar, K. Westerholt, H. Zabel, L. R. Tagirov, Yu. V. Goryunov, 
N.N. Garif'yanov, LA. Garifullin, Phys. Rev. B 61, 3711 (2000). 

[11] V. V. Ryazanov, V. A. Oboznov, A. S. Prokofiev, S. V. Dubonos, Pis'ma Zh. Eksp. 
Teor. Fiz. 77, 43 (2003) [JETP Lett. 77, 39 (2003)]. 

[12] A. Rusanov, R. Boogaard, M. Hesselberth, H. Sellier, J. Aarts, Physica C 369, 
300 (2002). 

[13] A.L Buzdin, B. Vujicic, M. Yu. Kupriyanov, Zh. Eksp. Teor. Fiz. 101, 231 

(1992) [Sov. Phys. JETP 74, 124 (1992)]. 
[14] E. A. Demler, G.B. Arnold, M.R. Beasley, Phys. Rev. B 55, 15174 (1997). 
[15] Yu.N. Proshin, M. G. Khusainov, Zh. Eksp. Teor. Fiz. 113, 1708 (1998); 116, 

1887 (1999) [JETP 86, 930 (1998); 89, 1021 (1999)]; 

M.G. Khusainov, Yu.N. Proshin, Phys. Rev. B 56, R14283 (1997); 62, 6832 
(2000). 

[16] L.R. Tagirov, Physica C 307, 145 (1998). 
[17] K.D. Usadel, Phys. Rev. Lett. 25, 507 (1970). 

[18] A. L Larkin, Yu. N. Ovchinnikov, in Nonequilihrium Superconductivity, edited 
by D.N. Langenberg and A.L Larkin (Elsevier, New York, 1986), p. 530, and 
references therein. 

[19] J. Rammer, H. Smith, Rev. Mod. Phys. 58, 323 (1986). 

[20] A. A. Golubov, M. Yu. Kupriyanov, V. F. Lukichev, A. A. Orlikovskii, Mikroelek- 
tronika 12, 355 (1983) [Sov. J. Microelectronics 12, 191 (1984)]. 

[21] M.Yu. Kupriyanov, V. F. Lukichev, Zh. Eksp. Teor. Fiz. 94, 139 (1988) [Sov. 
Phys. JETP 67, 1163 (1988)]. 



Proximity effect in SF systems 



77 



[22] M. Abramowitz, I. A. Stegun, Handbook of Mathematical Functions (Dover, New 
York, 1974). 

[23] P. M. Morse, H. Feshbach, Methods of Theoretical Physics (McGraw-Hill, New 

York, 1953), Vol. 1. 
[24] V.V. Ryazanov, private communication. 

[25] G. Sarma, J. Phys. Chem. Solids 24, 1029 (1963); see also D. Saint-.James, 
G. Sarma, E.J. Thomas, Type II Superconductivity (Pergamon, Oxford, 1969), 
p. 159. 

[26] R. P. Feynman, A. R. Hibbs, Quantum Mechanics and Path Integrals (McGraw- 
Hill, New York, 1965). 

[27] F.S. Bergeret, K.B. Efetov, A.I. Larkin, Phys. Rev. B 62, 11872 (2000); 
F.S. Bergeret, A.F. Volkov, K.B. Efetov, Phys. Rev. Lett. 86, 4096 (2001). 

[28] A. Kadigrobov, R.I. Shekhter, M. Jonson, Europhys. Lett. 54, 394 (2001); Fiz. 
Nizk. Temp. 27, 1030 (2001) [Low Temp. Phys. 27, 760 (2001)]. 

[29] A.F. Volkov, F.S. Bergeret, K.B. Efetov, Phys. Rev. Lett. 90, 117006 (2003). 

[30] L.N. Bulaevskii, A.I. Rusinov, M. Kulic, J. Low Temp. Phys. 39, 255 (1980). 

[31] A.I. Buzdin, A.V. Vedyayev, N. V. Ryzhanova, Europhys. Lett. 48, 686 (1999). 

[32] M. G. Khusainov, Yu. A. Izyumov, Yu. N. Proshin, Pis'ma Zh. Eksp. Teor. Fiz. 
73, 386 (2001) [JETP Lett. 73, 344 (2001)]. 

[33] I. Baladie, A. Buzdin, Phys. Rev. B 67, 014523 (2003). 

[34] G. Deutscher, F. Meunier, Phys. Rev. Lett. 22, 395 (1969). 

[35] J.Y. Gu, C.-Y. You, J. S. Jiang, J. Pearson, Ya.B. Bazaliy, S.D. Bader, Phys. 
Rev. Lett. 89, 267001 (2002). 

[36] I. Baladie, A. Buzdin, N. Ryzhanova, A. Vedyayev, Phys. Rev. B 63, 054518 
(2001). 

[37] V. L. Berezinskii, Pis'ma Zh. Eksp. Teor. Fiz. 20, 628 (1974) [JETP Lett. 20, 
287 (1974)]. 

[38] A.I. Larkin, Pis'ma Zh. Eksp. Teor. Fiz. 2, 205 (1965) [Sov. Phys. JETP Lett. 
2, 130 (1965)]. 

[39] I. A. GarifuUin, D. A. Tikhonov, N. N. Garif'yanov, L. Lazar, Yu. V. Goryunov, 
S.Ya. Khlcbnikov, L. R. Tagirov, K. Westerholt, H. Zabel, Phys. Rev. B 66, 
020505(R) (2002). 



Chapter 3 

Josephson effect in SFS junctions 



3.1 Critical current in SFIFS junctions 
3.1.1 Introduction 

Josephson structures involving ferromagnets as weak link material are presently a 
subject of intensive study. The possibility of the so-called "7r-state", characterized 
by the negative sign of the critical current Ic, was theoretically predicted in SFS 
Josephson junctions [TIElEllllSlElCIIHlElinillllllllllSI- The first experimental 
observation of the crossover from 0- to 7r-state was reported by Ryazanov et al. 
|14j and explained in terms of temperature-dependent spatial oscillations of induced 
superconducting ordering in the diffusive F layer. 

More recently a number of new phenomena was predicted in junctions with more 
than one magnetically ordered layer. First, the crossover to the 7r-state was predicted 
in Ref. for the case of parallel magnetizations in SFIFS Josephson junctions even in 
the absence of the order parameter oscillations in thin F layers. Second, the possibility 
of the critical current enhancement by the exchange field in SFIFS junctions with 
thin F layers and antiparallel magnetization directions was discussed in the regimes 
of small S layer thicknesses pS] and bulk S electrodes pSl- However, in the case of 
thin S layers, Ref dealt with an idealized model in the tunneling limit, which lead 
to a divergency of the critical current at zero temperature; moreover, the physical 
mechanisms of the discovered enhancement effect were unclear. At the same time, 
for the bulk S case an approximate method was used in Refs. [G.TG] and a part of the 
results was obtained beyond its appHcability range (this will be shown below — see 
discussion of Fig. I3.2|l . 

The above intriguing scenario motivated us to attack the problem of the Joseph- 
son effect in SFIFS junctions by self-consistent solution of the Usadel equations for 
arbitrary thicknesses of the F layers and arbitrary barrier transparencies. Below we 
show that the O-tt transition in the case of parallel orientation of the exchange fields 
h or enhancement of Ic by h in the antiparallel case with thin F layers occur when the 
effective energy shift in the ferromagnets (due to the exchange field) becomes equal 
to a local value of the effective energy gap induced into the F layers. Under this 
condition a peak in the local density of states (DoS) near the SF interfaces is shifted 
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Figure 3.1: Types of SFS systems considered in the present chapter. 
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to zero energy. In the models with DoS of the BCS type this leads to a logarithmic di- 
vergency of Ic in the antiparallel case at zero temperature, similarly to the well known 
Riedel singularity of the ac supercurrent in SIS tunnel junctions at voltage eV = 2A. 
We also describe the general numerical method to solve the problem self-consistently 
and apply it for quantitative description of the O-tt transition and Ic enhancement in 
SFIFS junctions. 



3.1.2 Model 

We consider the structure of SFIFS type, where I is an insulating barrier of arbitrary 
strength (see Fig. 13. Ih ). We assume that the S layers are bulk and that the dirty 
limit conditions are fulfilled in the S and F metals. Although our method is applicable 
in the general situation of different ferromagnets and superconductors, for simplicity 
below we illustrate our results in the case when equivalent S and F materials are used 
on both sides of the structure (although the directions of the exchange field in the 
two F layers may be different), both F layers have the thickness dp, and the two 
SF interfaces have the same transparency. At the same time, we do not put any 
limitations on dp and the transparency. 

The Usadel functions G, F, F obey the normalization condition + FF = 1, 
which allows the following parameterizations in terms of the new function $: 



V/W2 _,.$(w„)**(-^^rO' 

F{LOn) 



F{UJ„) 
F*{-LOn). 



(3.1) 



The quantity w„ — ujn + ih corresponds to the general case when the exchange field 
h is present. However, in the S layers h — Q and we have simply uJn = ^n- 

We choose the x axis perpendicular to the plane of the interfaces with the origin 
at the barrier I. The Usadel equations (E] in the S and F layers have the form 



es 



ttTc d 
LOnGs dx 
ttT, d 



ijJnGp dx 



- -A, 



=0, 



(3.2) 
(3.3) 



where Tc is the critical temperature of the superconductors, A is the order parameter 
(which is nonzero only in the S layers), w„ is the Matsubara frequency, and the 
coherence lengths ^ are related to the diffusion constants D as Cs(f) = •\/-Ds(f)/27I'Tc. 
The order parameter satisfies the self-consistency equations 



Ain^ + .rV^^^^ii^^o. 



(3.4) 



We restrict ourselves to the cases of parallel and antiparallel orientations of the ex- 
change fields h in the ferromagnets. 
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The boundary conditions at the SF interfaces {x — Tdp) have the form (see 
Ref. ri9' for details) 



e^JL^.^.k^l^,, (3.5) 



dx ujn dx 

d 
< 

dx 



±7i. — -'fi^^GsU---), (3.6) 



RbA ps£.s 
with 7B = — — , 7 — — , 
PFt.F PfKf 

where Rb and A are the resistance and the area of the SF interfaces; ps{F) is the 
resistivity of the S (F) layer. At the I interface {x = 0) the boundary conditions read 

— TT^Fl — — ■7^^F2, (o.'J 

(jJnl OX UJn2 OX 

,b/-^^<^f. = GfJ^-^), (3.8) 

•ii RbjA 
with 7b J = — , 

PfS,f 

where the indices 1, 2 refer to the left and right hand side of the I interface, respec- 
tively. 

In the bulk of the S electrodes we assume a uniform current-carrying supercon- 
ducting state 

^ ^ Aoexp(»[T^/2-f2m..x]) 

where m is the electron's mass, Vs is the superfluid velocity, and ip is the phase 
difference across the junction. 

The supercurrent density is constant across the system. In the F part it is given 
by the expression 



2ep ^ 



(3.10) 



while analogous formula for the S part is obtained if we substitute a;„ cjn. This 
expression, together with the boundary condition Ij3.8ll and the symmetry relation 
F{—bJn, h) = F(ujn, —h), yields the formula for the supercurrent across the I interface: 

I = ^T.^^Mni{~hi)FF2{h2)] (3.11) 

[the functions F are related to $ via Eq. H3.1|l ]. 

3.1.3 Limit of thin F layers 



Let us consider the limit of thin F layers: c?f ^ niin {^F,VDF/2h). Under the 
condition jb/i ^ 1 we can neglect the suppression of superconductivity in the su- 
perconductors. We assume further that the transparency of the I barrier is small. 
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7s,/ ^ inax(l,7B), hence the SF bilayers are decoupled. In this case we can set 
Vs — and expand the solution of Eq. H3.3II in the F layers up to the second order 
in small spatial gradients. Applying the boundary condition 113. 6|l . we obtain the 
solution in the form similar to that in SN bilayer ^20. 



F2 



1 +7SJ\/Wnl,2/7rTcG's 

dp 



Aoexp(Tj<rf?/2), 



(3.12) 



with 7BM = IB 



6 



This solution corresponds to the 7m = limit of Refs. J2l lEl therefore the 
analytical results II3.13|I - II3.18|I presented below are explicitly or implicitly contained 
in those works. We employ them to discover the physical mechanisms leading to 
the effects of the O-tt transition and the enhancement of the critical current by the 
exchange field. 

Substituting Eq. H3.12II into the expression for the supercurrent l|3.11|l we obtain 
I{ip) — IcSiuip. For the parallel orientation of the exchange fields, hi = h2 = h, the 
critical current is 

A 



27rT ^ Al 



iJn>0 



while for the antiparallel orientation, hi 



^0 



^2 


+ 52' 




we obtain 




1 



Here we have introduced the following notations 
A 



52 



B = 



{iBAlh) + iBili^n (2 + IBM^) , 
2LL>nhjBM^ (1 + IBM^) , 



(3.13) 



(3.14) 



(3.15) 
(3.16) 



n 



TrTr 



At H/ttTc — 1/jBM and small w„, the expression under the sum in Eq. (13.1411 
behaves as l/w„, thus at low T the critical current diverges logarithmically: /c"'' « 
ln(rc/T). This effect was pointed out in Refs. [l5lfT6j . 

The above results become physically transparent in the real energy E represen- 
tation. Performing the analytical continuation in Eqs. 113. l|l . Ij3.12ll by replacement 
ujn —iE, we obtain the expression for the DoS per one spin projection (spin "up") 
vf{E) = Rc Gi;^(£;) in the F layers: 



lypiE) 



Re 



E 



E^ 



A2 

^0 



(3.17) 



E + jBMiE-h) 



ttTc 
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which demonstrates the energy renormaUzation due to the exchange field. Equation 
pj7| yields 

.^(0)^Re_Jg j^^/"^- , (3.18) 

\/ilBMh/TTTc)'^ - 1 

which shows that at h/nTc — 1/jbm the singularity in the DoS is shifted to the 
Fermi level. Exactly at this value of h/nTc the maximum of Ic""^ is achieved due to 
the overlap of two E~^/'^ singularities. This leads to logarithmic divergency of the 
critical current lj^H4|l in the limit T ^ 0, similarly to the well known Riedel singularity 
of nonstationary supercurrent in SIS tunnel junctions at voltage eV = 2Ao, where 
the energy shift is due to the electric potential. At the same value of the exchange 
field h/nTc = "^I^BM the critical current changes its sign (i.e., the crossover from 
to TT contact occurs) for parallel magnetizations of the F layers [see Eq. II3.13II ]. 

We emphasize that the scenario of the O-tt transition in our case differs from those 
studied before where the 7r-shift of the phase was either due to spatial oscillations of 
the order parameter in F layers (see, e.g., Ref. Ql] and references therein) or due to 
the proximity- induced phase rotation in S layers jH]. In our case the phase does not 
change in either layer; instead, it jumps at the SF interfaces. This scenario is most 
clearly illustrated in the Hmit of large h where Eqs. 113. l|l . II3.12|I yield Fp oc — lAsgn/i 
whereas Fs oc A; thus the phase jumps by 7r/2 at each of the SF interfaces, providing 
the total TT-shift between Fpi{—h) and Fp2{h) [it is the phase difference between these 
two functions that determines the supercurrent according to Eq. Ij3.11ll ]. Physically, 
these jumps are not a property of the SF interfaces as such. Indeed the condition 
for the O-TT transition is determined by the ^bm parameter, which depends on both 
the interface transparency and the thickness of the F layer. Probably, the 7r/2 phase 
is acquired due to multiple passage of quasiparticles, reflected at the SF interface, 
through the F layer (while only a very small phase is acquired at a single passage). 

The considered effects take place only for sufficiently low I-barrier transparency. 
Indeed, it follows from Eq. Ij3.12ll that Gpi^Jn) oc 1/ ^JujZ for small ujn under condition 
h/irTc — '^Mbm- As a result, the boundary condition H3.8|l yields that at 

< mm , (3.19) 



TtTc \dFlBJ IBJ 

the solutions Ij3.12ll are not valid, since in this frequency range the effective trans- 
parency of the I interface (the parameter GfiGf2/ib,i) increases and the spatial gra- 
dients in the F layers become large, similarly to the KuHk-Omelyanchuk case |2T1[22] . 
At these frequencies, the Green functions G, F (and hence the contribution to the 
critical current from these frequencies) are /i-independent. As a result, the barrier 
transparency parameter 7b./ provides the cutoff of the low-temperature logarithmic 

singularity of /i"'' at h/iTTc = IHbm [see Eq. (jlTTijl ]. According to Eq. l(3?T9|l . the 
critical current saturates at low temperature T* = Tcmiii{£,F /dp^Bj, 1b/ib,i)- We 
note that any asymmetry in the SFIFS junction will also lead to the cutoff of /c°' 
divergency. The above estimates are done for the case of low barrier transparency, 
S,F/dF7B,i <C 1 and Jb/jbj <C 1. The opposite regime of high transparency requires 
separate study. 
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Figure 3.2: Enhancement of the critical current (antiparallel magnetizations, sohd 
Hues) and the O-tt transition at which Ic changes its sign (parallel magnetizations, 
dashed lines) in the SFIFS junction at T/Tc = 0.05, jbm = 1, and jm = 0. Inset: 
the same for large values of 7m- Note that the limit dp/^p = is taken at fixed 
7bm = jBdp/S.F- Since it is jbm that determines the strength of the proximity 
effect in the limit of thin F layers, Ic demonstrates dependence on h. 



3.1.4 General case 

For arbitrary F-layer thicknesses and the interface parameters, we solved the bound- 
ary problem (j;^.l|l - (j;^.9|l numerically using an iterative procedure. Starting from trial 
values of the complex order parameters A and the Green functions Gs,f we solve 
the resulting Hnear equations and the boundary conditions for the functions ^s,f- 
After that we recalculate Gs,f and A. Then we repeat the iterations until conver- 
gency is reached. The self-consistency of calculations is checked by the condition of 
conservation of the supercurrent II3.10|I across the junction. We emphasize that our 
method is fully self-consistent: in particular, it includes the self-consistency over the 
superfluid velocity Vs, which is essential (contrary to the constriction case) in the 
quasi-one-dimensional geometry. 

Figure l3!2l shows Idh) dependencies calculated at T = 0.05 Tc from the numerical 
solution of the boundary problem for the fixed value of jbm = 1 and a 

set of different F-layers thicknesses and the SF interface parameters 7. The normal 
junction resistance is Rn = Rb,i + '^Rb + 2ppdF/A. The curves dp/^F = are the 
limits of vanishing dp/^F ratio at fixed 7bm and are calculated from Eqs. Ij3.i;^|l . 
Ij3.14ll . For thin F layers the results depend only on the combination jm — idp/iF- 
The enhancement of Ic and the crossover to the 7r-state are clearly seen for the 
antiparallel and parallel orientations, respectively. In accordance with the estimates 
given above, these effects take place for the values of the exchange field h close to 
ttTc. The enhancement disappears with increasing gradients in the F layers since the 
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Figure 3.3: Enhancement of the critical current (antiparallel magnetizations, sohd 
Hues) and the O-tt transition at which Ic changes its sign (parallel magnetizations, 
dashed lines) in the SFIFS junction: influence of temperature and barrier trans- 
parency. The dotted line corresponds to T/Tc = 0.01 and ^F/dp^Bj ~ 0; the 
parameters for other curves are given in the figure. 

solution Eq. Ij3.12ll loses its validity. This is illustrated in Fig. 13. 21 bv increasing the 
thickness dp or jm- 

Thus in the case of large the enhancement effect is absent, in contrast to the 
statement of Ref. [l^. This contradiction is due to the fact that the calculation in 
the large 7m limit in Ref. is valid only at T ^ Tc/j^.j, whereas the effect of 
enhancement exists only at small T. Therefore the enhancement effect at large 7m and 
small T in Ref. was obtained beyond the applicability range of the approximate 
method. At the same time, at small 7m, the Ic{h) dependencies in Refs. 
are qualitatively correct but quantitatively inaccurate because in the approximate 
analytical calculation the correction over small 7m was taken into account only in 
the Green function but not in the order parameter A. The accurate calculation 
requires to consider this correction self-consistently (similarly to Ref. [23|)- 

Influence of temperature and barrier transparency on the critical current anomaly 
is shown in Fig. 13.31 One can see that, in accordance with the above estimate, the 
cutoff of /c"'' singularity is provided by finite temperature or barrier transparency. 
Namely, with the decrease of the barrier strength parameter jb,i the peak magnitude 
starts to drop when the ratio dFjB,i/£,F becomes comparable to T/Tc- With further 
decrease of dp'^sj/^.F the singularity disappears, while the transition to the 7r-state 
shifts to large values of h. 

Figure ini riemonstrates the DoS in the F layers for one spin projection, calculated 
numerically in the limit of small I-barrier transparency. At ft. = we reproduce the 
well-known minigap existing in SN bilayer. At finite h the gap shifts in energy 
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Figure 3.4: Normalized density of states for spin "up" in the F layer for various 
exchange fields. 



(asymmetrically) and the peak in the DoS reaches zero energy at h/nTc = I/^bm- 
One can see that even for a small value 7m — 0.05 the peaks are rather broad, 
therefore the singularity in I'f'^ is suppressed by very rapidly. 

In the practically interesting limit of finite F-layer thickness (see Fig. 13.511 the 
numerical calculations demonstrate monotonic suppression of Ic with increase of the 
exchange field h for antiparallel magnetizations of the F layers and the O-tt crossover 
for the parallel case. One can see from Fig. 13.51 that for given temperature and 
thickness of the F layers it is possible to find the value of the exchange field at which 
switching between parallel and antiparallel orientations will lead to switching of Ic 
from nearly zero to a finite value (or to switching between and tt states) . 

The case of parallel F-layers magnetizations in the absence of the I barrier cor- 
responds to the standard SFS junction where the O-tt transition is possible due to 
spatial oscillations of induced superconducting ordering in the F layer. The ther- 
mally induced O-tt crossover in SFS junction was observed in Ref. ^J, where simple 
theory based on the linearized Usadel equations was also presented. Here we obtain 
such a crossover (see the inset in Fig. 13.511 from the fully self-consistent solution in 
the range of the exchange fields corresponding to that of Ref. Jlj . 

3.1.5 Conclusions 

In this section, we have developed a general method to solve the Usadel equations 
in SFIFS junctions self-consistently. Using our method, we have investigated theo- 
retically the Josephson current in SFIFS and SFS junctions. To clarify the physical 
mechanisms behind the observed effects, we analytically considered the limiting case 
of thin F layers. 
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Figure 3.5: Critical current in the general case: switching effect. T/Tc = 0.5, the 
solid and dashed Hnes correspond to the antiparallel and parallel orientations of mag- 
netizations, respectively. Inset: thermally induced O-tt crossover in the parallel case. 

If the magnetizations of the two F layers in the SFIFS junction are antiparallel, 
then the critical current can be enhanced by the exchange field, as it was predicted 
before in similar structures. We demonstrate that this effect is similar to the Riedel 
singularity that takes place in SIS junction due to shifting of the density of states by 
voltage; in the SFIFS junction, the shift is due to the exchange field. The logarithmic 
divergence of the maximal critical current is cut off by finite temperature or barrier 
transparency. 

If the magnetizations are parallel, then the junction can undergo the transition to 
the TT-state even in the limit of thin F layers. We demonstrate that this effect can be 
due to the effective 7r/2 shifts at the two SF interfaces. 
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3.2 Nonsinusoidal current-phase relation in SFcFS 
and SIFIS Josephson junctions 

3.2.1 Introduction 

The relation between the supercurrent / across a Josephson junction and the differ- 
ence (f between the phases of the order parameters in the superconducting banks is 
an important characteristic of the structure. The form of /(<^) dependence is essen- 
tially used for analyzing the dynamics of systems containing Josephson junctions 24J. 
Studying /((/?) also provides information on pairing symmetry in superconductors 25j. 

In structures with tunnel-type conductivity of a weak Hnk (SIS) the current- 
phase relation is sinusoidal, I{ip) = IcSin(p with Ic > 0, in the whole temperature 
range below the critical temperature. At the same time, in point contacts (ScS) 
and junctions with metallic type of conductivity (SNS) strong deviations from the 
sinusoidal form take place at low temperatures T [22] with the maximum of I{(p) 
achieved at tt/2 < ipmax < tt. 

The situation drastically changes if there is magnetoactive material in the region 
of weak link. The transition from the 0-state {Ic > 0) to the 7r-state {Ic < 0) has been 
theoretically predicted in a variety of Josephson structures containing ferromagnets 
PElEllilllElCIIHlElIiniiniliainiand experimentally observed in SFS and SIFS 
junctions [14ll2(i[[?7]. In the general case modifications of I{ip) do not reduce to the 
change of the critical current sign in the sinusoidal Josephson relation. It was shown 
that the presence of a ferromagnet may result in a nonsinusoidal shape of I{(p). In the 
diffusive regime, this situation occurs in long SFS junctions with ideally transparent 
interfaces [Hlini- However, in this case the effects take place only in a narrow interval 
of very low temperatures (due to smallness of the Thouless energy), while in the 
present section we shall consider short-length structures where the effects are more 
pronounced and exist practically in the whole temperature range. 

In this section we investigate anomalies of the I{ip) relation in several types of 
SFS structures which allow analytical solution: the SFcFS point contact with clean 
or diffusive constriction as a weak link, and the double-barrier SIFIS junction; the 
ferromagnetic layers are assumed to be thin, and the magnetization is homogeneous 
throughout the F part of the system. In particular, we show that the maximum 
of I{ip) can shift from 7r/2 ^ (fimax < tt to < v'max < 7r/2 as a function of the 
exchange field in the ferromagnet. Previously, current-phase relation of this type was 
theoretically predicted either if superconductivity in the S electrodes was suppressed 
by the supercurrent in the SNS structure jlHl EHI EDI or in the vicinity of T = in 
long SFS junctions |H1E1- 

The outHne of the section is as follows. We start with studying the SFcFS structure 
composed of two SF sandwiches linked by a clean Sharvin constriction with arbitrary 
transparency D. We show that the energy-phase relation of this junction can have 
two minima: aX ip — and Lp = tt (while the energy of the junction in the pure 0- 
or TT-state has a single minimum — a±Lp = OoTLp = Tr, respectively); such situation 
was discussed in Ref pi] for a clean SFS junction. As a result, the I{lp) dependence 
can intersect zero not only a± ip = and ip — tt but also at an arbitrary value po 
from the interval < tpo < tt. The salient effects which occur in junctions with 
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clean constriction survive averaging over the distribution of transmission eigenvalues 
and thus occur also in diffusive point contacts. Physically, the properties of SFS 
structures are explained by splitting of Andreev levels due to the exchange field; to 
demonstrate this, we study the spectral supercurrent (this physical mechanism was 
formulated for long SFS junctions at low temperatures in Refs. (HUHI)- Finally, we 
show that the same mechanism provides shifting of the I{(p) maximum to (^c < 7r/2 in 
the double-barrier SIFIS junctions which can be more easily reaHzed in experiment. 



3.2.2 SFcFS with clean constriction 

We start with a model structure (see Fig. 13.1b ) composed of two superconducting 
SF bilayers connected by a clean constriction with transparency D (the size of the 
constriction a is much smaller than the mean free path I: a I). We assume that 
the S layers are bulk and that the dirty Hmit conditions are fulfilled in the S and F 
metals. For simplicity we also assume that the parameters of the SF interfaces 7 and 
7b obey the condition 

7<max(l,7B). (3.20) 

We shall consider the symmetric structure and restrict ourselves to the limit of thin 
F layers: 

dp < min (^i., y/DF/2h}j . (3.21) 

Under condition (I3.2()|) . we can neglect the suppression of superconductivity in the 
S electrodes by the supercurrent and the proximity effect, and reduce the problem to 
solving the Usadel equation in the F layers: 



d 



2 

dx 



^'4x'^ 



_ ^GF<i>F - 0, (3.22) 



with the boundary conditions at the SF interfaces {x = Tdp)'- 

^,,k^JL^,.aJiL-il), (3.23) 

UJn OX \ OJn iOn J 

Gs^^4^=, <i>s(TdF) = Aoexp(Tz(p/2). 

Here Aq is the absolute value of the order parameter in the superconductors. 

Under condition II3.21II . the spatial gradients in the F layers arising due to the 
proximity effect and current are small. Then we can employ the solution H3.12|l . 
obtained in Sec. 13.11 

^Fi,F2 = $0 exp(T?(^/2), $0 = ^Ao, (3.24) 

where 



x/u}"^ + A^ dp 

W ^UJn+^nlBM^. ^ ^ iBM^lBT^, (3.25) 
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and the indices 1 and 2 refer to the left- and right-hand side of the constriction, 
respectively. 

The supercurrent in the constriction geometry is given by the general expression 

EH 

47rT {FiF2~ FiF2)/2 
1= Im > f 5—^^ -^Mi ----r, (3.26 

eRN ^^2-D[l-GiG2-{FiF2+FiF2)/2y ^ ' 

where Rn is the normal-state resistance of the junction. Inserting Eq. Ij3.24ll into 
this expression, we obtain 

_ 2'kT ^ A^sin(^ 

i^o^' + ^o [l-Z?sin^(^/2)]- ^^-'^^ 

Finally, the current-phase relation takes the form 

, , 27rr CA^sinoj 

C = A2 [1 - D sin^ (<^/2)] - {jBM^f + (1 + jBMnf , 

At small ojn, the function C [and hence /(</?)] changes its sign at finite phase 
difference 

o • /l - {iBMh/irTcY 
Lpc = 2 arcsm W (3.29) 

if the exchange field is in the range \ - D < {-iBMhl'KT,)'^ < 1. The results for 
I{lp) are shown in Figs. 13.61 13.71 and can be understood (similarly to Refs. 0121) by 
considering the spectral supercurrent density j{E) that is defined as the quantity 
which determines the supercurrent according to the formula 



I = 7nr- / dE f{E) b> (E) + ,n (E)] , (3.30) 

where f{E) — l/{e^/'^ + 1) is the Fermi distribution function. The physical meaning 
of j| (E) and ji (E) is the current carried by the quasiparticles with energy E and 
spin up or down (denoted as t or |) if the corresponding states are occupied, while 
f{E) describes their occupation.^ To calculate the spectral supercurrent density, we 
perform the analytical continuation in Eq. H3.27II : the result is 

/ „N An sin uy 
ME)^~lm 2 ^ (3.31) 

Al[l-Dsm\ip/2)] - E + -fBAiiE - ah)^/^^/TrT, 



^In the paper A. A. Golubov, M. Yu. Kupriyanov, Ya. V. Fominov, Pis'ma Zh. Eksp. Teor. Fiz. 
75, 709 (2002) [JETP Lett. 75, 588 (2002)] we used a different form of expression for the total 
current, which can be written in our present notations as 



1 /• oo f E \ 



4ei?jv 

However, here we prefer the expression in the form 13.301 because its physical meaning is clearer. 
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Figure 3.6: Current-phase relation in clean SFcFS junction with ideally transparent 
constriction (D = 1) at T/Tc — 0.01, ■jbm — 1 for different values of the exchange 
field h. Inset: spectral supercurrent density aX ip = 27r/3 for h/nTc = (solid line) 
and h/nTc = 0.4 (dotted line, the arrows near the peaks denote the corresponding 
spin direction). 



where ct = ±1 corresponds to spin up or down. Now we should recall that this formula 
originates from the retarded Green functions, hence we should add the infinitesimal 
imaginary part iO to the energy E. As a result, j{E) is given by a sum of delta- 
functions S{E — Ea) where Ea are the energies of the Andreev bound states. At 

7bm = 0, the well-known result Ea = ±Ao^J 1 — D sm'^{(p/2) is reproduced, while at 
finite 7sAf the exchange field splits each bound state into two (see inset in Fig. l3.6|l . 

At small temperature, the integration in Eq. II3.30|I is effectively restricted to the 
negative energy domain. The reversing of the supercurrent sign is then explained 
by the fact that at ip — ipc the positive peak for spin down crosses zero leaving the 
domain E < 0, and simultaneously the negative peak for spin up moves from the 
region E > Q into the region E < 0. 

The sign-reversal of the supercurrent (the O-tt transition) can also be achieved 
at a fixed h due to nonequilibrium population of levels. This phenomenon has been 
studied in long diffusive SNS l^^iij and SFS junctions [HIE]. 



3.2.3 SFcFS with diffusive constriction 

To get the I{(p) relation for the diffusive point contact [/ ^ a <C min(^i?, ^/Dpj2h)] 
we integrate Jg p{D)I{D)dD, where I{D) is given by Eq. (I3.27|l for the clean- 
constriction case (note that Rn oc in this equation) and p{D) is Dorokhov's 
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Figure 3.7: Current-phase relation in clean SFcFS junction at T/Tc = 0.01, jbm — 1, 
h/wTc = 0.8 for different values of the barrier transparency D. 



density function p{D) = 1/2D\/1 — D [21]. Finally, we arrive at the result 

X 47rr „ Ao cos(<^/2) / Ao sin(ip/2) \ 

= Rc > ' arctan ' . 

(3.32) 

This expression coincides with the direct solution of the Usadel equations (solving the 
Usadel equations, we should take into account that the gradients in the constriction 
region are large). In the 7bm = limit, the current does not depend on the charac- 
teristics of the F layers, and Eq. II3.32|I reproduces the Kulik-Omelyanchuk formula 
for the diffusive ScS junction j2Tl 122] . 

Calculation of I{}p) using the above expression yields results similar to those for 
the clean point contact, however the transition from 0- to 7r-state becomes less sharp 
(see Fig.EHl- 

Temperature dependence of the critical current in this case shows thermally- 
induced O-TT crossover with nonzero critical current at the transition point, in agree- 
ment with results of Refs. [2l El Ej- This is a natural result since the barrier 
transparency is high and the current-phase relation is strongly nonsinusoidal. 



3.2.4 SIFIS 

Now we turn to the double-barrier SIFIS junction (I denotes an insulating barrier) — 
see Fig. l3.1h . This structure is easier for experimental implementation than an SFcFS 
junction. 

We assume that condition H3.2flll is satisfied; then we can neglect the suppression 
of superconductivity in the S electrodes by the supercurrent and the proximity effect. 
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Figure 3.8: Current-phase relation in diffusive SFcFS point contact at T/Tc = 0.01, 
JBM = 1 for different values of the exchange field h. 



In this case the system is described by Eqs. (I3.22|1 - II3.23|I . although now instead of two 
F layers connected by a constriction we have a continuous F layer (at —dp < x < dp)- 
We also assume that the F layer is thin [condition (I3.21|l ] and that ^ dp/£_p, 
hence the spatial gradients in the F layer are small. Then we can expand the solution 
of Eqs. (|3.22|I - II3.2^^|I up to the second order in small gradients, arriving at 

/ t^riGs An sin((j9/2) x , . 

= $oCos(v?/2)+i^^^^ (3.33) 

Gp = =, (3.34) 

v/(^2+ci>2cos2(^/2) 

with <I>o defined in Eq. Ij3.24ll [in the final result H3.33|l we retained only the first 
order in gradients — this accuracy is sufficient for calculating the current]. 

Inserting the solution H3.33II . H3.34II into the general expression II3.10|I for the 
supercurrent, we obtain 

I(m) = ——Re > , ° ^ (3.35 

(our assumptions imply that i?jv ~ ^Rb)- This result demonstrates that the SIFIS 
junction with thin F layer is always in the 0-state.^ Nevertheless, I{ip) is strongly 
modified by finite h (see Fief. l3.9|l . especially at small temperatures. Figure clearly 
demonstrates that an increase of h results not only in suppression of the critical 



^In the case under discussion, when the F layers are thin and the interface parameters obey 
condition 13.201 . the phase of the order parameter is constant in the S part and almost constant in 
the F part, however it jumps at the two SF interfaces (see. Sec. 1,3.11 . The two jumps compensate 
each other in the SIFIS junction with a single F layer, whereas in the SFcFS junctions they add up 
at the weak link thus opening a possibility for the Tr-state. 



Josephson effect in SFS junctions 



95 




Figure 3.9: Current-phase relation in double-barrier SIFIS junction at T/Tc — 0.02, 
7bm = 1 for different values of the exchange field h. The value ip = 27r/3 will be 
used in Fig. Kui 



current, but also in the shift of the I{(p) maximum from (^max ~ 1-86 at ft. = to the 
values smaller than 7r/2. In the limit of large exchange fields, h/nTc » I/'^bm, the 
I{ip) dependence returns to the sinusoidal form. 

The physical origin of these results becomes clear in the real energy E representa- 
tion. Performing the analytical continuation in Eq. H3.35II by replacement a;„ — > —iE, 
we obtain the spectral supercurrent density: 

jAE) = - Im (3.36) 



yAf~PJA2cos2(^/2)-i?2 



E = E + -fBM{E-ah)y " 



nTr 



Equation Ij3.3(ill implies that at 



(^e = 2arccos(^;^^^ , (3.37) 

singularities in j{E) are shifted to the Fermi level. At (p > Lpc, the negative singularity 
in j| {E) crosses the Fermi level and appears in the E < domain (recall that at low 
temperatures only j{E) at E < contribute to the current), whereas the positive 
peak in ii(E) leaves the E <Q domain (this process is illustrated in Fig. KlKijl . As 
a result, the contribution to the supercurrent from low energies changes its sign, and 
the supercurrent I{ip) becomes suppressed at > (see Fig. I3.9|l . However, at 
energies with larger absolute values E ^ — Aq, modifications in i[E) are weak, and 
the resulting I[ip) does not change its sign. This mechanism is similar to the one that 
takes place in long SFS junctions at low temperatures [HI EI- 
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Figure 3.10: Spectral supercurrent density in diffusive double-barrier SIFIS junction 
with thin ferromagnetic interlayer at 7_bm = = 27r/3 for two values of the 

exchange field h. The chosen value of ip corresponds to Lpc at jBAih/nTc = 0.5, and 
the figure demonstrates that the positive peak for the spin down projection disappears 
from while the negative peak for the spin up projection appears in the E < domain. 
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3.2.5 Conclusions 

In this section, we have studied the current-phase relation in several types of Joseph- 
son junctions with thin ferromagnetic interlayers: SFcFS junction with ballistic or 
diffusive constriction and SIFIS planar junction. The current-phase relation can be 
highly nonsinusoidal in these junctions due to splitting of Andreev bound states by 
the exchange field. 

In particular, the maximum of /(</?) can be achieved at (p < 7r/2. In SFcFS 
junctions, the supercurrent can change its sign at arbitrary value of the phase between 
and n. When the maximum absolute value corresponds to negative current, the 
junction is in the 7r-state. 

Another consequence of this sign change is that the energy-phase relation for the 
junction has two minima: at = and ip = tt. The phenomenon can be used for 
engineering qubits. 
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Chapter 4 



Decoherence due to nodal 
quasiparticles in c?-wave Josephson 
junctions 

4.1 Introduction 

In addition to the fundamental importance, Josephson junctions between d-wave 
superconductors are of great interest for the rapidly developing field of quantum 
computations. In this field, the mathematical aspects (so to speak, writing programs) 
are rather well elaborated, while the quantum computer itself has not been realized 
so far, and there are only few successful realizations of a qubit (quantum bit). It 
is not clear at present, which of the proposed qubit implementation will prove most 
successful from the practical point of view. 

Solid-state proposals, and in particular superconducting devices, have a number 
of advantages, e.g., scalability and variability ll| (2j. Particularly interesting are 
the so-called quiet qubits, which are intrinsically degenerate, i.e., do not require 
any external source for maintaining the degeneracy. Such qubits can be realized 
in systems involving d-wave superconductors pl!^, which are of the so-called "phase 
qubit" type (the information is encoded by the phase difference (p across the junction). 
Recently, it was experimentally demonstrated that a double-well dependence of the 
energy versus the phase difference (inside a single period) is indeed realized in the 
Josephson junctions between d-wave superconductors [5^ (see Fig. I4.1|l . 

For any qubit implementation, there are processes that hamper its successful op- 
erating — the so-called decoherence processes, which destroy a quantum state of the 
qubit. Although the quiet phase qubits are rather well isolated from the environ- 
ment, there are intrinsic mechanisms of decoherence even at low temperatures. The 
quantum tunneling of the phase between the two minima leads to fluctuating voltage 
across the junction, which excites quasiparticles. The dissipative current across the 
interface arises, leading to a finite decoherence time r^p. 

The knowledge of is essential for estimating the efficiency of the qubit: short 
decoherence time makes the qubit senseless from the practical point of view, while a 
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Figure 4.1: Schematic dependence of the Josephson energy U on the phase difference 
ip. The barrier of the height 2Ej separates two nontrivial minima. The spHtting of 
the lowest energy level due to the tunnehng across the barrier is denoted At. 



long enough decoherence time opens the way for quantum correction algorithms that 
in principle allow to perform an infinitely long computation 0. 

The relevance of the quasiparticle processes at low temperatures is specific for d- 
wave superconductors. In the conventional s-wave case, the quasiparticle transport 
below the gap is suppressed. At the same time, in gapless anisotropic superconductors 
the gap vanishes in certain directions (the nodal directions), hence the low-energy 
quasiparticle appear. In the present chapter, we consider a DID Josephson junction 
(D = d-wave superconductor, I = insulator), and study the decoherence due to nodal 
quasiparticles (quasiparticles moving along the nodal directions) . 

4.2 Decoherence time (general strategy) 

Theoretical description of the quantum dynamics of a tunnel junction between two 
s-wave superconductors was developed in Ref. [3 (see Ref. jH] for a review); in par- 
ticular, the effective action for the phase difference ip was obtained. Later this de- 
scription was generalized to the case of rf-wave superconductors in Refs. PI EH- The 
effective action for ip is similar to the general case considered by Caldeira, Leggett et 
al. [n|E|, who studied influence of dissipation on quantum tunneling in macroscopic 
systems. The dissipation was described as being due to the interaction with a bath 
of oscillators (the environment). The "strength" of the environment, depending on 
the frequency ui, is characterized by the spectral function J{oj). In the Josephson 
junction, the environment is represented by the quasiparticles, and the spectral func- 
tion is given by hI{hLtj/e)/e, where / is the dissipative quasiparticle current taken at 
"voltage" huj/e f?,. 

A system Hving in a double- well potential and described by an extended coordinate 
can be "truncated" to the two-state system (spin 1/2) with the two states (cz = ±1) 
corresponding to the minima of the potential (see Fig. l4.1|l . The theory of dissipative 
two-state systems is thoroughly elaborated ^2] for the cases when the spectral func- 
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tion behaves as J(w) oc up to some high-frequency cutoff. The situations when 
s = 1, s > 1, and < s < 1 are called ohmic, superohmic, and subohmic, respectively. 
In this language, the dissipation due to nodal quasiparticles in the Josephson junction 
is superohmic, as we demonstrate below. 

What is the decoherence in such a system? Assume that during the time t < the 
system is held in the right well (i.e., at = 1). At t = the constraint is released, 
and we consider the expectation value of the system coordinate: P{t) — {<7z{t)). 
Below we shall encounter the superohmic case at zero temperature. Then ^21 

P{t) = cos{Att/h) exp(-t/T^) (4.1) 

— the cosine describes coherent oscillations between the two wells (At is the tun- 
nel splitting of levels, see Fig. I4.1|l while the exponential leads to their incoherent 
damping. 

The decoherence time is expressed in terms of the spectral function |^ . Re- 
turning from the general theory to the particular case of the Josephson junction, we 
write the corresponding result as 



4e Airh 



(4.2) 



where Sip is the distance between the potential minima and Rq — nh/e'^ ~ IBkil 
is the quantum resistance. Comparing the decoherence time with the characteristic 
time of oscillations between the wells, h/At, we obtain the quality factor 

Q = I^^ (43) 

which must be large for successful operating of the qubit. 

In the DID junction, the tunnel splitting At is much smaller than the order pa- 
rameter A, hence T^p is determined by the quasiparticle current at low "voltage". 



4.3 Quasiparticle current 

We consider the grain-boundary Josephson junction between two quasi-two- 
dimensional dx2-y2-waYe superconductors with cyHndrical Fermi surfaces. The ori- 
entations of the superconductors are characterized by the angles between the a-axes 
and the normal to the interface (the x axis) — see Fig. 14.21 According to Ref. 
we consider the mirror junction, in which the misorientation angles on both sides are 
equal in magnitude but opposite in sign, a/ — a (we take —45° ^ a < 45° because 
all physically different situations in the mirror junction are reahzed in this interval). 
The order parameter depends on the direction (parametrized by the angle 9) and the 
distance to the interface: 

Al,«(x, 0) = AL,fl(x)e*'^^'« cos {20 T 2a) , (4.4) 

where the indices L and R refer to the left- and right-hand side of the junction, 
respectively. 
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The two generate states of the DID junction correspond to two minima of the 
junction energy and carry no current across the interface (the Josephson current is 
proportional to the derivative of the Josephson energy with respect to the phase; the 
derivative is zero at the minima). At the same time, physically, the two states differ 
by the direction of spontaneous currents flowing along the interface |18l I14j . These 
currents produce fluxes that can interact with the environment and make the d-wave 
qubit not completely quiet (21 ■ From this point of view, the mirror orientation of 
the junction is the most preferable one, because in the mirror junction there is no 
total current along the interface (the currents are nonzero but compensated). In 
principal, the presence of the spontaneous currents implies the phase change along 
the interface. However, the currents appear as a result of the Andreev reflections, 
which are of the second order over transparency. In the tunneling limit, the currents 
and the corresponding phase change along the interface are small, and A has the form 
l|4.4ll with constant ipL and (pR. 

The quasiparticle current in the tunneling limit at low temperatures, /cgT ^ huj, 
is given by 

I(hu/e) = ^^ i^(g)j;os0 1^ dEv{E- hiu, 0) v {E, 9) , (4.5) 



ei?jv J-tt/2 D 

£» = / deD{e) COS0. 

J-tt/2 

Here Rn is the normal-state resistance of the interface, iy{E, 6) is the density of states 
(DoS) at the interface, normalized to the normal-metal value, and D(9) is the angle- 
dependent transparency of the interface. We have not labelled the DoS by the indices 
L and R because vl[E,9) — vr{E,9) in the mirror junction. 

Below we calculate the nodal contribution to the current 114. 5|l at ftw ^ Aq, where 
Ao — A(±oo) is the bulk amplitude of the order parameter. The angle integration 
contributing to the current is then Hmited to narrow angles around the nodal direc- 
tions, where the low-energy DoS is nonzero (as we shall see below, the width of the 
angles is 69 = hxj/Ao). 

To calculate the DoS, we employ the quasiclassical approach. The quasiclassical 
matrix Green function 

a.{i _g ,«) 

obeys the Eilenberger equation |15[ll6j and satisfies the normalization condition = 
1. It can be parametrized as 

_ 1-ab _ 2a - _ 2b 

^"l + afc' -'^l + ab' ■'^l + ab' ^ ' 

then the normalization condition is automatically satisfied. The equations for the 
new functions a{x, 6) and b{x, 9) take the form of the Riccati equations [Tt] : 

hvF COS 6* ^ - 2iEa + A*a^ - A 0, 
dx 

hvF cos9^ + 2iEb - Ab^ + A* = 0, (4.8) 
dx 
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Figure 4.2: DID junction of mirror orientation a/ — a. The positive lobes of the order 
parameter are shaded. An electron e moving along a truly nodal direction 9i of the left 
superconductor, tunnels into an induced nodal direction of the right superconductor. 
^R{Si) 7^ 0, therefore the electron experiences the Andreev reflection; the hole h 
returns to the interface, and after reflection at the interface escapes into the bulk 
along the truly nodal direction —9i. In this process, the total current into the bulk of 
the right superconductor is composed of the Cooper pair along 6i and the hole along 
-9i. 

where v-p is the absolute value of the Fermi velocity vp, and 9 denotes the angle 
between vp and the x axis. 

In the tunneling limit, the DoS is calculated at an impenetrable interface. Let 
us consider, e.g., the right superconductor (the right half-space). We need to find 
the low-energy DoS in two cases: 1) in the vicinity of a nodal direction, so that 
E,A{9) <C Ao, 2) at a gapped direction, so that E <C A(6'). In the first case, the 
spatial scale £,e = ?iWF cos0/|£+| on which the quasiclassical Green functions vary 

(we denote £± = \J — |A(oo, ±6')|^), is much larger than the coherence length 
^ = hvp/2nkBTc on which variations of A occur. This allows us to regard A as 
constant when integrating Eqs. H4.8|l over x. In other words, the functions a and b 
at low energies do not feel the suppression of A near the interface, because it takes 
place on a small scale. In the second case, the spatially dependent parts of a and b 
are proportional to E/A{9) <C 1 and hence small. Thus a and b at the interface are 
equal to their bulk values, as if A was constant. 

Thus we can regard A{x,9) as equal to the bulk value Ao{9) — A{oo,9). The 
integration of the functions a and b over x in Eqs. H4.8II is stable only in the directions 
determined by the sign of 008 6*. At cos9 > 0, the function b{x,9) is stably integrated 
from a; = oo to the interface {x = 0), hence 

6(0,0)^6(00,0) = *:^^^. (4.9) 

At the same time at cos 9 > 0, the function a is stably integrated from the interface 
to a; = oo. Therefore to find a(0, 9), we consider the trajectory directed along n — 9. 
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Since cos(7r — 9) < 0, the function a is stably integrated from a; = oo to the interface. 
Finally, the direction n — 6 is converted to 6 upon reflection at the specular interface: 

a(0, 9) = a(0, n ~ 9) ^ a(oo, n - 9) ^ i ^~^7'f"^ . (4.10) 
As a result, the DoS = Reg at the interface is 

The gap in the spectrum is ^^(e') = min (|Ao(6')| , |Ao(-6')|). 

The DoS is symmetric, i^{9) = i'{—9), because the Green functions are continuous 
upon reflection. Thus in each superconductor there are two "truly" nodal directions 
9i [i = 1, 2) in the interval — 7r/2 < 9 < 7r/2, and also two "induced" nodal directions 
—9i. Near a nodal direction Eg{9) — 2Ao|6' — 9i\. Along a truly nodal direction, the 
gap vanishes and the DoS is the same as in the normal metal, v{E) = 1. For an 
"induced" nodal direction this is so only near the interface. 

In the left superconductor, the truly nodal directions are 6*1,2 — a± 45°. Due to 
the mirror symmetry, the truly nodal directions of the right superconductor coincide 
with the induced nodal directions of the left one, and vice versa. In total, there 
are four nodal directions in the junction, which are symmetric with respect to the 
interface normal. 

In this situation, the transport is due to the processes of the following type. An 
electron moving along a truly nodal direction 9i of the left superconductor, tunnels 
into an induced nodal direction of the right superconductor (see Fig. I4.2|l . However, 
the electron cannot escape into the bulk of the right superconductor because Aji{9i) ^ 
0. Therefore the electron experiences the Andreev reflection ^Sl; the hole returns to 
the interface, and after reflection at the interface escapes into the bulk along the 
truly nodal direction —9i. In this process, the total current into the bulk of the right 
superconductor is composed of the Cooper pair along 9i and the hole along —9i, which 
is overall equivalent to the transfer of one electron. 

The nodal contribution to the current Il4.5|l appears only due to integrating in the 
vicinity of the nodal directions where Eg < hw. The DoS near the nodal directions at 
small energies can be found from Eq. II4.11|I . Below we distinguish the general case 
when Aq{9) ^ ±Ao{—9), and two special cases: Ao{9) = Ao{—9) (at a = 0°) and 
Ao{9) = -Aoi-9) (at a = 45°). 

At a = 0°, the truly nodal and induced nodal directions coincide in each super- 
conductor, and Eq. H4.11|l yields the BCS-like DoS: 

\E\ 

t/Qo (g, 61) = Rc ' ' . (4.12) 

^E^~\Ao{9)\' 

At Of = 45°, the truly nodal and induced nodal directions again coincide, and Eq. 
Ij4.11ll yields the DoS of the inverse BCS type: 



iy^,.iE,9)=Re^ — . (4.13) 



Decoherence due to nodal quasiparticles in d-wave Josephson junctions 



107 



2.5 








A„(-e) - A„(e) 


2.0 








1 Ao(-6) - 3 Ao(e) 

A„(-6) = 8 A„(e) 


1.5 




N 




V 1 










I \ 


> 1.0 












^ \ ■ 
\'. 
\' 




I / 

7 


0.5 




\ 




[' 


0.0 











-10 -5 5 10 



Figure 4.3: Density of states following from Eq. H4.11|l . The energy is normalized to 
Ao{9), while Ao{—9) is varied. 



Finally, if \a\ » huj/Ao and 45° — \a\ » hw/Ao (i.e., a is not too close to 0° and 
±45°), then Ao{9) in the essential angle of the width S9 — hui/Ao around a nodal 
direction is much smaller than Ao{—9). Then in the region of energies that contribute 
to the quasiparticle current, |Ao(0)| < \E\ < hw <ti \Ao{—9)\, the DoS is again given 
by the inverse BCS formula: 



^Je^-\Ao{9)\ 



2 



iygiE,9^9,)=Re^ — . (4.14) 



Fisfure IT^ demonstrates the DoS at different angles 9, which are parametrized 
by different ratios Ao(-6l)/Ao(6i). At A o{-9) = Ao(6'), the DoS has the BCS-like 
square-root singularity near Eg [see Eq. 14.1211 ]. At Aq{—9) ^ Ao{9), the DoS has 
the inverse-BCS behavior near Eg [see Eq. H4.14II ]. 

Inserting Eqs. H4.12ll - H4.14|l into Eq. H4.5|l . we obtain: 

I{hu; e) = -^r^^^^ > ~ , (4.15) 

eRN Ao ,ft2 ^ 

where 6*1^2 = Q;±45° and ^ is a number, which depends on the orientation of crystals: 
A{0°) « '0.46, ^(45°) « 0.19, and A{a) = 2A{45°) « 0.37 when a is not too close to 
0° or ±45°. 

In Refs. (niE|, the quadratic current-voltage characteristic, / oc lj^, was obtained 
for the case of aligned nodal directions (i.e., for the a/ a orientation). 



108 



Chapter 4 



4.4 Estimate 



Equations l|4.3|l . H4.15II yield: 



Q 



( 



5(p^A{a) J2 D{di)cosei/D J Rq At' 



27r \ Rn Ao 




(4.16) 



1=1,2 



To proceed further, we need to estimate the tunnel splitting At (see Fig. l4.1|l . This 
is easily done due to the formal equivalence between quantum mechanics of a Joseph- 
son junction described by the phase (p and a particle described by the coordinate x 
[19] . To employ this analogy, we should substitute 



where m is the particle mass and Ec = e'^/2C is the charging energy of the Josephson 
junction (C is its capacitance). For the estimate, we assume that the second harmonic 
dominates in the energy-phase relation, U{tp) = Ej{l + co%2ip), and the energy of the 
levels is small compared to Ej. Then the attempt frequency is wq = 'iy/TEjEc/fi, 
the tunnehng action is calculated between the points ip — —n/2 and 7r/2, and we 
obtain: 



To obtain a numerical estimate, we take the characteristics of the junction as in 
the experiment of Il'ichev et al. 0. The capacitance of the junction is C ~ 10~^^ F 
|2(1] . hence Ec/ks ~ O.IK. The characteristic Josephson energy is of the order of 
several Kelvin. For estimate, we take 2Ej/kB = 7K. The resistance of the interface 

is i?7v ~ 50O [sg. 

As a result, At/ks ~ 2.5 • 10"'' K. Finally, we estimate S(p ~ tt, Ao/fc^ ^ 200 K, 
and assume a thin (5-functional barrier with transparency D{9) — Dqcos'^ 9, then the 
quality factor is Q ~ lO'^-^ 10*. Here we have retained only the order of magnitude for 
Q, because we cannot expect a higher accuracy in the case when important character- 
istics of the junction (e.g., C and Ej) are known only by the order of magnitude. We 
also made an essential assumption that the second Josephson harmonic dominates, 
which is the most favorable case, hence the obtained estimate is "optimistic". 

The latter assumption can be realized under special conditions, while in a more 
common situation the first and the second harmonics are of the same order. Es- 
timates for this case were done in a recent work [2l], where the characteristics of 
mesoscopic junctions between high-Tc superconductors were experimentally studied 
and theoretically analyzed. A characteristic value of At ^ 0.1 K was reported under 
the conditions that correspond to Rn 100 Jl. Assuming such parameters for the 
mirror junction, we obtain Q ~ 10-^ 10^, which is a more realistic estimate than that 
for the experiment of Ref . ^ . 

The above estimates for Q are very different. At the same time, a general con- 
sequence of Eq. II4.16|I is that the quality factor grows as the splitting At becomes 




(4.18) 
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smaller. We note in this respect, that the values of the critical current (and hence the 
Josephson energy) measured in Refs. (HE], are much smaller than expected. If the 
critical current is enhanced to the expected value, then At decreases, which finally 
leads to an increase of Q. 

If a 7^ 0°, the low-energy quasiparticles are presented not only by the nodal 
quasiparticles, but also by the midgap states (MGS) with zero energy (22]. In the 
case of specular interface and clean superconductors, considered in this chapter, the 
DoS corresponding to the MGS is proportional to 5[E), hence the MGS on the two 
sides of the interface do not overlap and do not contribute to the current at a finite 
voltage. 

In the asymmetric case, when ul ^ ±Q!_r (precisely speaking, when | |q;l| — |q;_r| | > 
hLo/Ao), the nodal directions of the left and right superconductors do not match each 
other. At first sight, this solves the problem of decoherence since the transport 
from nodal to nodal direction is suppressed. However, a more important transport 
"channel" arises: between the nodal directions and the MGS. This leads to a stronger 
decoherence than in the symmetric case. 

In the mirror junction, the MGS contribute to the quasiparticle current if they 
are split and/or broadened [23- • To take into account the contribution of the MGS 
into decoherence, the present approach should be considerably modified. This issue 
requires a separate study. 

The low-energy quasiparticles (and dissipation due to them) can be suppressed 
due to finite size of the d-wave superconductors. This issue is schematically discussed 
in Sec. 14.51 below. 



4.5 Finite size effects 

The effect of finite size on the low-energy quasiparticles in DID junctions can be most 
easily studied in the quasi-one-dimensional geometry assuming finite thickness L (the 
dimension in the x direction) of the rf-wave superconductors. Then each trajectory is 
successively reflected from the I interface and from the outer surface of the layer. This 
leads to periodicity of the order parameter profile along a trajectory (see Fig. 14.411 . 
The finite size leads to two main effects: 

1. For nodal quasiparticles (this corresponds to A(— 0) = in Fig. 14.411 . the Andreev 
quantization arises, similarly to the SNS junction |18ll24j . The system of equidistant 
energy levels appears, and the level spacing is ttvfCOsO/L pHI, while the lowest 
level is iTVpcos9/2L [21]. Thus, the energy gap is formed. 

2. The MGS, which exist if A(6') and A(-6') have different signs, split. This effect, 
similar to the level splitting in a double-well quantum potential, is due to the 
overlap of the localized states across the lowest of the barriers (A(— 0) in Fig. l4.4|l . 
To estimate the corresponding splitting, below we consider a simplified model shown 
in Fig. 14.51 
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Figure 4.4: Schematic illustration of a periodic order parameter profile along a trajec- 
tory in a rf-wave layer of width L. We have chosen the trajectory with the incidence 
angle 0, hence the part of the trajectory between two successive reflections has the 
length L/ cos6. Actually, at a 7^ 0° the order parameter is suppressed near the in- 
terfaces, however, to capture the qualitative effects we neglect this suppression and 
consider the model of piecewise constant A. 



4.5.1 Tunnel splitting of the zero Andreev level 

A classical quantum-mechanical problem is splitting of the level in a symmetric 
double-well potential. A similar effect takes place in the case of two coupled su- 
perconducting point contacts. Let us consider the splitting of the zero-energy levels 
in two coupled point contacts with the phase difference of tt across each contact (see 
Fig.|131l. 

We shall use the technique of the Bogolyubov - de Gennes equation: 

In the Andreev approximation we separate rapid oscillations of the solution: 

= (^2) ^M^kFx\ (4.20) 

and retain only the first derivatives of the slow functions uq and vq. The Bogolyubov 
- de Gennes equation then takes the form 




d. . . (4.21) 



The boundary conditions are the continuity of uq and vq at the interfaces. The 
sign of the momentum is a good quantum number, and we can first consider the 



Decoherence due to nodal quasiparticles in d-wave Josephson junctions 



111 



A2 




A2 




-Ai 





Figure 4.5: The profile of the order parameter, corresponding to two point contacts. 

anzatz Il4.2fl|l . which corresponds to the momentum directed to the right. In the 
general case of arbitrary phase difference, the energy levels for the momentum directed 
to the left can then be obtained after inverting the sign of the phase. At the same 
time, in the case of real A (when the phase difference is either or tt), the energy 
levels are degenerate with respect to the momentum direction. 
The solutions in the three spatial regions are 



voj + VA?-^V \E- i^Al - 



= D[^ , a<z, (4.22) 



uo _ 

voJ ^ \E- i^A'i - E^ 



where 



/Al,~E^ 

91,2 - . (4.23) 

Substituting these solutions into the boundary conditions at the two interfaces (a; = 
and X = a), we obtain four linear equations for the four coefficients A, B, C, and D. 
The spectrum is determined from the condition that the corresponding determinant 
is zero. The straightforward algebra leads to the equation 



{E' + A1A2) tanh(gia) - ^(A^ - E^){Al - E^) = 0. (4.24) 

In the limit aAi/wp ^ 1, the splitting is small, and the levels are 

2A1A2 / aAi\ 
S = ±— — — exp . (4.25) 



Ai + A2 \ vf 

4.6 Conclusions 



In this chapter, we have proposed an approach that allows to calculate the decoher- 
ence time due to nodal quasiparticles in the DID junctions, which can be used as 
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phase qubits. The dissipation in the junctions of mirror orientation is weaker than in 
the asymmetric ones. We find the superohmic dissipation with s = 2 in the mirror 
junction, which becomes weak at small splitting of the ground state. The superohmic 
case is most favorable (compared to ohmic and subohmic) for possible qubit applica- 
tions. For available experimental data, wc estimate the quality factor. Finally, wc 
illustrate the possibility to suppress quasiparticle decoherence due to finite size of the 
d-wave superconductors. 
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Appendix A 

Interpretation of r as escape times 



The quantities ts and tn introduced in Eq. H1.31|l may be interpreted as escape times 
from the corresponding layers. The arguments go as follows. 

If the layers are thin, then the diffusion inside the layers is "fast" and the escape 
time from a layer is determined by the interface resistance. The time of penetration 
through a layer or the interface is determined by the corresponding resistance: Rs{n) 
or Rb, hence the diffusion is "fast" if Rs(n) ^ Rb- 

Let us use the detailed balance approach, and consider an interval of energy dE. 
In the S layer, the charge in this interval is Qs = evs dE Ads- Let us define the 
escape time from the S layer ts, so that the current from S to N is equal to Qs/ts- 
On the other hand, this current can be written as dE/eRB, hence 



and we immediately obtain 



Q^_dE_ 



ts = {A.2) 
E>S 



Similarly, we obtain the expression for the escape time from the F layer t^r. As a 
result, the relations between the quantities t defined in Eq. H1.31|l and the escape 
times t are simply 

TS = 2ts, TN = 2tN. (A.3) 



115 



Appendix B 

Analytical results for SF bilayers 



B.l Limit of thin S layer 

(i) When ds ^ £,s and h >• nTcs, problem H2.ir)|l - H2.13ll can be solved analytically. 
The first of the above conditions implies that A can be considered constant, and F"*" 
weakly depends on the spatial coordinate; so F~^{x, w„) — 2A/uJn + A{uJn) cosh{ks[x — 
ds]). The boundary conditions determine the coefficient A; as a result 



2A 



AsiiOn) + W{iOn) 



(B.l) 



where ks, As, and W are defined in Eq. (12. 1211 . Finally, the self-consistency equation 
for Tc takes the form 



In 



ReV' 



1 , 7^s 



1 



2 2dsjB+Bi, T, 



where Bh does not depend on uJn due to the condition h ^ ttTcs- 

kh 



B, = 



kh£,F ta,nh{khdF) 



1 



cS 



In the limiting case of perfect interface (75 =0) and thick F layer {dF 
ljB.2ll reproduces the result of Ref. jT|. 



(B.2) 

(B.3) 
00), Eq. 



(ii) If the F layer is also thin, dF ^ \/DF/2h, Eq. IIB.2p is further simplified: 



In 



cS 



Re?A 



1^ 
2 Ts 



1 



-i + TFh 



2^T, 



where ts, tf are defined similarly to Eq. Ijl.31|l : 



TS 



2dsRBA 
PsDs 



TF 



2dFRBA 
PfDf 



(B.4) 



(B.5) 
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and have the physical meaning of the escape time from the corresponding layer. They 
are related to the quantities 7, 75 used in the body of the paper as 

IB 1 ds I dp , . 

rs = T^-T' '^P = "fB^^-^- (B.6) 

(iii) If the S layer is thin, ds <C ■^s, and the SF interface is opaque, jb ^ 00, the 
critical temperature of the bilayer only slightly deviates from Tcs- In this limit Eq. 
I|B.1|I applies with W — 7/75 ^ 1, and we finally obtain: 

r. = T.5-^. (B.7) 

Interestingly, characteristics of the F layer {dp, h, etc.) do not enter the formula. 
In particular, this formula is valid for an SN bilayer (2] (where N is a nonmagnetic 
normal material, h — 0) because Eq. ljB.7|l was obtained without any assumptions 
about the value of the exchange energy. 



B.2 Cooper limit 



In the Cooper limit, when both layers are very thin [ds <C ^jDsJ^^^) dp ^ 
min(-\/£'F/2w_D, \jDp/2h), with lod the Debye energy of the S metal] and the in- 
terface is transparent, the bilayer is equivalent to a homogeneous superconducting 
layer with internal exchange field. This layer is described by effective parameters: 
the order parameter A*^''*^^ the exchange field h'^'^^\ and the pairing constant A'-°'^^. 
In this subsection we develop the ideas of Ref. [HI, demonstrate a simple derivation 
of this description, and find the limits of its applicability. 

The Usadel equations 112.111 . 112. 2|l for the two layers can be written as a single 



equation: 

Dp9{-x) +DsO(x) d'^F 
2 



\uJn\F - ihsgn{un)0{~x)F + A0{x) = 0, (B., 



where 9 is the Heaviside function [9{x > 0) = 1, 9{x < 0) = 0]. The self-consistency 
equation (I2.3|l can be rewritten as 

A(a;) = \9{x)ttT F{x, c^„), {B.9) 

where A is the pairing constant. 

First, we consider the ideal SF interface: 75 = [see Eq. (I2.48|l ]. then F{x) is 
continuous at the interface and nearly constant across the whole bilayer, i.e., Fs{x) ~ 
Fp{x) = F. Applying the integral operator to Eq. IIB.8|l : 

^ dx + —r^ — — I ' dx (B.IO) 



i^sds + vpdp 7_dp vsds + vpdp Jo 
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(here ly is the normal-metal density of states), and cancelling gradient terms due to 
the boundary condition II2.47|I . we obtain the equations describing a homogeneous 
layer: 

-\un\Fiun) - ih^''''^ sgn(w„)FK) + A^^*^) = 0, (B.ll) 
^(cff) ^ ^(off)^y J2 FK), (B.12) 

with the effective parameters (see also Ref. 0): 

/j(eff) ^ -^^h, A(<=ff) = ^^A, (B.13) 
\(off) . ^(cff) exp(C) / 1 

where C is Euler's constant and Tj:^^^ is the critical temperature of the layer in the 
absence of ferromagnetism (i.e., at /i^''^^ — 0). The critical temperature is determined 
by the equation 

ln^=Re^ - - . (B.14) 



Tc \2 2ttTcJ V2, 

Actually, the description in terms of effective parameters l|B.13p is appHcable at 
an arbitrary temperature (i.e., when the Usadel equations are nonlinear) and has 
a clear physical interpretation: the superconducting (A, A) and ferromagnetic {h) 
parameters are renormalized according to the part of time spent by quasiparticles 
in the corresponding layer. This physical picture is based on interpretation of r as 
escape times, which we present in the next subsection. 

Now we discuss the applicability of the above description for a nonideal interface 
(7b 7^ 0). In this case F is nearly constant in each layer, but these constants are 

'2 
S 



different: Fs{x) « Fs + Csix-ds f , Ff{x) « Fp + Cpix+dF? , where \Fs\ > \Cs\d 



and \Ff\ S> ICfI^I-. Using the Usadel equation ljB.811 and the boundary conditions 
1233, (E3H|), we find the difference SF = Fg ~ Fp: 

A 

SF = p -. (B.15) 



1 

Ts 



1 

1 



Tp {\LUn\ + ihsgnUJn)_ 

Finally, the homogeneous description is valid when \SF/F\ <^ 1 [with F determined 
by Eq. ljB.ll|l ]. which yields: 

niax(ft,, cj£)) max(T5, Ti?) <C 1 (B.16) 

(here w„ ~ w^i has been taken as the largest characteristic energy scale in the quasi- 
homogeneous bilayer). 
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Appendix C 

Applicability of the single-mode 
approximation for calculating Tc in SF 
bilayers 



C.l Case study 

As pointed out in Sec. 12.1.31 the single-mode approximation (SMA) is applicable only 
if the parameters of the bilayer are such that W [see Eq. Il2.12|l ] can be considered 
w„-independent. An example is the case when 7^ ^ \Bh\, hence W = 7/75. 

The condition ^ \Bh\ can be written in a simpler form; to this end we should 
estimate \Bh\. We introduce the real and imaginary parts of khi kh — k'/^ + ik'j^, and 
note that /c^ > /c^'. Then using the properties of the trigonometric functions and the 
estimate tanha; ~ min(l, x) we obtain 

\Bh\^[k',^^FtSinh{k'„dF)r\ (C.l) 
and finally cast the condition 75 ^ \Bh\ into the form 



-L « n.in I ^max ^) ; | max -A- ) } , (C.2) 

where the ratio Tc/Tcs originates from (jJu/t^TcS with w„ ~ ttTc as the characteristic 
energy scale in the bilayer. 

If condition ljC.2|l is satisfied, then the SMA is vaHd and Tc is determined by the 
equations 

Otanfl7^^=^. (C.4) 

\ is J IB 

These equations can be further simplified in two Hmiting cases which we consider 
below. 
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IB is 

in this case Eq. (IC.4II yields Vt^ ^ and Eq. ljC.3ll takes the form 



7-B ds 



In 



which reproduces the 75 ^ | Hmit of Eq. IIB.211 . 

(2)^^»1: 

IB is 



■4, 



in this case Eq. \CAl yields ri||- = |, and Eq. ll(L3ll takes the form 

TaS 



In- 



2 

^2 r 



Tcs_ 



(C.5) 



(C.6) 



Equations IIC.3ll - ljC.6|l can be used for calculating the critical temperature Tc and 
the critical thickness of the S layer dcs below which the superconductivity in the SF 
bilayer vanishes (i.e., Tj. = 0). 



C.2 Results for the critical temperature 

In the limit when Tc is close to Tcs, Eqs. ll(L5ll . iKLfill yield 



Te = I 1 - 



Tf_i_is_ 

4 7_B ds 



it — <C mm 



IB 



is ' ds 



and 



T, = T 



cS 



1 - 



4 ds 



•f '^■S ^ /, IB 

II — » max 1, — 

is V 7 



(C.7) 



(C.8) 



Using relations IB.fijl one can check that result l|(L7|l is equivalent to Eq. 



C.3 Results for the critical thickness 

The critical thickness of the S layer dcs is defined as the thickness below which there 
is no superconductivity in the SF bilayer: Tc{dcs) — 0.^ When Tc — *■ 0, Eq. ljC.3ll 
yields — l/-y/2exp(C) (where C ~ 0.577 is Euler's constant), and Eq. lj(".4|l takes 
the form 

, ' tan ( , ' ^] = ^. (C.9) 
v/2cxp(C) \ v/2exp(C) is J 1b 

ExpHcit results for dcs can be obtained in limiting cases: 

— = 2exp(C) — if — < 1, 

£S 7-B 7-B 4S 



(C.IO) 



^When calculating dcS, we assume that the phase transition is of the second order. However, this 
issue may require a separate study. In principle, under some circumstances, the order of the phase 
transition in SF systems can change from second to first one — see Sec. 12.1.51 
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and 



dcS 



exp(C) 



7 ds 
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Spatial dependence of the order 
parameter in SF bilayers 



According to the self-consistency equation, in the S layer the proximity order param- 
eter F{x, r = 0) is proportional to A(x): 



where A is the pairing constant which can be expressed via the Debye energy: 



A(a;) can be found as the eigenvector of the matrix L — 1 ln(T'c5/Tc) [see Eq. II2.33II ]. 
corresponding to the zero eigenvalue. 

After that we can express F{x, t = 0) in the F layer via A(x) in the superconduc- 
tor. The Green function Fpix^ujn) in the F layer is given by Eq. 112. 7|l . with C(ujn) 
found from the boundary conditions: 



A{x) 



(D.l) 




(D.2) 




(D.3) 



The Green function at the S side of the SF interface is 



Fs{0,i^n) 



F+{0,UJn)+Fg{0,UJ^) 

2 



(D.4) 



The symmetric part F, 



is given by Eq. II2.31|I . The antisymmetric part is 



C (ujn) cosh {ks[x 



ds]), 



(D.5) 



with C (oJn) found from the boundary conditions: 



^7 Im Bh 1 f/(0,a^n) 



(D.6) 



AsllB + Bh\'^+ -f{-fB + ReBh) J cosh(fcsds) ' 
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Finally, the proximity order parameter in the F layer is the Fourier transform [see 
Eq. jinU] of 



FF{x,UJn) 

Bh 



1 



i"f Im Bh 



AsllB + Bh\^ + -f{-fB + ReBh) 
cosh (^kh[x + dp] 



7b + Bh J cosh{khdF) 



GiO,y;u„)A{y)dy. 



(D.7) 



Summary 



The research described in this thesis deals with thermodynamic and transport prop- 
erties of various superconducting junctions. The systems under consideration are 
composed of conventional s-wave superconductors in contact with normal metal or 
ferromagnet, and of anisotropic rf-wave superconductors.^ Various aspects of the 
proximity and .Joscphson effects in such systems are studied theoretically. 

When a superconductor contacts a normal metal, a number of phenomena known 
as the proximity effect takes place. The two materials influence each other on a 
spatial scale of the order of the coherence length in the vicinity of the interface. In 
particular, the superconducting correlations between quasiparticles are induced into 
the normal metal, because the Cooper pairs penetrating into the normal metal have a 
finite lifetime there. Until they decay into two independent electrons, they preserve 
the superconducting properties. Alternatively, the proximity effect can be viewed as 
resulting from the fundamental process known as the Andreev reflection, at which an 
electron impinging from the normal metal onto the interface with the superconductor 
is reflected back as a hole, while a Cooper pair is transferred into the superconductor. 

The Josephson effect is a macroscopic quantum effect that occurs when two super- 
conductors are connected via a weak link, which can be either a nonsuperconducting 
material (insulator or metal) or a geometrical constriction. The superconducting con- 
densate in each of the two weakly coupled superconductors is described by its wave 
function and the corresponding phase. The essence of the Joscphson effect is the 
appearance of the supercurrent (nondissipative current) between the superconductors 
at finite phase difference ^p between them. The supercurrent can flow in the absence 
of voltage, and in the simplest case is given by the sinusoidal .Josephson relation 
/ = /cSin(/?. The quantity Ic (or, more generally, the maximal supercurrent) is called 
the critical current. 

Superconductivity in thin SN bilayers 

In Chapter 1, the proximity effect in thin SN bilayers is studied. 

Although the investigation of the proximity effect in SN systems was started 
about forty years ago, the technology allowing to produce and measure experimental 
samples of mesoscopic dimensions was achieved relatively recently. In particular, it 
became possible to study SN structures consisting of thin layers (having thickness 

^The notations: S — s-wave superconductor, D — d-wave superconductor, N — normal metal, F 
— ferromagnetic metal, I — insulator, c — constriction. 
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smaller than the coherence length). Such structures behave as a single superconduc- 
tor with nontrivial properties. From practical point of view, the proximity structures 
can be used as superconductors with relatively easily adjustable parameters, in par- 
ticular, the energy gap and the critical temperature. The parameters of the proximity 
structures can be tuned, e.g., by varying the thicknesses of the layers. This method 
has already found its apphcation in superconducting transition edge bolometers and 
photon detectors for astrophysics. 

While most of theoretical works on SN proximity structures have focused on the 
hmit of ideally transparent interfaces, the experimental progress requires the advances 
in theory that take into account arbitrary interface transparency. This crucial pa- 
rameter determines the strength of the proximity effect and at the same time is not 
directly measurable. 

In Chapter 1, the theory of superconductivity in thin SN sandwiches (bilayers) 
in the diffusive hmit is developed, with particular emphasis on the case of very thin 
superconductive layers, dg ^ (In ■ The proximity effect in the system is governed by 
the SN interface resistance (per channel) Pb- The energy gap, the order parameter, 
the critical temperature, the density of superconducting electrons, the parallel and 
perpendicular critical magnetic fields of the bilayer are investigated as functions of 
pB- 

The case of relatively low resistance (which can still have large absolute values) can 
be completely studied analytically. This case corresponds to the situation when the 
diffusion of quasiparticles between the layers is fast compared to the characteristic 
time A^^. The theory describing the bilayer in this limit is of Bardeen-Cooper- 
Schrieffer (BCS) type but with the minigap (in the single-particle density of states) 
E'g <C A substituting the order parameter A in the standard BCS relations; the orig- 
inal relations are thus severely violated. In the opposite limit of an opaque interface, 
the behavior of the system is in many respects close to the BCS predictions. Over the 
entire range of pB, the properties of the bilayer are found numerically. It is demon- 
strated that the gap is a nonmonotonic function of the interface resistance, reaching 
its maximum in the region of moderate resistances. Also in the region of moderate re- 
sistances, a jump of the parallel critical field (due to a redistribution of supercurrents 
in the bilayer) is discovered. Finally, it is shown that, due to geometrical symmetry, 
the results obtained for the bilayer also apply to more comphcated structures such as 
SNS and NSN trilayers, SNINS and NSISN systems, and SN superlattices. 

Proximity effect in SF systems 

In Chapter 2, the critical temperature in SF bilayers and FSF trilayers is studied. 

Compared to SN structures, the physics of SF systems is even richer. In contrast 
to the SN case, the superconducting order parameter does not simply decay into the 
nonsuperconducting metal but also oscillates. This behavior is due to the exchange 
field in the ferromagnet that acts as a potential of different signs for two electrons in 
a Cooper pair and leads to a finite momentum of the pair (similarly to the Larkin- 
Ovchinnikov-Fulde-Ferreh state in bulk materials) . This oscillations reveal itself in 
nonmonotonic dependence of the critical temperature of SF systems as a function 
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of the F layers thickness. At the same time, in most of the works investigating this 
effect, the methods to calculate Tc were approximate. An exact method to calculate 
Tc at arbitrary parameters of the system was lacking. The need for such a method 
was also motivated by a recent experiment that did not correspond to the previously 
considered approximations and limiting cases. 

In Chapter 2, numerical methods are developed to exactly calculate the critical 
temperature of a dirty SF bilayer at arbitrary parameters of the structure (thick- 
nesses of the layers ds and dp, interface transparency). The methods arc applied to 
study the nonmonotonic behavior of the critical temperature versus thickness of the 
F layer. Good agreement with experimental data is demonstrated. In hmiting cases, 
analytical results for the critical temperature and the critical thickness of the S layer 
are obtained. 

Another interesting effect in SF systems takes place if the magnetization of the 
ferromagnet is inhomogeneous. Then the triplet superconducting component can arise 
in the system. The triplet component corresponds to pairing between electrons with 
the same spin projection (while in the conventional case the Cooper pairs are formed 
by electrons with opposite spin projections). Recently, it was demonstrated that the 
triplet component also arises in the case of several homogeneous but noncollinearly 
oriented ferromagnets. However, the conditions at which the superconductivity is 
not destroyed in this system were not found. The simplest system of the above 
type is an FSF trilayer. The answer to the question about the conditions for the 
superconductivity to exist can be obtained when studying the critical temperature 
of the system. At the same time, a method to calculate Tc in a situation when the 
triplet component is generated, was lacking. 

Therefore, the methods developed for SF bilayers in the first part of Chapter 2, are 
then applied to study the critical temperature of the dirty FSF trilayer at arbitrary 
parameters of the system (mutual orientation of the magnetizations a, thicknesses of 
the layers, interface transparency). Changing of characteristic types of the T^dp) 
dependence while varying the mutual orientation is demonstrated. In the general case, 
it is analytically proven that Tc{q) is a monotonic function. In interesting limiting 
cases, analytical results for the critical temperature of the system and the critical 
thickness of the S layer are obtained. The necessary conditions for the existence 
of the odd (in energy) triplet superconductivity in multilayered SF structures are 
formulated. 

A possible practical application of FSF structures is a spin valve, a system that 
switches between superconducting and nonsuperconducting states when the relative 
orientation of the magnetizations is varied. Although the superconductive spin valve 
is not yet experimentahy reahzed, the work in this direction has already started. 

Josephson effect in SFS junctions 

The Josephson effect in the systems containing ferromagnets (e.g., structures of the 
SFIFS, SIFIS or SFcFS type) also has a number of pecuharities. Among them are 
the transition from the ordinary (0-state) to the so-cahed 7r-state (in other words, the 
inversion of the critical current sign or the additional tt phase shift in the Joseph- 
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son relation) when the two ferromagnets are ahgned in parallel, the enhancement 
of the critical current by the exchange field in the case of antiparallel orientation, 
and nonsinusoidal current-phase relation. The interest to SFS junctions with non- 
trivial current-phase relation is in particular due to their possible employment for 
engineering logic circuits of novel types (both classical and quantum bits). 

In Chapter 3, the Josephson effect in dirty SFIFS junctions is studied by a self- 
consistent method at arbitrary system parameters (F layers thickness, interface trans- 
parencies) . When the magnetizations of the ferromagnets are antiparallel the effect of 
the critical current enhancement by the exchange field is observed, while in the case 
of parahel magnetizations the junction exhibits the transition to the 7r-state. In the 
limit of thin F layers, these peculiarities of the critical current are studied analytically 
and explained qualitatively. The mechanism of the O-tt transition, at which the phase 
jumps by 7r/2 at the two SF interfaces, is discovered. It is demonstrated that the 
effect of the critical current enhancement by the exchange field (at the antiparallel 
magnetizations ahgnment) is similar to the Riedel singularity in SIS junctions. The 
logarithmic divergence of the maximal critical current is cut off by finite temperature 
or finite transparency of the tunnel interface. The effect of switching between and 
TT states by changing the mutual orientation of the F layers is demonstrated. 

Nonsinusoidal dependence of the Josephson current on the phase difference in 
SFcFS point contacts and planar double-barrier SIFIS junctions is studied in the 
hmit of thin F layers. It is demonstrated that in SFcFS junctions, the current can 
cross zero not only at = and = tt but also at an intermediate value. This 
implies that the energy of the junction has two minima — at Lp = Q and = tt, hence 
the and tt states of the junction coexist. If the minima are of the same depth, the 
system can be used as a quantum bit (qubit). The physical mechanisms leading to 
highly nontrivial /((^) dependence are identified by studying the spectral supercurrent 
density. 

Decoherence due to nodal quasiparticles in d-wave 
Josephson junctions 

Another interesting type of nonuniform superconducting systems is a junction between 
superconductors of nontrivial symmetry. The anisotropic superconductors with the 
d-wave symmetry of the order parameter are widely discussed because this symmetry 
is realized in the high-temperature superconductors. A possibility to implement a 
so-called qubit (quantum bit) based on d-wave junctions was proposed theoretically. 
Quantum bit is, simply speaking, a quantum mechanical system with two states 
(which can be imagined as spin 1/2). While a classical bit can be either in one state 
or in the other, a qubit can also be in a superposition of the two states. If a quantum 
computer is built of such qubits, it would have the advantage of natural computational 
parallelism that can enormously speed up certain types of computational tasks. A 
possibility to implement a qubit based on DID junction stems from the fact that the 
energy of such a junction as a function of the phase difference can have a double- 
well form with two minima. This degeneracy of the ground state arises due to the 
nontrivial symmetry of the superconductors. Due to tunneling between the wells. 
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the ground state splits, and the two resulting levels effectively form the quantum- 
mechanical two-state system. At the same time, the gapless nature of the (i-wave 
superconductors leads to appearance of low-energy quasiparticles which can destroy 
the quantum coherence of the qubit and hence hamper its successful functioning. 
Calculation of the corresponding decoherence time is necessary for estimating the 
efficiency of the proposed qubits. 

In Chapter 4, the decoherence time due to the low-energy nodal quasiparticles 
in the Joscphson junction between two rf-wave superconductors is calculated. The 
decoherence is due to an intrinsic dissipative process: quantum tunneling between the 
two minima of the double-well potential excites nodal quasiparticles which lead to in- 
coherent damping of quantum oscillations. In DID junctions of the mirror orientation 
(a/ — a), the contribution to the dissipation from the nodal quasiparticles is found 
to be superohmic, hence the quality factor (the number of coherent oscillations that 
occur before the decoherence overcomes) grows as the tunnel splitting of the ground 
state decreases. The quality factor is estimated for available experimental data on 
DID junction. The suppression of low-energy quasiparticles due to finite size of the 
d-wave superconductors is discussed. 
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De onderzoeksthema's van dit proefschrift zijn de thermodynamische eigenschappen 
en transporteigenschappen van verscheidene supergeleidende contacten. De 
onderzochte systemen zijn samengesteld uit conventionele s-type supergeleiders 
in contact met een normaal metaal of ferromagneet en uit anisotrope rf-type 
supergeleiders.^ Een aantal aspecten van het Josephson effect en het effect van de 
nabijheid van een supergeleider in dergelijke systemen is theoretisch onderzocht. 

Als een supergeleider in contact wordt gebracht met een normaal metaal vindt 
een aantal verschijnselen plaats dat onder de noemer nabijheidseffect valt. De 
twee materialen beinvloeden elkaar op het grensvlak over een lengteschaal van 
de ordegrootte van de coherentielengte. De supergeleidende correlaties tussen 
quasideeltjes worden in het normale metaal gei'nduceerd vanwege de eindige 
verblijftijd van de binnendringende Cooperparen. Tot het moment van uiteenvallen 
in twee onafhankelijke elektronen blijven de supergeleidende eigenschappen bewaard. 
Aan de andere kant kan het nabijheidseffect worden gezien als een manifestatie van 
het fundamentele proces dat Andreev reflectie wordt genoemd, waarbij een elektron 
vanuit het normale metaal op het grensvlak met de supergeleider wordt gereflecteerd 
als een gat terwijl een Cooperpaar wordt doorgegeven aan de supergeleider. 

Het Josephson effect is een macroscopisch kwantumeffect dat optreedt wanneer 
twee supergeleiders verbonden zijn via een zwakke koppehng, hetgeen een niet- 
supergeleidend materiaal (isolerend of metallisch) kan zijn of een geometrische 
constrictie. Het supergeleidende condensaat in elk van de twee supergeleiders 
wordt door een golffunctie beschreven met elke een fase. De essentie van het 
Josephson effect is het optreden van een superstroom (niet-dissiperende stroom) bij 
een eindig faseverschil if tussen de supergeleiders. Superstroom kan zonder aangelegd 
potentiaalverschil stromen en is in het meest eenvoudige geval gegeven door de 
sinusoi'dale Josephson relatie / = /c sin 99. De grootheid Ic (in het algemeen de 
maximale superstroom) wordt de kritieke stroom genoemd. 

Supergeleiding in dunne SN-bilagen 

In hoofdstuk 1 is het nabijheidseffect in dunne SN-bilagen bestudeerd. 

Hoewel het onderzoek naar het nabijheidseffect in SN-structuren reeds veertig jaar 
geleden is begonnen, is de technologie om experimented structuren met mesoscopische 

^De notaties zijn: S — *-type supergeleider, D — d-type supergeleider, N — normaal metaal, F 
— ferromagnetisch metaal, I — isolator, c — constrictie. 
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afmetingen te prepareren en te meten pas recent bereikt. In het bijzonder werd het 
mogelijk SN-structuren te bestuderen bestaande uit dunne lagen (met een dikte die 
kleiner is dan de coherentielengte) . Dergelijke structuren gedragen zich als een enkele 
supergeleider met niet-triviale eigenschappen. De SN-structuren kunnen, praktisch 
gczicn, wordcn gcbruikt als supcrgclcidcrs met rclaticf makkclijk aan tc passcn 
eigenschappen zoals de 'gap' en de kritieke temperatuur. De eigenschappen kunnen 
bijvoorbeeld worden aangepast door de laagdiktes te varieren. Deze methode wordt 
rccds gcbruikt in bolometers die zijn gcbaseerd op dc ovcrgang van de supergeleidende 
naar normale toestand en fotondetectoren voor astrofysische doeleinden. 

Hoewel de meeste theoretische studies naar SN-structuren zich hebben 
geconcentreerd op de limiet van volledig transparante grensvlakken, vereist de 
experimented ontwikkeling een vooruitgang in de theorie voor een willekeurige 
transparantie van het grensvlak. Deze cruciale parameter bepaalt de sterkte van 
het nabijheidseffect maar is niet direct meetbaar. 

De theorie van de supergeleiding in dunne SN-structuren (bilagen) in de limiet van 
diffuus transport is ontwikkeld in hoofdstuk 1, met een specifieke nadruk op structuren 
met een hele dunne supergeleidende laag, ds ^ (In- Het nabijheidseffect in het 
systeem wordt bepaald door de weerstand van het SN-grensvlak (per kanaal) pB- De 
'gap', de ordeparameter, de kritieke temperatuur, de dichtheid van supergeleidende 
elektronen en de parallelle en loodrechte kritieke magneetvelden van de bilaag zijn 
onderzoclit als functie van Pb- 

Het geval van een relatief lage weerstand (hetgeen nog steeds absoluut hoge 
waardes toelaat) kan geheel analytisch worden bestudeerd. Dit geval komt overeen 
met de situatie waarin de diffusie tussen quasideeltjes snel is vergeleken met de 
karakteristieke tijd A~^. De theorie die de bilaag in deze limiet beschrijft is de 
Bardeen-Cooper-Schrieffer (BCS) theorie waarin nu de ordeparameter A vervangen 
wordt door de 'minigap' (in de toestandsdichtheid van de quasideeltjes) Eg <^ A, 
hetgeen dus een ernstige schending van de oorspronkelijke vergelijkingen veroorzaakt. 
Het systeem in de tegenovergestelde limiet van een beperkt transparant grensvlak 
gedraagt zich echter in veel opzichten ongeveer naar de voorspellingen uit de BCS 
theorie. De eigenschappen voor een bilaag zijn numeriek bepaald voor alle Pb- 
Het is afgeleid dat de 'gap' een niet-monotone functie van de weerstand van het 
grensvlak is, waarbij het maximum wordt bereikt in het gebied met middelmatige 
weerstandswaarden. In dit gebied met middelmatige weerstandswaarden is een sprong 
ontdekt in het parallelle kritieke veld (vanwege een herverdeling van superstromen 
in de bilaag). Tenslotte is aangetoond dat, vanwege geometrische symmetrie, de 
resultaten voor een bilaag ook gelden voor ingewikkeldere structuren zoals SNS- en 
NSN-trilagen, SNINS- en NSISN-structuren en SN-superstructuren. 

Het nabijheidseffect in SF-structuren 

In hoofdstuk 2 is de kritieke temperatuur van SF bilagen en FSF-trilagen bestudeerd. 

Vergeleken met de SN-structuren is de fysica van SF-structuren zelfs nog 
gevarieerder. In tegensteUing tot in SN-structuren, dringt de supergeleidende 
ordeparameter niet slechts eenvoudigweg het normale metaal binnen, maar kan 
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daarbij ook oscilleren. Dit gedrag wordt veroorzaakt door het 'exchange' veld in 
de ferromagneet, hetgeen zich gedraagt alsof er voor elk van de twee elektronen 
in een Cooperpaar een potentiaal met tegengestelde polariteit bestaat, waardoor 
het paar een impuls krijgt (vergelijkbaar met de Larkin-Ovchinnikov-Fulde-Ferrell 
tocstand in materialcn). Dczc oscillatics manifcstcrcn ziclizclf in hot nict-monotonc 
verband tussen de kritieke temperatuur Tc van SF-systemen en de laagdikte van de 
F-laag. In de meeste studies naar dit effect zijn de methoden om Tc uit te rekenen 
slcclits bcnadcringcn cn ccn prcciczc mcthodc om Tc uit tc rckcncn als functic van dc 
parameters van het systeem ontbrak. De behoefte aan een dergelijke methode kwam 
ook voort uit recente experimented bevindingen die niet overeenkomen met de tot 
nu toe beschouwde limietgevallen en benaderingen. 

In hoofdstuk 2 zijn numerieke methodes ontwikkeld om de kritieke temperatuur 
van een SF-bilaag in de diffuse limiet exact uit te rekenen als functie van de 
parameters van het systeem (laagdiktes ds en dp, grcnsvlak transparantie). De 
methoden zijn toegepast in een studie naar het niet-monotone verband tussen de 
kritieke temperatuur en de laagdikte van de F-laag. Een goede overeenstemming met 
de experinientele gegevens is hierljij Ijereikt. In speciale limietgevallen zijn analytische 
uitdrukkingen gevonden voor de kritieke temperatuur en de kritieke dikte van de S- 
laag. 

Een ander interessant effect vindt plaats in SN-structuren wanneer de magnetisatie 
van de ferromagneet niet homogeen is. In dit geval kan een triplet component in de 
supergeleiding ontstaan. De triplet component komt overeen met de paring tussen 
elektronen met dezelfde projectie van de spin (terwijl in het conventionele geval de 
Cooperparen worden gevormd door elektronen met tegengestelde spin projecties). 
Het is recent aangetoond dat de triplet component ook kan ontstaan in systemen 
met homogcne, maar niet-colineair georienteerde ferromagneten. De voorwaarden 
waaronder de supergeleiding niet wordt tenietgedaan in zulke structuren waren echter 
nog niet gevonden. Het meest eenvoudige voorbeeld van een dergelijke structuur is 
de FSF-trilaag. Het antwoord op de vraag naar de voorwaarden voor het bestaan van 
supergeleiding kan worden beantwoord door de kritieke temperatuur van het systeem 
te onderzoeken. Een methode om Tc te berekenen in structuren waarbij een triplet 
component wordt gegenereerd, ontbrak echter nog. 

Vanwege dit feit, zijn de methodes voor de SF-bilaag uit hoofdstuk 2 gebruikt 
om de kritieke temperatuur van een diffuse FSF-trilaag te bepalen als functie van 
de parameters van het systeem (wederzijdse orientatie van de magnetisaties a, 
laagdiktes en grensvlak transparantie). Het is aangetoond hoe de functies Tddp) 
veranderen wanneer de orientaties worden gewijzigd. In het meest algemene geval is 
het analytisch bewezen dat Tc{a) een niet-monotone functie is. In de interessante 
hmietgevallen zijn analytische uitdrukkingen voor de kritieke temperatuur en de 
kritieke dikte van laag S verkregen. De noodzakelijke voorwaarden voor het 
bestaan van oneven (in energie) triplet supergeleiding in multilaags SF-structuren 
is geformuleerd. 

Een mogelijke praktische toepassing van FSF-structuren is de spin-transistor, een 
systeem dat schakelt tussen de supergeleidende en de niet-supergeleidende toestand 
wanneer de relatieve orientatie van de magnetisaties wordt veranderd. Hoewel de 
supergeleidende spin transistor experimenteel nog niet is gerealiseerd, is het onderzoek 
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in deze richting wel begonnen. 

Het Josephson effect in SFS contacten 

Het Josephson effect in structuren die ferromagneten bevatten (bijvoorbeeld SFIFS-, 
SIFIS-, SFcFS-structuren) heeft sen aantal bijzondere eigenschappen, waaronder 
de overgang van de gewone toestand (0-toestand) naar de zogenaamde 7r-toestand 
(met andere woorden: het teken van de kritieke stroom wordt geinverteerd, of: een 
additionele faseverschuiving van tt vindt plaats) wanneer de ferromagneten parallel 
gcoricntccrd zijn, dc tocnamc van dc kritiokc stroom door hot 'exchange' veld wanneer 
de ferromagneten anti-parallel georienteerd zijn en een niet-sinusoi'dale stroom-fase 
relatie. De interesse in SFS-contacten met een niet-triviale stroom-fase relatie 
komt met name voort uit de mogelijke toepasbaarheid in het realiseren van nieuwe 
elektronische circuits (met klassieke of kwantum bits). 

In hoofdstuk 3 is het Josephson effect bestudeerd in diffuse SFIFS-structuren 
met willokourig gekozen parameters (laagdiktc van F, grcnsvlak transparantie) door 
middel van een consistente methode. Wanneer de magnetisaties van de ferromagneten 
anti-parahel zijn georienteerd, wordt een toename van de kritieke stroom door het 
'exchange' veld gevonden, terwijl in het geval van een parallelle orientatie het contact 
een overgang naar de 7r-toestand kan hebben. Deze bijzonderheden van de kritieke 
stroom zijn analytisch bestudeerd en kwalitatief verklaard in de limiet waarin de 
laagdikte van F klein is. Het mcchanisme van dc overgang van naar tt, waarbij de 
fase met 7r/2 verspringt op de twee SF-grensvlakken, is ontdekt. Het is aangetoond 
dat de toename van de kritieke stroom door het 'exchange' veld (bij een anti- 
parallelle orientatie) gelijkenis vertoont met de Riedel singulariteit in SIS-contacten. 
De logaritmische divergentie van de maximale kritieke stroom wordt afgezwakt door 
de eindige temperatuur of de niet-ideale transparantie van het grensvlak. Het 
verscliijnsel van schakelen tussen de 0- en de vr-toestand door het veranderen van 
de wederzijdse orientatie van de ferromagnetische lagen is aangetoond. 

Het niet-sinusoidale verband tussen de Josephson stroom en het faseverschil in 
SFcFS-puntcontacten en planaire dubbele-barriere SIFIS-contacten is bestudeerd in 
het limietgeval waarbij de F laag dun is. Het is aangetoond dat de stroom in SFcFS- 
contacten niet alleen de nulwaarde kan passeren bij = en = tt, maar ook 
bij een tussenliggende waarde. Dit betekent dat dc cncrgic van het contact twee 
minima heeft (op = en = tt) en de 0- en de 7r-toestand tegelijk bestaan. 
Wanneer de minima een gelijke diepte hebben, kan het systeem worden gebruikt 
als een kwantumbit (qubit). Het fysische meclianisme dat voor de niet-triviale I{lp) 
relatie zorgt, is ontrafeld door het bestuderen van de spectrale superstroomdichtheid. 

Decoherentie door quasideeltjes in de nodale richting 
van d-type Josephson contacten 

Een ander interessant type niet-uniform supergeleidend systeem is een contact 
tussen supergeleiders met een niet-triviale symmetrie. Veelal worden anisotrope 
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supergeleiders bestudeerd met de symmetrie van de ordeparameter van het li-type 
omdat deze symmetrie wordt gerealiseerd in hoge-temperatuur supergeleiders. De 
mogelijkheid om een zogenaamd qubit (kwantumbit) gebaseerd op d-type contacten 
te implementeren is reeds theoretisch voorgesteld. Een qubit is in feite een 
kwantummechanisch systccm met twee toestandcn (dit kan wordcn voorgesteld 
als spin 1/2). Terwijl een klassieke bit in de ene of in de andere toestand is, 
kan een qubit zich ook in een superpositie van de twee toestanden bevinden. 
Een kwantumcomputcr gebaseerd op qubits zou dan hot voordccl hcbbcn dat een 
natuurlijke vorm van parallel rekenen kan worden gebruikt, wat bepaalde soorten 
numerieke taken enorm zou kunnen versnellen. De mogelijkheid een qubit te 
implementeren met behulp van DID-contacten komt voort uit het feit dat de energie 
van dit contact een dubbele potentiaalput met twee minima vormt als functie van 
het faseverschil. Deze gedegenereerde grondtoestand ontstaat vanwege de niet-triviale 
symmetrie van de supergeleiders. De grondtoestand splitst vanwege tunnelproccsscn 
tussen de energieminima en de twee resulterende energietoestanden vormen een 
kwantumechanisch twee-toestandensysteem. Het ontbreken van de 'gap' in d-type 
supergeleiders leidt tot de aanwezigheid van laag-energetische quasideeltjes die de 
kwantumcoherentie van de qubit kunnen onderdrukken en daarmee een succesvol 
gebruik bemoeilijken. 

In hoofdstuk 4 is de tijd van de dccolicrentic door de laag-energetische 
quasideeltjes uitgerekend voor Josephson contacten met twee d-type supergeleiders. 
De decoherentie wordt veroorzaakt door een intrinsiek dissipatief proces: een 
kwantumtunnelproces tussen de twee minima van de dubbele potentiaalput slaat 
nodale quasideeltjes aan, wat aanleiding geeft tot een incoherente demping van de 
kwantumoscillaties. De bijdrage aan de dissipatie van de nodale quasideeltjes in DID 
contacten met een gespiegeldc orientatic (q/ — q) is super-ohmisch. Hicrdoor stijgt de 
kwaliteitsfactor (het aantal coherente oscillatie dat plaatsvindt voordat decoherentie 
optreedt) bij een afnemend verschil tussen de grondtoestanden. De kwaliteitsfactor 
is gescliat voor de beschikbare experimented gegevens. De onderdrukking van laag- 
energetische quasideeltjes door een beperkte laagdikte van de d-type supergeleiders 
is toegelicht. 
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